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Abstract 



We analyse the effects of thermal quasiparticles in leptogenesis using hard-thermal-loop-resummed 
propagators in the imaginary time formalism of thermal field theory. We perform our analysis in 
a leptogenesis toy model with three right-handed heavy neutrinos A'^i, A'^2 and N^. We consider 
decays and inverse decays and work in the hierarchical limit where the mass of A'^2 is assumed to 
be much larger than the mass of A'^i , that is M2 ^ Mi . We neglect flavour effects and assume that 
the temperatures are much smaller than M2 and M3. We pay special attention to the influence of 
fermionic quasiparticles. We allow for the leptons to be either decoupled from each other, except 
for the interactions with neutrinos, or to be in chemical equilibrium by some strong interaction, 
for example via gauge bosons. In two additional cases, we approximate the full hard-thermal-loop 
lepton propagators with zero-temperature propagators, where we replace the zero-temperature 
mass by the thermal mass of the leptons mi{T) in one case and the asymptotic mass of the 
positive- helicity mode ^/2mi{T) in the other case. We calculate all relevant decay rates and CP- 
asymmetries and solve the corresponding Boltzmann equations we derived. We compare the final 
lepton asymmetry of the four thermal cases and the vacuum case for three different initial neutrino 
abundances; zero, thermal and dominant abundance. The final asymmetries of the thermal cases 
differ considerably from the vacuum case and from each other in the weak washout regime for zero 
abundance and in the intermediate regime for dominant abundance. In the strong washout regime, 
where no infiuences from thermal corrections are commonly expected, the final lepton asymmetry 
can be enhanced by a factor of two by hiding part of the lepton asymmetry in the quasi-sterile 
minus-mode in the case of strongly interacting lepton modes. 



Zusammenfassung 



Wir analysieren die Effekte von thermischen Quasiteilchen in Leptogenese, wobei wir Propagatoren 
verwenden, die durch harte thermische Schleifen resummiert sind. Wir arbeiten ini Imaginar- 
zeit-Formalismus der thermischen Feldtheorie. Unsere Analyse wird in einem Beispielmodeh von 
Leptogenese niit drei rechtshandigen Neutrinos A^i, A'^2 und N^ durchgefiihrt. Wir betrachten 
Zerfahe und inverse Zerfalle und nehmen den hierarchischen Grenzfah an, in dem die A'^2-Masse 
wesentlich groBer als die A^i -Masse ist, das heif5t M2 ^ Mi. Wir vernachlassigen Flavoureffekte 
und nehmen an, dass die Temperaturen sehr viel kleiner sind als M2 und M3. Wir legen beson- 
deres Augenmerk auf den Einfluss der fermionischen Quasiteilchen und lassen sowohl eine vollige 
Entkopplung der Leptonmoden bis auf die Neutrinowechselwirkung zu, als auch eine starke Kop- 
plung der Moden, beispielsweise durch Eichbosonen. In zwei zusatzlichen Fallen nahern wir den 
vollstandigen, durch harte thermische Schleifen resummierten Leptonpropagator durch normale 
Vakuumpropagatoren an, die eine Masse erhalten, die in einem Fall der thermischen Leptonmasse 
mi{T) und in einem anderen Fall der asymptotischen Masse der positiven Helizitatsmode, \/2m^(r) 
entspricht. Wir berechnen alle relevanten Zerfallsraten und CP- Asymmetrien und losen die 
hergeleiteten Boltzmann-Gleichungen. Wir vergleichen die resultierende Leptonasymmetrie der vier 
thermischen Falle und des Vakuumfalls fiir drei verschiedene Anfangswerte der Neutrinoverteilung; 
verschwindende, thermische und dominante Anfangsverteilung. Die Leptonasymmetrien der ther- 
mischen Falle unterscheiden sich in gewissen Parameterregionen stark vom Vakuumszenario als 
auch untereinander; namlich im schwachen Washout-Regime fiir verschwindende Anfangsverteilung 
und im intermediaren Regime fiir dominante Anfangsverteilung. Aufierdem vergrofiert sich die 
finale Leptonasymmetrie im starken Washout-Regime, wo typischerweise keine thermischen Effekte 
erwartet werden, um einen Faktor von etwa zwei, wenn ein Teil der Leptonasymmetrie auf die 
quasisterile Minus-Mode im Falle stark wechselwirkender Leptonmoden iibertragen wird. 



Introduction 

And the angel of the presence spake [...]: 
Write the complete history of the creation, 
how in six days the Lord God finished all 
His works and all that He created [...]. [5] 

Leptogenesis 2, 1 

The question of the origin of ah things that we observe and that are present in our life and 
surroundings, which in its last consequence is nothing else than the question of the origin of mankind 
itself, has always fascinated us and driven us to search for answers in science, religion, philosophy 
and the arts. On the scientific side, physics as the study of the laws of nature (f] cpuaixr] "nature") 
and within physics, cosmology (6 xba\i.oQ "order") as the science of the order and the evolution of 
the universe, address this question and have their own formulation of it. What is the origin of the 
matter that is the building block of all things we observe, including ourselves? 

The matter in nature consists of electrons, which belong to the leptons (Xstitoc; "small"), and 
the much heavier protons and neutrons, which belong to the baryons (papuc; "heavy") and are in 
turn made up of quarks. The theory that describes how the smallest ingredients of matter, the 
elementary particles, interact, is the standard model of particle physics (SM), which has been tested 
to great accuracy by experiment. According to the SM, matter particles, quarks or leptons, can 
only be created in pairs together with their antiparticles, that is, antiquarks and antileptons. If 
we assume that the early universe was indeed without form and void^ [6] , that is in the language 
of particle physics, there was no excess of one particle species over the other, there would have 
to be an equal amount of particles and antiparticles today. More specifically, since annihilation 
of particles and antiparticles proceeds at fast rates, no structures like atoms, molecules, galaxies, 
stars, planets, DNA, cells and finally living organisms could have formed and we would observe^ a 
universe populated almost exclusively by photons and the slowly interacting neutrinos. There are 
only two possible ways out of this obviously wrong scenario: Either we assume an excess of matter 
over antimatter as an initial condition of the big bang or we find some mechanism which does not 
strictly obey the conservation of baryon and lepton number and can create such an asymmetry 
dynamically in some early phase of the evolution of the universe. 



^In fact, according to standard cosmology, the early universe was by no means void, but rather a vibrant soup of 
all particles that we know and possibly many more species interacting rapidly with each other. 
^Or rather, not observe. 
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It might seem tempting to assume an excess of particles over antiparticles as an initial condi- 
tion of the universe, or more specifically the excess of baryons over antibaryons, since this is the 
asymmetry we measure on cosmological scales. The fact that this approach has to be refused is 
not even mainly because it would be a highly unsatisfactory approach from a scientific point of 
view, but the assumption itself clashes with another important theory in early universe cosmology, 
inflation. There is broad consensus that, in order to cure serious problems of cosmology, the early 
universe must have undergone such a phase of rapid expansion, which is so fast that it would dilute 
any baryon asymmetry we could realistically impose as an initial condition of the young universe. 
We are thus bound to find a mechanism that creates a baryon asymmetry dynamically, we need a 
theory for baryogenesis. 

Among the many baryogenesis theories that solve the problem of the matter-antimatter asym- 
metry, we focus on a variant called leptogenesis [7], which is a particularly attractive model since 
it simultaneously solves two problems: The creation of a baryon asymmetry via the detour of a 
lepton asymmetry on the one hand, and the explanation of why the neutrinos have such a small 
mass compared to all the other particles of the SM via the seesaw mechanism [8-11] on the other 
hand. In short, one adds heavy, right-handed neutrinos to the SM, which interact with the SM 
neutrinos and suppress their mass. In the early universe, these heavy neutrinos decay into leptons 
and Higgs bosons and create a lepton asymmetry, which is lateron converted to a baryon asymmetry 
by anomalous SM processes [12,13]. 

Ever since the development of the theory 25 years ago, the calculations of leptogenesis dynamics 
have become more refined and many effects and scenarios that have initially been neglected have 
been considered^. Notably the question how the hot and dense medium of SM particles influences 
leptogenesis dynamics has received increasing attention over the last years [15-22]. At high tem- 
perature, particles show a different behaviour than in vacuum due to their interaction with the 
medium: they acquire thermal masses, modified dispersion relations and modified helicity proper- 
ties. All these properties can be summed up by viewing the particles as thermal quasiparticles with 
different behaviour than their zero-temperature counterparts, much like the large zoo of single- 
particle and collective excitations that are known in high density situations in solid-state physics. 
At high temperature, notably fermions can occur in two distinct states with a positive or negative 
ratio of helicity over chirality and different dispersion relations than at zero temperature. 

Thermal effects have been considered by references [15-22]. Notably reference [16] performs an 
extensive analysis of the effects of thermal masses that arise by resumming propagators using the 
hard thermal loop (HTL) resummation within thermal field theory (TFT). However, the authors 
approximated the two fermionic helicity modes with one simplified mode that behaves like a vacuum 
particle with its zero-temperature mass replaced by a thermal mass^. Due to their chiral nature, 



■^For an excellent review of the development in this field, we refer to reference [14]. 

^Moreover, an incorrect thermal factor for the CP-asymmetry was obtained, as has been pointed out in refer- 
ence [23]. 
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there are serious consequences to assigning a chirality breaking mass to fermions, hence the effects 
of abandoning this property should be examined. Moreover, it seems questionable to completely 
neglect the negative-helicity fermionic state which, according to TFT, will be populated at high 
temperature. We argue in this study that one should include the effect of the fermionic quasipar- 
ticles in leptogenesis calculations and possibly in other early universe dynamics, since they behave 
differently from zero-temperature states with thermal masses, both conceptually and regarding 
their numerical influence on the final lepton asymmetry. We do this by analysing the dynamics 
of a leptogenesis toy model that includes only decays and inverse decays of neutrinos and Higgs 
bosons, but takes into account all HTL corrections to the leptons and Higgs bosons, paying special 
attention to the two fermionic quasiparticles. In a slightly different scenario, we assume chemical 
equilibrium among the two leptonic modes, thereby simulating a scenario where the modes interact 
very fast. As a comparison, we calculate the dynamics for two models where we approximate the 
lepton modes with ordinary zero-temperature states and modified masses, the thermal mass m^{T) 
and the asymptotic mass of the positive-helicity mode, \/2m£{T). 

The thesis is structured as follows: In chapter 1, we present a short overview over leptogenesis 
and explain the standard dynamics. We also discuss the limitations of our approach, where we 
do not include flavour and resonant effects or effects from a possible abundance of N2 or N^. 
Chapter 2 is devoted to a brief and comprehensive introduction into TFT in the framework of the 
imaginary time formalism (ITF), where we present the necessary ingredients for the further analysis, 
in particular frequency sums for fermions and bosons. We pay special attention to the resummation 
of hard thermal loops, which form the basis for the description of thermal quasiparticles. In 
chapter 3, we present the toy model of leptogenesis and calculate decays and inverse decays. At 
high temperature, the thermal mass of the Higgs bosons becomes larger than the neutrino mass, 
such that the neutrino decay is no longer possible and is replaced by the decay of Higgs bosons 
into neutrinos and leptons. We discuss in detail the conceptual and numerical differences of the full 
two-mode approach to the one-mode approach and the vacuum result. The CP-asymmetry for the 
different approaches is the main topic of chapter 4. The CP-asymmetry in the two-mode approach 
consists of four different contributions due to the two possibilities for the leptons in the loops. 
We present some useful rules for performing calculations with the fermionic modes and compare 
the analytical expressions for the CP-asymmetries in different cases. We restrict ourselves to the 
hierarchical limit where the mass of Ni is much smaller than the mass of A''2 , that is M2 ^ Mi . The 
temperature dependence of the CP-asymmetry is discussed in detail for the one-mode approach, 
the two-mode approach and the vacuum case. The differences between the asymmetries and the 
physical interpretation of certain features of the asymmetries are explained in detail. Chapter 5 
deals with the evaluation of the Boltzmann equations. We derive the equations and explicitly 
perform the subtraction of on-shell propagators for our cases in appendix E. We compare our four 
thermal scenarios, wich are the decoupled and strongly coupled two-mode approach and the one- 
mode approach with thermal and asymptotic mass, to the vacuum case. We show the evolution 
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of the abundances for three different initial conditions for the neutrinos, that is zero, thermal and 
dominant abundance. We explain the dynamics of the different cases in detail and find considerable 
differences both of the thermal approaches to the vacuum case and of the two-mode cases to the 
one-mode cases. We summarise the main insights of this work in the Conclusions and give an 
outlook on future work and prospects. 

In appendix A, we present Green's functions at zero temperature, while in appendix B, we derive 
the analytical solution for the lepton dispersion relations. Some quantities relating to leptogenesis 
in the vacuum case are derived in appendix C. In appendix D, we present analytical expressions 
for the CP-asymmetry contributions of the two cuts through {N',i'} and {N',(p'}^, which we did 
not consider in chapter 4, since we are working in the hierarchical limit. Appendix E is devoted to 
the detailed description of a correct subtraction of on-shell propagators in our scenario. 



We shamelessly stole our notation for the cuts in the vertex contribution from reference [22]. 
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Leptogenesis 



1.1 The Matter- Antimatter Asymmetry 

The matter-antimatter asymmetry of the universe is usuaUy expressed as 

UB -riQ 



r] = 



n^ 



il.l) 



where ns, n^, and n-y are the number densities of baryons, antibaryons, and photons, respectively, 
and the subscript imphes that the value is measured at present cosmic time. There might be an 
excess of leptons over antileptons as well, but its contribution to the energy density of the universe 
is small compared to the contribution of the baryons. 

The photon density is proportional to T^ and the temperature of the universe is inferred via 
the cosmic microwave background radiation (CMB), which shows an almost perfect blackbody 
spectrum. The baryon density can be inferred in two ways: First, from the abundances of the 
light elements D, ^He, ^He and ^Li, which are a direct probe of the primordial abundances formed 
during the big bang nucleosynthesis (BBN) phase at redshifts z '^ lO"*^". Of these, the deuterium 
abundance depends most sensitively on the baryon-to-photon ratio rj, while the other elements 
exhibit a weaker dependence. A measurement of the abundances gives [24] 

?7 = (5.8 ± 0.7) X 10"^° (1.2) 

at 95% confidence level, see figure 1.1. Second, there exist even more stringent constraints on ij, 
which originate in the observation of the CMB anisotropies. The anisotropics reflect the acoustic 
oscillations of the photon-baryon fluid at the time of decoupling, which in turn depend on the 
baryon-to-photon ratio at this time, at a redshift z ~ 1000. A high baryon density enhances the 
odd peaks in the angular power spectrum relative to the even peaks, as can be seen in figure 1.2. 
The measurement of the CMB anisotropies from the 7- year Wilkinson Microwave Anisotropy Probe 
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Figure 1.1: The observed abundances of light elements compared to the standard BBN predic- 
tions [24]. The smaller boxes indicate 2cr statistical errors only, the larger ones the combined 2cr 
statistical and systematic errors. 
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Figure 1.2: The dependence of CMB temperature anisotropies on the baryon abundance fib = 
Pb/ Pcrit compared with data [26], where pciit is the critical energy density of the universe. The 
relation to 77 is given by 77 = 2.74 x 10"^ ^b h'^, where h = iJo/(100 km s'^ Mpc^^) = 0.72 ± 0.08 is 
the present Hubble parameter. 



(WMAP) data gives [25] 



rj= (6.15 ±0.25) x lO" 



-10 



(1.3) 



The agreement of these two indicators at extremely different redshifts is an important success 
of standard big-bang cosmology. One might be tempted to assume this asymmetry as an initial 
condition of the universe, but as we hinted at in the introduction, there are at least two grave 
arguments against such an assumption. The first argument is that such an initial condition would 
be a highly fine-tuned one, since it would imply that for every 100 million quark-antiquark pairs, 
there would have been one additional quark. The second argument comes from the broad agreement 
that the universe has undergone an inflationary expansion in its very early phase that solves some 
otherwise unresolvable problems from cosmological observations. The inflationary phase would 
wash out any asymmetry that might have existed after the big bang. Thus, we need a theory that 
creates the baryon asymmetry dynamically after inflation. Such theories are commonly referred to 
as baryogenesis theories. 



1.2 Sakharov's Legacy 



Sakharov has formulated three necessary conditions [27] that baryogenesis theories have to fulfil : 



1. 



^There are exotic models that do not fulfill all of these conditions, but still produce a baryon asymmetry, such as 
Dirac leptogenesis for example, which does not require a violation of lepton number [28] . 
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• Baryon nuniher B non-conservation: 

It is immediately clear that processes which create a baryon asymmetry do not conserve 
baryon number B. 

• C and CP symmetry violation: 

Processes that are invariant under the discrete transformations of charge conjugation C or the 
product of charge conjugation and parity reversal, CP, will not create a baryon asymmetry. 
This is because in these cases, the B violating processes that create baryons proceed at 
the same rate as their C- and CP-conjugated processes which create an equal amount of 
antibaryons. Thus we need C- and CP-symmetry violation. 

• Deviation from thermal equilibrium: 

In chemical equilibrium, the entropy is maximal if the chemical potentials associated with 
non-conserved quantum numbers, such as B in our case, vanish. Since also the masses of 
quarks and antiquarks are equal by CPT-invariance, the phase space densities 

1 + eVP +™?/^ 
are equal in thermal equilibrium and thus also the number densities. 

1.3 Why Does the Standard Model Fail? 

In principle, all three Sakharov conditions are met in the standard model of particle physics (SM). 

1.3.1 B non-conservation 

Baryon and lepton number, B and L, are accidental symmetries in the SM and are conserved at tree 
level. However, at the one-loop level, B and L are not conserved in anomalous diagrams where the 
B and L currents couple to two electroweak gauge bosons through a fermion triangle [29-32]. In the 
vacuum structure of the electroweak theory, there are non-perturbative transitions from one vacuum 
to a different vacuum with differing B and L numbers [33]. These instanton transitions change 
B and L by three each and conserve B — L. At zero temperature, the instantons are suppressed 
by the instanton action, i.e. exp(— Svr^/^i^), where g is the SU(2) coupling constant, so that their 
rate is negligible. However, there are static, but unstable field configurations which correspond to 
saddle points of the field energy of the gauge-Higgs system, so-called sphaleron configurations [12]. 
These saddle points can be reached thermally and mediate a thermal transition from one vacuum to 
another vacuum. At temperatures above the electroweak scale, T > Tew ~ 100 GeV, the sphaleron 
processes are fast and in equilibrium [13]. Thus, they will also wash out any previously existing B 
or L asymmetry that is not linked to a B — L asymmetry. 



1.3 Why Does the Standard Model Fail? 



1.3.2 C and CP violation 

The weak interactions violate C niaxinially since only left-handed particles and right-handed an- 
tiparticles couple to the gauge bosons, whereas their charge conjugated states, left-handed an- 
tiparticles and right-handed particles do not couple via SU{2)l gauge interactions. Only a charge 
conjugation together with a parity transformation {CP) converts left-handed particles into their 
right-handed antiparticles. However, even CP-symmetry is broken in the SM in B- and i^-meson 
systems through a CP-violating phase in the quark mixing matrix [34]. If appropriately nor- 
malised [35], this CP violation is still several orders of magnitude too small as to produce an 
asymmetry of the order i] ~ 10^^". This is one reason why baryogenesis does not work in the 
SM without any additional assumptions. Therefore, any reasonable baryogenesis theory needs 
additional sources for CP violation. 



1.3.3 Deviation from thermal equilibrium 

In the early universe, there may be a departure from thermal equilibrium during the electroweak 
phase transition [36,37], which is in principle suitable to create a baryon asymmetry. The deviation 
from equilibrium occurs in particle interactions through the bubble walls between the broken and 
the unbroken phase. In order to obtain an irreversible asymmetry, the potential barrier between 
the two phases has to be large enough, that is the phase transition has to be strongly first order. 
This imposes constraints on the Higgs potential, which in turn relate to an upper bound on the 
Higgs mass. However, the experimental lower bound on the Higgs mass is too high as to allow for 
this kind of phase transition and we arrive at the second reason why baryogenesis does not work in 
the SM. The theory requires an additional mechanism to obtain a sufficient departure from thermal 
equilibrium. 

1.3.4 Ways out: baryogenesis theories 

We see that successful baryogenesis needs two new ingredients: First, new sources of CP violation 
and second, a different mechanism for a departure from thermal equilibrium^. Moreover, any 
mechanism that creates a B or L asymmetry at higher temperatures also has to violate B — L, 
otherwise this asymmetry will be washed out by the sphaleron interactions. Several possibilities 
to meet these requirements have been investigated, such as leptogenesis [7], grand unified theory 
(GUT) baryogenesis [40-47], electroweak baryogenesis [38,39], the Affieck-Dine mechanism [48,49] 
and other, more exotic variants (see, e.g. reference [50]). Leptogenesis is the model that we focus 
on in this work. 



^Or additional degrees of freedom that allow for a first order phase transition, like in the minimal supersymmetric 
standard model (MSSM) with a light stop particle [38,39]. 
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1.4 Leptogenesis 

1.4.1 The unbearable lightness of neutrino masses and the seesaw as a way out 

The attractiveness of leptogenesis [7] arises from the feature that in addition to creating a baryon 
asymmetry, it simultaneously solves a seemingly unrelated puzzle, which is the smallness of neutrino 
masses. If one turns the argument around, by employing a seesaw mechanism [8-11] in order 
to explain why neutrinos have a non-zero, but extremely small mass, we naturally arrive at a 
mechanism for the generation of the baryon asymmetry without additional effort. 

There is experimental evidence that at least two neutrinos have a non-zero mass which is several 
orders of magnitude smaller than the masses of the charged fermions. Oscillation experiments have 
established two mass-squared differences between the neutrino mass eigenstates and global fits 
give [24] 

Amii =m\-m\ = (7.59 ± 0.20) x 10"^ eV^ = m^„i, 

(1.5) 
jAmlal = \r(i\-m\\ = (2.43 ± 0.13) x 10"^ eV^ = m^^^ . 

Neutrinos are the only SM fermions for which a Majorana mass term is in principle allowed since 
they do not carry a U{1) charge, 

^rn, = -j^cArrAapVp + h-c , (1.6) 

where Va are the neutrino fields and the superscript c denotes charge conjugation. The subscripts 
a and (3 denote the neutrino flavour and [m]^/? is the Majorana mass mixing matrix. This operator, 
however, is not invariant under the SU{2)l gauge group, so the simplest operator which respects the 
symmetry of this gauge group and reduces to equation (1.6) upon spontaneous symmetry breaking 
is the dimension five operator {i^ci4>){^i3'P)i where 

is the lepton doublet with fiavour a and 

frh+\ 

(1.8) 




the Higgs doublet. This operator in turn is not renormalisable, so it must be the effective remnant 
of new physics that is realised at a higher energy scale, in order not to spoil the renormalisability of 
the theory. Thus, if there is new physics above the electroweak scale, it will induce this dimension 
five operator at lower energies unless some symmetry prevents it. This observation is a strong 
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argument in favour of Majorana mass terms. Moreover, Dirac mass terms of the form mVnUL 
require the addition of right-handed neutrinos at low energy that are singlets under the SM gauge 
group and whose existence could only be inferred via the exclusion of a low energy "Majorana" 
type mass. If the dimension five operator is induced via tree-level interactions with heavy particles 
at a mass scale M, which can be much higher than the electroweak breaking scale MeW) this will 
automatically lead to a light neutrino mass scale of M^-^/M for Yukawa couplings of order one, 
which explains the smallness of neutrino masses. Since these heavy particles can be viewed as 
suppressing the mass of the neutrinos, the mechanism is called seesaw mechanism. 

Since the heavy particles have to couple to a Higgs doublet and a lepton doublet, there are 
three prominent possibilities, which are called seesaw type I-III: 

• Type I: The heavy particles are S'C/(2)2,-singlet fermions [8-11]; 

• Type II: The heavy particles are S't/(2)-triplet scalars [11,51-54]; 

• Type III: The heavy particles are S'?7(2)-triplet fermions [55-57]. 

The type I seesaw is the simplest and the framework in which leptogenesis is usually implemented, 
so we concentrate on this type. 

We add two or three singlet fermions iVj to the SM, sometimes referred to as "right-handed 
neutrinos", which are assumed to have rather large Majorana masses Mj, close to the scale of some 
possibly underlying grand unified theory (GUT), -Bgut ~ 10^^, 10^^ GeV. The additional terms of 
the Lagrangian can be written in the mass basis of the charged leptons and of the singlet fermions 
as 

c = i m^-i^'N, - x^amreabt) -^Yl Mi^Ni + h.c. , (1.9) 

i 

where the Higgs doublet (j) is normalised such that its vacuum expetation value (vev) in 

{<p) = r] (1.10) 

is f ~ 174 GeV and Xia is the Yukawa coupling connecting the Higgs doublet, lepton doublet and 
heavy neutrino singlet. The indices a and b denote doublet indices and €ab is the two-dimensional 
total antisymmetric tensor which ensures antisymmetric S'C/(2)-contraction. 

The effective mass matrix [m]^/? of the light neutrinos as defined in equation (1.6) can be written 
as 

[m\a(s = Yl KkM^^XfSkV^ . (1.11) 

k 

It can be diagonalised as 

Dm = U^[m]U, (1.12) 
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where Dm = d\a,g{mi^m2-,'m^) and U is the leptonic mixing matrix, also cahed Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix. It is a 3 x 3 unitary matrix and therefore depends, in general, 
on six phases and three mixing angles. Three of the phases can be removed by redefining the phases 
of the charged lepton doublet fields. Doing this, the matrix U can be conveniently parameterised 
as 

C/ = f/-diag(l,e'",e^'^), (1.13) 

where a and /3 are called Majorana phases. If the neutrinos had Dirac mass terms, the Majorana 
phases could be removed by redefining the phases of the neutrino fields. The matrix U can be 
parameterised in a way similar to the Cabbibo-Kobayashi-Maskawa (CKM) matrix, 

(C13C12 C13S12 siae""^ \ 

-C23S12 - S23Sl3Cl2e"^ C23C12 - S23Sl3.Sl2e"^ ■S23C13 , (1-14) 

S23C12 - C23Sl3Cl2e"^ -■S23C12 - C23Sl3Sl2e"^ C23C13 / 

where Qj = cos Oij and Sij = sin Oij and 6ij are the angles of rotations in flavour space, which 
connect the flavour basis with the mass basis. 

Thus there are twelve physical parameters at low energies in the leptonic sector of the SM if we 
add the mass matrix of equation (1.11): the three charged lepton masses rrie, m^, m^, the three 
neutrino masses mi, m2, fn^,^ and the three angles and three phases of the PMNS matrix U . Of 
these parameters, seven have been measured, me-, m^, rrir, Am^i, lAmggI, O12 and ^23- There exists 
an upper bound on the angle ^13, however, the mass of the lightest neutrino and the three phases 
of U are experimentally not accessible at the moment. 

In the case of three heavy neutrinos, there are nine additional parameters in the high-energy 
theory, amounting to 21 parameters in total. These high-energy parameters cannot be measured 
at experimentally accessible scales, but are nevertheless important for leptogenesis. Moreover, if 
one assumes only two right-handed neutrinos [58], there are fourteen parameters in total and one 
of the light neutrinos is massless such that its corresponding phase also vanishes. This amounts to 
ten physical parameters at low energy and four additional parameters at high energy. Such "two- 
right-handed-neutrino" (2RHN) models are attractive since they have strong predictive power. We 
will concentrate on the case of three heavy neutrinos in this work, since it allows for the possibility 
that all light neutrinos have mass and the number of generations is the same as at low energy. 

1.4.2 Sakharov and leptogenesis 

We examine whether and how the Sakharov conditions are fulfilled in leptogenesis. The heavy 
neutrinos only couple to the lepton and Higgs doublets, thus the processes Nj — t- li(f) and Nj — )■ ii(j) 
are the only possible tree-level decays. Since the Lagrangian in equation (1.9) violates L, there is 
more than one possibility to assign a lepton number to the heavy neutrinos. One usually assigns to 
them a lepton number of zero, then the decays violate lepton number and also, very importantly. 
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B — L hy one unit. A lepton and B — L asymmetry can thus be generated at high temperatures, 
which win be converted to a baryon asymmetry by the sphaleron processes at the electroweak scale. 

Charge conjugation C is violated maximally in the SM and the Yukawa couplings X^k can have 
CP-violating phases. However, CP violation can only arise via an interference of the tree-level 
decay and its higher order corrections, most importantly the one-loop contribution. 

Concerning the third Sakharov condition, the heavy neutrinos will decouple from thermal equi- 
librium when the expansion rate of the universe is faster than the interaction rate, i.e. if the decay 
rate F of the neutrinos is smaller than the Hubble rate H. They will in this case nevertheless decay 
into lepton and Higgs doublet, but the decay will be out of equilibrium. 



1.4.3 A simple model 

We see that the Sakharov conditions can in principle be fulfilled in leptogenesis, it remains to be 
determined, which lepton and baryon asymmetry will be generated for which parameters of the 
high energy theory. 



Creating a lepton asymmetry 

In order to introduce the calculations one has to perform in leptogenesis, we will make some 
simplifying assumptions and present a model which serves as an example for determining the 
leptogenesis dynamics. Firstly, we assume that the masses of the heavy neutrinos are strongly 
hierarchical, i.e. Mi <^ M2 <^ M3. This corresponds to the hierarchical masses of the standard 
model fermions and simplifies the calculation. If the reheating temperature is larger than the masses 
of the two heavier neutrinos, one cannot neglect the effect of producing the A'^2 and N^. Thus, in 
order to simplify calculations, we assume a reheating temperature which is lower than these masses 
but larger than Mi in order to allow for a thermal production of the leptogenesis protagonists. 
Furthermore, we will, again for the sake of simplicity, also neglect flavour effects, even though they 
play an important role. We assume that the lepton states ^at^ , into which the Ni decay, will keep 
coherence of their flavour mixing until there are no A^i processes. This is only true if leptogenesis 
happens at a temperature larger than about 10^^ GeV, below which processes mediated by the 
r- Yukawa coupling become fast [59,60]. However, flavour effects depend on the Yukawa couplings 
and do not necessarily have to be important below these temperatures. Moreover, the thermal 
effects which are examined in this thesis can be extended to include flavour effects if appropriately 
modified. This work aims at giving an insight into the effect of quasiparticles, where adding flavour 
might be an important second step in the future. We refer the reader who wishes to learn about 
flavour effects and the possible influence of A'^2 and N^ to the review by Davidson et al. [14]. He or 
she will find an extensive list of references therein. 
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Figure 1.3: The CP-asymmetry A^q-^i- The graph in the middle is the self-energy contribution, 
the graph on the right the vertex contribution. In the self-energy diagram, the Higgs boson and 
lepton can be particles or antiparticles. 



The CP-asymmetry 



The CP-asymmetry in a lepton flavour i that arises in the decay of a heavy neutrino of generation 
j is defined as 

The tree-level amplitude cannot be different for the CP-conjugated process, but if one considers 
the one-loop level, the interference terms in the squared sum of the matrix elements for tree level, 
A^O; and one-loop level, A^i, can be CP-asymmetric, i.e. M.o M.^ / M.o M.I, where Aii denotes 
the CP-conjugated matrix elements. 

There are two contributions to the one- loop amplitude Mi as displayed in figure 1.3. One is the 
self-energy contribution where a virtual N^ is exchanged in the s-channel. The second one is the 
vertex-correction graph. There is only a CP-asymmetry if the virtual N^ is a different generation 
than Nj, i.e. A'^2 or N^ for the decay of iVi. The CP-asymmetry is proportional to the imaginary 
part of the one-loop diagram which in turn corresponds to assuming an on-shell condition for the 
loop propagators. In the self-energy diagram, one can only put the Higgs boson and lepton on the 
mass shell and A^^ is necessarily off-shell since M^ ^ Mj. In the vertex correction, it is kinematically 
not possible to put the neutrino in the loop on its mass shell since it is heavier than the Higgs 
boson and the lepton. Thus, again, the Higgs boson and the lepton are on-shell. We will see in 
chapter 4 that at finite temperature, also the neutrino can be on its mass shell. 

The calculation of the CP-asymmetry is worked out in Appendix C.2. Including all contribu- 
tions, it is written as 



1 



Svr (AtA)jj 



Im 



Kji^''^)jk>^ik gi^kj) + 



Svr (AtA)jj 



Im 



X*j{VX)kjKk 



1 



•^kj 



where 



xkj = mI/m] 



(1.16) 



(1.17) 
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and 

1 -. ., ^, fl+x 



9{x) 



+ 1 - (l + x)ln 



^^^ --^ (1.18) 



— X \ X 

and we have summed over the lepton generations in the loop. The second term in equation (1.16) 
corresponds to the self-energy diagram with particles in the loop as opposed to antiparticles. It is 
lepton flavour changing but lepton number conserving. Moreover, it does not exhibit an imaginary 
part if we sum over the final state lepton flavors i, so that 

e = X]^^ = ^(AtAV"^™ 1 ''^^^''^'' fdi^kj)- (1.19) 

The dominant contribution to the CP-asymmetry in A^i decays, which we consider, comes from A'^2 
in the loop since M3 ^ M2. In our further analysis, we will neglect the effect of M3 and set j = 1 
and k = 2. In the limit of hierarchical right-handed neutrinos, it is possible to derive an upper 
bound on the CP-asymmetry in A'^i decays, the so called Davidson-Ibarra-bound [61], 



3 Ml, , 3 Ml 



e^^^("^3-"ii)«— ^^Am2,^, (1.20) 



since in this limit, the light neutrino spectrum is expected to be hierarchical as well. 



Anticipating the discussion of the next two sections, it is possible to parameterise the resulting 
baryon asymmetry by three numbers as 

r] = deK, (1.21) 

where ci is a factor which accounts for the dilution of the asymmetry through the expansion of the 
universe and the distribution into different particle species. It is of the order ~ 0.01. The efficiency 
factor K parameterises the dynamics of leptogenesis and how much of the asymmetry is washed out 
by other processes and is for a vanishing initial lepton asymmetry given by 

n = ^^^^. (1.22) 

It is calculated by solving the Boltzmann equations that govern the evolution of the particle species 
and usually of the order ~ 0.1. Since we have to explain an asymmetry of r/ ~ 5 x 10~^*^, we can, 
for typical efficiency factors, require a CP-asymmetry of e '^ 10~^. Together with the Davidson- 
Ibarra-bound, this translates to a lower bound on the mass of the lightest heavy neutrino of about 

Ml > ^ J!!^ 6^109GeV. (1.23) 
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The gravitino problem 

If one assumes that the heavy neutrinos are produced thermally, a reheating temperature of 
^RH > Ml is needed in order to produce a sufficient amount of neutrinos and this can lead to the 
so-called gravitino problem, which we discuss by closely following reference [14]. In certain models 
of supersymmetry (SUSY), a high reheating temperature leads to an overproduction of gravitinos. 
The gravitinos are long-lived and will decay into lighter particles if they are not the lightest su- 
persymmetric partner (LSP). If too many gravitinos decay during or after BBN, the decays will 
destroy the agreement between predicted and observed light element abundances [62-65]. There 
are several possibilities to avoid such a scenario, both on the SUSY side and on the leptogenesis 
side. The possibilities on the SUSY side include the following: 

1. The gravitinos decay before BBN. Heavy enough gravitinos arise in anomaly-mediated sce- 
narios [66], but not in standard gravity mediated scenarios. 

2. Late-time entropy production can dilute the gravitino abundance, but it also dilutes the B — L 
asymmetry [67] 

3. The gravitino is the LSP (and the dark matter particle), in which case the bound on Trh is 
less restrictive [68,69]. The gravitino is indeed the LSP in gauge-mediated scenarios. 

We assume that either SUSY is not realised in nature, which does not pose a problem to our 
leptogenesis model, or that one of the above scenarios solves the gravitino problem. This way, it is 
possible to consider temperature regimes T > Mi. 



Converting a lepton to a baryon asymmetry 

The L asymmetry, which is generated at high temperature, is partially converted into a B asym- 
metry by the electroweak sphaleron processes, which are in equilibrium above ~ Tew- Taking into 
account the chemical potentials of all particles, it can be worked out (e.g. reference [14]) that 

Yb = -^Yl = cYb-l , (1.24) 



where 



Yx = ^ (1.25) 

s 



is the number density nx of species X over the entropy density s and c = Yb/Yb-l accounts for 
the sphaleron conversion oi a, B — L asymmetry into a B asymmetry. In the SM, c = 12/37 for a 
smooth electroweak phase transition. 
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Solving Boltzmann equations 

The intuitive requirement for the neutrino decays to be out of equihbrium is that they decay slower 
than the universe expands, that is the decay rate 

Svr Sttv^ ^ ' 

is smaller than the Hubble rate H, where 

is the conveniently defined effective neutrino mass, which is of the order of the light neutrino mass 
scale. In the language of this mass, the out-of-equilibrium condition reads 



v" 



Ml 



= m* ~ 1.1 X lO^^eV, (1.28) 

T=Mi 



where m* is called equilibrium neutrino mass. The parameter region where this condition is satisfied 
is called strong washout regime. However, leptogenesis is also possible when fhi > m* . This can 
be seen if one solves the set of Boltzmann equations [70, 71] 

^ = -iD + S)iNr,,-N^^\), (1.29) 

^^^ = -eD{Nr,,-N';,\)-WNB-L, (1-30) 

where Nx is the number density of the species X per comoving volume which contained one photon 
when T S> Mi. The temperature is contained in z = Mi/T and {D,S,W) = {T£,,Ts,Tw)/{Hz), 
where Td is the decay rate, Ts the sum of all A'^i scattering rates, and Tw the sum of all processes 
that wash out the generated L asymmetry, such as inverse decays for example. We derive a simple 
version of equations (1.29) and (1.30) in appendix C.3. 

Solving these equations typically gives efficiency factors of k ~ 0.1. The baryon asymmetry can 
be approximated as 

where the first factor is the equilibrium Ni number density divided by the entropy density at 
T ^ Ml. The number of relativistic degrees of freedom </* is given by 

i=bosons i=fcrmions 

where i denotes species with mass rrii <^ T and the factor 7/8 arises from the difference in Fermi 
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and Bose statistics [72]. At the temperature of leptogenesis, all SM particles have negligible masses, 
so 5f* = 106.75. The expressions r] and Yb are related via photon and entropy density today as 

■So 7.1 



CHAPTER 2 



Thermal Field Theory 



The topic of this work is the role of the equihbrium quantum effects that are imphed by the 
presence of a hot, dense medium for leptogenesis models. The influence of the medium is studied 
by means of an effective, statistical quantum field theory, which takes into account the temperature 
of the surrounding medium, hence called thermal field theory (TFT). In this chapter, we give an 
introduction into TFT and explain the methods and formalisms we use later on. The reader who 
wishes to learn more about TFT is referred to the books by Kapusta [73], LeBellac [74] and Das [75]. 
Our presentation follows the more intuitive approach of Thoma's lecture notes [76]. 



2.1 Green's Functions at Finite Temperature 

Thermal field theory comprises all three basic branches of modern physics, namely quantum me- 
chanics, relativity and statistical physics. We want to derive Feynman rules and diagrams at finite 
temperature. To this end, we consider the two-point Green's function or propagator. We would 
like to find an expression for this quantity at finite temperature. For simplicity, we consider the 
case of a scalar field (p. The propagator at T = is worked out in Appendix A and is the vacuum 
expectation value of the time-ordered product of two fields at spacetime points x and y, 

iAF = {0\T^.(l)y\0). (2.1) 

At finite temperature, the vacuum expectation value of an operator A has to be replaced by the 
quantum statistic expectation value of the corresponding statistical ensemble, 

(0|^|0)^(^)p^tr(p^), (2.2) 
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where p is the density operator of the statistical ensemble. For the canonical ensemble, which we 
will consider, it is given by 

P = ^e-^^ , (2.3) 

where 

and H is the Hamiltonian operator of the system with the discrete eigenvalues and eigenstates 

H\n) = En\n). (2.5) 

The partition function is 



Z = tr(e~'^^), (2.6) 



so we can write 



{A), = |tr(^e-^^) = I j;(n|^|n)e-^^" , (2.7) 

n 

where the sum is over all thermally excited states \n) of the system, which are eigenstates of 
the Hamiltonian H. Since the heat bath is a distinguished reference frame for our situation, our 
calculations are not Lorentz-invariant. We will always work in the rest frame of the heat bath, 
which is the preferred rest frame at finite temperature. 

Calculating the statistical expectation value for the scalar propagator, we have 

\^l>\x -y) = \ j;(n|r{</>(x)</>(y)}|n)e-^^" . (2.8) 

n 

Using the Fourier representation of ^ in equation (A. 3) in the propagator for the case xq > yo gives 
1 r d^fc d^/fc' 1 1 



iAl>'{x-y) = ^l 



Zj (2vr)3/2 (2vr)3/2 (2^,)i/2 (2c.^)i/2 ^^_^^ 



Y^ e-^^- {n\ [a(k)e"^ ^'^ + a) (k)e ' ^'^l [a(k')e-' ^''^ + a^ (k')e ^ ^'-^l | 



n) 



The states 

|n) = |ni(ki),n2(k2),...,ni(ki)) (2.10) 

are orthonormalised states with nj bosons of momentum kj. Acting with the creation and destruc- 
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tion operators on \n) according to equation (A. 4) results in 

d^fc Sk' 1 1 



'^?"<^-^'4/(l 



2 



)3/2 (2^)3/2 (2^,)l/2 (2^/ )i/: 

X Y, e-'^-^-ii^lk) + l]<5^(k - k')e"'-^-^+'-^'-^ + n(k)(53(k - \^)e'^-^-''^'-y] 

n 

\ ( A^k 1 



(27r)3 2wfe 



^e-'^-^"{[n(k) + l]e-'^-(^-^) + n(k)e'^-(^-^))} 



(2.11) 



We use 



^ E-(k)e-^^" = ^^^ ^ Mu.,) , (2.12) 

n 

where £'„ = X^fcWfcn(k) and fsi^k) is the Bose-Einstein distribution. Using this relation, we find 

(^^{[1 + /B(^fe)]e-^(^-^) + /B(u;fc)e'^(--^)} . (2.13) 

We see that for T = 0, which implies that fs = 0, we arrive at the vacuum result of equation (A. 7). 
We can interpret this result as follows: The zero-temperature part describes the usual propagation 
of a particle from y to x. However, on top of spontaneous creation at y, there is also induced 
creation at x (proportional to ub) and absorption at x (oc ub) due to the presence of the thermal 
particles in the bath. The propagator in equation (2.13) is the propagator of a free field without 
interactions, and forms the starting point for perturbation theory when we add interactions. 



2.2 Imaginary Time Formalism 

The propagator (2.13) is a useful quantity, but not very helpful in constructing Feynman rules. We 
need a representation which consists of a 4-dimensional i^-integration, so that we can formulate 
Feynman rules in momentum space. To achieve this, we continue the propagator analytically to 
imaginary times t with < t = it < f3 and sum over the discrete energies 

ko = 2TTmT, (2.14) 

the so-called Matsubara frequencies instead of integrating, so that 



/ 



i-'^E- p-i'^) 
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The propagator can be written as 

.A?>«W=irj:/^-^e---. (2.16) 



(27r)3 K2 _ ^2 

One can also motivate the introduction of imaginary time in the following way: for r = it — t- /?, 
the Boltzmann factor exp(— /?//) has the form of the time evolution operator exp{—iHt). As a 
consequence, thermal propagators become periodic in /3, 

A|;>0(x -y) = A|;>0(x, y, r, 0) = A^^x, y, r, /3) , (2.17) 

where we have taken Tx = t and Ty = 0. In general 

A|;>0(r) = A^>0(T + n/3) (2.18) 

holds for any integer n. 

The two consequences of this relation are: 

1. The time r is restricted to the interval [0,/?], the so-called Kubo-Martin-Schwinger- or KMS- 
condition. 

2. The Fourier integral over k^ in vacuum quantum field theory (QFT) becomes a Fourier series 
over the Matsubara frequencies ko = 2iTmT. (For fermions we have ko = (2n + l)i7rT, as 
described in section 2.4. 

In real life, it is usually necessary to integrate over more than one propagator. In these cases, 
it is difficult to perform the summation over the zeroth component of the loop momentum ko. 
A convenient way out of this problem is to use the so-called Saclay representation, which is a 
mixed representation, performing the Fourier transformation in time only. In the following, we use 
A = A^^^ and write the propagator in the Saclay representation, leads to 

A(T,a;,) = -TY^^-^'^MK), (2.19) 

fco 

where the Fourier coefficients are given by 

rl3 
A{K) = - dre'=«^A(T,a;fc). (2.20) 

We can perform the sum (2.19), 

AiT,0Jk) = -^ {[1 + fBiu;k)]e-'''^ + /ij(^fc)e"'="} , (2.21) 
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which agrees with equation (2.13). 

It is often convenient to write the propagator as the sum 

AiT,ojk) = 5] ^[1 + fB{su;k)]e~'^''^ = J] A,(T,a;fc), (2.22) 

s=±l ^^ s=±l 

where we allow for negative energies sw^ in /^(sa;^). In frequency space, the two parts also 
decompose, 

A(i^) = y A,(K) = y — - — i , (2.23) 

where the relations 

A,(T,a;fc) = -T5^e-*^°^A,(K), 

As{K) = - dre'=«^A,(r,Wfc) (2.24) 

Jo 

hold for the propagator parts as well. 



2.3 The Scalar Field 

We write down the Feynman rules for a simple interaction theory, the neutral scalar field with a 
0^ interaction, given by the Lagrangian 



^{d,<P){d,<P)-l, 



C = - {d^cp) {d^cp) --rr?^^-g^^\ (2.25) 



where is a real field and ra the mass of the corresponding particles. We can derive the Feynman 
rules for the interaction of these fields. One can follow the operator formalism with which we 
started in the derivation of the finite temperature Green's functions in section 2.1 and which relies 
on the interaction picture and a thermal equivalent of Wick's theorem at zero temperature. As 
usual in QFT, there is a second approach, which is to derive the Feynman rules from the path 
integral formalism. Both derivations can be found in the literature [73-75]. 
The Feynman rules for this theory read as follows: 



1. The propagator is given by 

A^ — m^ 
with /cg = 27ri nT. 



--^V'm = ^P^-3. (2-26) 
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Figure 2.1: The (/)'*-tadpole. 



2. In loop integrals we make the replacement 






3. The vertex reads as in vacumn — i4!(7^. 

4. Symmetry factors, e.g. 1/2 for the tadpole, are the same as in vacuum. 

As a simple example of a loop diagram, we compute the tadpole of the (/)^-theory shown in 
Fig. 2.1. According to the above Feynman rules, we get 

n = il(-i4!,^)ir^/^.A(ir). (2.28) 



The Saclay representation (2.20) proves useful in this calculation. 



n = Vlg" 1-0^1^ dr A(t, a;,)r ^ e^^o^ , (2.29) 



since the Matsubara frequency fco, over which we have to sum, appears in the propagator only in 
the exponent, which makes the summation simple at the expense of introducing another integral 
over T. The sum reduces to a J-function, 



oo 



T Y, e^o(^-^') = ^ Yl e2^^"^(^-^') = T6{T{t - r')) = S{t - r') . (2.30) 

n=—co n=—oo 

Thus we get 

n = 1252|^A(0,u;fc) (2.31) 

= 65^ //^-[1 + 2/^(0;,)]. 
The first term in the sum [1 + 2fB{oJk)] leads to an ultraviolet (UV) divergence, but is identical 
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to the corresponding zero-temperature term. The second term results in a finite integral since /b 
falls off exponentially fast for large momenta. We see that the tadpole can be decomposed into a 
vacuum and a finite temperature part, 

n = n^=o + n^>o , (2.32) 

and it is sufficient to renormalise the divergence of the vacuum term. The finite temperature part 
can be integrated analytically if m = and in this case yields the simple result 

n^>o = g^T^ . (2.33) 

In fact, one can show that generally, renormalisation at zero temperature is sufficient to remove 
all UV divergences of the theory at finite temperature. We do not prove this statement here, but 
the interested reader can find a proof in chapter 3.5 of Bellac's book for example [74]. However, 
the property can be understood in a more intuitive way: temperature does not modify the theory 
at distances much smaller than 1/T and thus, the ultraviolet divergences are the same as at zero 
temperature. Infrared (IR) divergences, on the other hand, are a different story and we will deal 
with them in chapter 2.5. 

In more complex calculations, we often have to perform the Matsubara sum over two or more 
propagators, so it is useful to write down these frequency sums for further reference: One such sum 
is given by 

r5^A(K)A(P-i^)=T5^5^A,,(K)A,,(P-K), (2.34) 

feo fco *1'«2 

where we can in principle allow the two bosons corresponding to the two propagators A(P — K) 
to have different masses mi and m2- We can perform this sum by using the Saclay representa- 
tion (2.20), transforming the Matsubara sum into a 5- function and executing the imaginary time 
integrals. We arrive at^ 

rV A.,(i^)A,,(P -K) = - '''' 1 + /M^ic^i) + /b(^2^2) (2.35) 



fco 



AUkl^q PO - SlUl - S2U}2 



where ooi = yk^ + mf and UJ2 = \/{p — k)^ + m^ are the energies of the fields. 



^The frequency-space propagators in Bellac's book [74] are defined with a minus sign relative to our convention 
^thiswork|-^-j _ —/xBciiaci-j^-j jj^ defiues the Fourier transformation for the mixed representation with a minus sign, 
so the mixed-representation propagators agree with this work. Moverover, the frequency sums involve exactly two 
propagators, so they agree with this work as well. 
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Another frequency sum, 

TY,koAsAK)^sAP-K), (2.36) 

fco 

can be evaluated by integrating the Fourier representation of A(r) by parts [74]. We arrive at 

T J: fcoA., WA,,(P - K) = _I^ ^ + fBis.-i) + fBis..,) ^2_3^^ 

^ 4:Ll!2 Po - SlLOl - S2UJ2 

fco 
so we see that the calculation amounts to replacing /cq by the respective propagator pole sw^. 



2.4 The Dirac Field 

We are considering spin 1/2 particles in a spinor representation iIj{x) with the free Lagrangian 
density 

Coire^c = i^{il''d^ - m)iP , (2.38) 

where we work in the chiral or Weyl representation of gamma matrices. 
We define a fermion Matsubara propagator through 

Sjsix -y) = {T{il,^{x)^s{y)})p , (2.39) 

where 7 and 6 denote spinor indices and p indicates that we are taking the statistical expectation 
value for a density operator p, for which we take the canonical ensemble (2.3). 

Similar to the scalar field, the fermion propagator obeys a KMS condition when going to imag- 
inary time t with r = it, however, the fermion propagator is antiperiodic in /3, 

S{x -y) = 5(x, y, T, 0) = -5(x, y, r, /3) , (2.40) 

where Tx = t and Ty = 0. The minus sign arises when changing the time order of ip and ^p because 
the spinors anticommute. In general 

5(r) = (-l)"5(r + n/3) (2.41) 

holds for any integer n. 

Because of the antiperiodicity, the Matsubara frequencies are given by 

fco = (2n + l)i7rr, (2.42) 
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where n is an integer. Analogous to equation (2.16), one can derive the Fourier representation of 
the Matsubara propagator, 

ko 



If we write the propagator as^ 



^W = J-^ = (^ - m) • A{K) , (2.44) 

then the mixed representation is given by 

A(r, Uk) = -TY^ e-'^'-AiK) , (2.45) 

fco 

where 

K{K) = - dre'^«^A(T,a;fc). (2.46) 

Jo 

As in the scalar case, one can perforin the sum in equation (2.45), 

A(r,a;fc) = t^ {[1 - friu^k)] e~"'=" - /F(^fc)e"'=0 (2-47) 



s=±l 



^— [1-/Hsa;fc)]e— '=- 



^ As(r,tjfc) 



where 

/H-.) = ^^ (2.48) 

is the Fermi-Dirac distribution and we allow for negative energies sujf^. As in the scalar case, the 
two parts decompose in frequency space as well, 

K=Y, Ks{K) = Y, -^-^ , (2.49) 



^Note that A(_ft') = 1/{K'^ —m?) is not the same as the scalar propagator A(A') because the Matsubara frequencies 
fco = (2n + l)i7rT are different for fermions. 
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where again the relations 



and 



A,(T,a;fc) = -T^e-^'"^A3(K) (2.50) 



Ks{K) = - dre'=«-A,(r,Wfc) (2.51) 

JO 



hold. 



It is straightforward to calculate the following four basic frequency sums as in Eqs. (2.35) and 
(2.37): 

1. Fermion-boson case: 

TY ^,,{K)K,,{P - K) = - '''' I + fBJsiu^) - MS2U.,) ^ 

^ 4uJkUJq po - SlWl - S2UJ2 

^ 4aj2 Po - siuji - S2UJ2 

feo 



2. Fermion-antifermion case: 



ry A,AK)A,AP - K) = - '''' ^ - fFis.u^i) - fFis2U.2) ^ 

^ AUkUJq Po - SiWi - S2UJ2 

TY^koA,AK)KiP - K) = _I^ ^-Ms,u.,)-Ms2U2) 

^ AUJ2 Po - SiUJi - S2UJ2 



We can obtain Eqs. (2.52) - (2.53) by replacing fB{soj) for a bosonic line by —fp{suj) for a fermionic 
line. The substitution /b(/co) ~^ ~fF{ko) is in fact a systematic rule for substituing a fermionic line 
for a bosonic line in an arbitrary Feynman graph. We will not prove this rule, but the interested 
reader is referred to the treatment by, e.g. Bellac [74]. 

For gauge fields, it can be shown [74] that the propagator reads 

DAK) = j^i^2 (2-54) 

in Feynman gauge. We will not explicitly employ the gauge field propagator in the future discus- 
sion, even though it is used in deriving the thermal masses of the Higgs bosons and leptons in 
references [77-80], to which we refer in section 3.4. We quote the Feynman gauge, other gauges 
can be found in reference [74]. 
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2.5 Hard Thermal Loop Resummation 

One might expect that with the formahsm and techniques of the preceding sections, it is possible 
to calculate all diagrams in all finite temperature field theories. However, using the perturbation 
theory as described, one encounters the following serious problems: 

1. IR divergences: A famous example is the energy loss of a heavy fermion in a plasma [81,82]. 

2. Gauge dependent results: An example of this is the gluon damping rate in a hot QCD plasma, 
which turns out to be different in different gauges [83,84] 

3. Power counting: It turns out that the resummation of infinitely many higher order diagrams 
can contribute at a lower order in perturbation theory than expected. 

In order to cure or at least alleviate these problems, the so called hard thermal loop (HTL) 
resummation technique has been invented in the late 80s and in the beginning of the 90s by Braaten 
and Pisarski [85,86]. Instead of using bare propagators and vertices, they suggested to use effective 
vertices and effective propagators, constructed by resumming certain diagrams, the so-called HTL 
self energies. In this way an improved perturbation theory has been established, which we briefiy 
present. 

2.5.1 HTL self energies 

Our first step is to isolate the diagrams that should be resummed into effective propagators. To 
this end, we distinguish between two scales. In a hot plasma, there are two momentum scales, 
the hard scale T and the soft scale gT, where we assume for now that the coupling constant g is 
much smaller than one. One could add more scales, such as g^T for example, but two scales will 
be sufficient for our discussion. The diagrams to be resummed are the HTL self energies, one-loop 
diagrams in which the external momenta are soft (oc gT) and all internal momenta hard (oc T). 

For the scalar self energy with the (j)'^ interaction in Fig. 2.1, the result in the HTL limit is the 
same as the full result of equation (2.33). However, in the case of a fermion, the bare self energy 
is gauge dependent [87] like the gluon self energy, whereas in the HTL limit we obtain a different, 
gauge independent result, which we present in the following. 

In all practical calculations of this thesis, the bare mass of fermions will be negligible compared 
to the temperature, so we study only the case of massless fermions. The plasma introduces the 
rest frame of the heat bath as a special Lorentz frame. In a general frame, the heat bath has 
four-velocity n" with u^Uq = 1. In the rest frame of the plasma, we can write n" = (1,0, 0, 0) and 
|i = 7o. The general expression for the self-energy in the rest frame of the thermal bath is given 
by [77] 

S(P) = -a(P)i^-6(P)|i, (2.55) 
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0-T 




Figure 2.2: The HTL fcrmion self-energy 



where the factors a and h are given by 

a{P) = -^ [tr (J^S) - poti (70S)] , (2.56) 

^(^) = 4^ [P'tr (70$]) -potr (/>$])] , 
where the traces are evaluated in the HTL approximation [74], and one finds 

Ti =tr(i^S) =4ml 

^2 = tr (7oZj) = 2mp-ln (2.57) 

P po-p + ie 

with the effective fermion mass 

, e^rVs for QED 
^F={ ^ ' _ (2.58) 

g T /16 for a Yukawa interaction £y = —gil^il^cj) . 

We will not need a resummed gauge boson propagator, therefore we will not quote the gauge 
boson self energy. Again, the interested reader will find information in the literature [74]. 



2.5.2 Effective propagators and dispersion relations 

We can construct effective propagators which lead to an improved perturbation theory by re- 
summing the HTL self energy diagrams. For a scalar field with a self energy H, the resummed 
propagator is given by the Dyson- Schwinger equation in Fig. 2.3. This diagrammatic equation 
reads 

iA* =iA-hiA(-in)iA-h... 
i i 

" A-i-n " K2 - m2 - n ' ^^-^^^ 
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Figure 2.3: The rcsummcd scalar propagator. 



where K is the momentum of the propagator and m the zero temperature mass. The dispersion 
relation for a particle is given by the pole of its propagator, and for the bare propagator we get 

K"^ -m"^ = k'^-uj'^{k)=0, (2.60) 

that is 



ko = u;{k) = VhF+^ . (2.61) 

However, at finite temperature we get a different effective propagator for a collective mode with 
an effective mass nies = V^ + H, where it is often possible to neglect the zero temperature mass 
m with respect to the self energy such that meg = vIT. The dispersion relation of the collective 
scalar particle is then given by 



ojcs = y k^ + mis . (2.62) 

Effective masses and dispersion relations as above, generated by the interaction with a medium, 
have been introduced in various areas in physics, such as the effective mass of an electron in a 
crystal or the reduced velocity of a photon in a medium. 

Considering fermions, we restrict ourselves to the case where the bare mass is negligible and, 
in a similar way as in equation (2.59), we get for the effective propagator 

S*(K) = — i , (2.63) 

where Shtl is given by Eqs. (2.55)-(2.58). It is very convenient to rewrite this propagator in the 
helicity-eigenstate representation [88,89], 

S*{K) = ^A+(i^)(7o - k • 7) + ^A_(i^)(7o + k • 7), (2-64) 
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Figure 2.4: The two dispersion laws for fermionic excitations compared to the standard dispersion 
relation o;^ = fc^ + m|n. 



where k = k/fc, 



A±iK) 



ko±k + 



■nil 



±1 



±kn 



2k 



— ln^° 



kn — k 



(2.65) 



and mp is the effective fermion mass defined in equation (2.58). 

This propagator has two poles, the zeros of the two denominators /S.±. The poles can be seen as 
the dispersion relations of collective excitations of the fermions that interact with the hot plasma, 



ko = uj±{k). 



(2.66) 



We have found an analytical expression for the two dispersion relations making use of the Lambert 
W function [2], which is calculated in appendix B. The dispersion relations are shown in Fig. 2.4. 

Note that even though the dispersion relations resemble the behaviour of massive particles and 
uj = mp for zero momentum fc, the propagator S*{K) (2.64) does not break chiral invariance. Both 
the self energy S(i^) (2.55) and the propagator S*{K) anticommute with 75. The Dirac spinors 
that are associated with the pole at /cq = oj+ are eigenstates of the operator (70 — k • 7) and they 
have a positive ratio of helicity over chirality, x = +1- The spinors associated with k^ = a;_, on the 
other hand, are eigenstates of (70 + k • 7) and have a negative helicity-over-chirality ratio, x = ~1- 
At zero temperature, fermions have x = +1- The introduction of a thermal bath gives rise to 
collective fermionic modes which have x = — 1- These modes have been called plasminos since they 
are new fermionic excitations of the plasma and have first been noted in references [77,78]. 

We can introduce a spectral representation for the two parts of the fermion propagator (2.65) 
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[90], 

A^iK)=rduJ^±^, (2.67) 

where the spectral density p±{Lj,k) [88,91] has two contributions, one from the poles, 

p^^%io,k) = Z±{uj,k)d{uj - uj±{k)) + Z^{uj,k)6{uj + uj^{k)) , (2.68) 

and one discontinuous part, 

pf'iu, k) = \ml{k^u:)^ _ ^ ^^^2 _ ^2) (2.69) 

{k{uj^k) -m|,[Qo(a;)=FQi(a;)]} + [\ tt m?p{l ^ x)] 

where x = uj/k, 9{x) is the heaviside function and Qq and Qi are Legendre functions of the second 
kind, 

Qo(x) = ^ In ^^ , Qi (x) = xQo(x) - 1 . (2.70) 

2 X — 1 

The residues of the quasi-particle poles are given by 

Z±{uj,k) = ^^^ '~ — , where Z+ + Z_ = 1 . (2.71) 

One can describe the non-standard dispersion relations uj± by momentum-dependent effective 
masses m±{k) which are given by 



m±{h) = yj oj'iik) - k"^ = ^j2Z{Lo,k)mF ■ (2.72) 

These masses are shown in Fig. 2.5. 

2.5.3 HTL resummation technique 

We have collected the necessary ingredients to build an improved perturbation theory at finite 
temperature. We noted that naive perturbation theory suffers from the problems of IR divergences 
and gauge dependent results. The reason for this is that the naive perturbative expansion is 
incomplete at T > 0. Infinitely many higher order diagrams can contribute to lower order in the 
coupling constant. These diagrams can be taken into account by resummation. 

We will not discuss the HTL resummation technique for gauge theories in detail (see, e.g. [74]), 
but we will present rules for perturbative calculations. One has to consider the self energies in the 
HTL approximation, such as H, S* and the gauge boson polarisation tensor H* ^. Due to Ward 



As mention above, we will not need the gauge boson propagator in this work, but we mention it in order to sketch 
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Figure 2.5: The momentum-dependent effective masses m±. 




Figure 2.6: The HTL vertex correction 



identities, fermion self energies are related to vertices, e.g. 



ie [S(Pi) - S(P2)] = (Pi + P2)/.r^(Pi,P2) . 



(2.73) 



Therefore one also has to consider the HTL correction to vertex contributions'^, as shown in fig- 
ure 2.6, where all internal lines are hard ~ T. The HTL propagators are constructed by resumming 
the self energies via the Dyson-Schwinger equation as explained above. The HTL vertices are given 
by adding the HTL correction to the bare vertex. Examples are shown in figure 2.7. When calculat- 
ing diagrams, we have to use effective propagators and vertices if all external legs are soft; otherwies 
bare propagators and vertices are sufficient. In this way, contributions of the same order in g are 
included, gauge independent results are obtained and the IR behaviour of the theory is improved. 



the HTL resummation technique. 

^Even though we wiU not encounter effective vertices in the course of this work, we add them in order to present 
a fuU picture of the HTL theory. 
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Figure 2.7: Examples of an effective HTL propagator and an effective vertex. 



After all, the HTL improved perturbation theory has been successfully applied to thermal QCD 
for the description of the quark-gluon plasma (see e.g. [92]) 

Having outlined the construction of the HTL resummation technique as a perturbation theory 
that cures many serious shortcomings of naive perturbation theory at finite temperature, some 
remarks concerning our use of HTL propagators seem necessary. First, even though it is sufficient 
to use bare propagators and vertices if one external momentum is hard, it is always possible to 
resum self-energies and thus capture effects which arise from higher-order loop diagrams and thus 
take into account the appearance of thermal masses and modified dispersion relations in a medium. 
In fact, since the effective masses we will encounter do typically not satisfy the condition nics ^ T 
but are rather in the range m^^/T ~ 0.1 - 1, the effect of resummed propagators is noticeable even 
when some or all external momenta are hard ~ T. In summary, we always resum the propagators 
of particles that are in equilibrium with the thermal bath, that is in our case the Higgs bosons 
and the leptons, in order to capture the effects of thermal masses, modified dispersion relation and 
modified helicity structures. This approach is justified a posteriori by the sizeable corrections it 
reveals, similar to the treatment of meson correlation fuctions in reference [93]. 



CHAPTER 3 



Decays and Inverse Decays 



3.1 The Quest of This Thesis 

In order to make any statement about the amount of baryon asymmetry that is produced by a 
baryogenesis or, in our case, leptogenesis model, one has to adopt a quantitative description of 
the dynamics that take place at this phase and result in generating the asymmetry. There are 
two ways to do this: One can either adopt the consistent quantum mechanical view of the system 
and calculate how quantum systems evolve with time when they are not in equilibrium. Or one 
can view the particle distributions as classical distributions and adopt the equations which govern 
them, that is, a set of Boltzmann equations. The set of equations that govern the dynamics 
of the first non-equilibrium approach are the so called Kadanoff-Baym equations [94] and this 
approach has received some attention in recent publications [17,19-21,23,95-100]. The easier and 
more traditional Boltzmann approach is possible when interactions are not fast and particles are 
sufficiently close to their equilibrium distribution, such that before and after each interaction, they 
can be seen as classical particles. The interaction itself is calculated by means of quantum field 
theory and appears in the collision term of the equation. 

Whichever viewpoint we adopt, we have to calculate the quantum mechanical amplitudes that 
enter the equations. Traditionally, amplitudes are calculated in vacuum and inserted into Boltz- 
mann equations. However, as densities and temperatures are high, it is important to calculate 
amplitudes at finite temperature and compare the result to the zero-temperature case. In this 
work, we adopt the Boltzmann view of particle evolution and calculate the corrections due to finite 
densities and temperatures. However, the amplitudes are related to the two-point functions one 
has to use when adopting the Kadanoff-Baym approach. 

The Boltzmann equations, which we stated for zero temperature in equations (1.29) and (1.30), 
contain two main quantities, which are calculated from amplitudes: First, the rates with which 
interactions take place, that is, the decay and inverse decay rates To and F/^), the scattering rates 
Ts and the washout rates T]^. Second, the difference between the L violating rates and their 
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CP-conjugates, the CP-asymmetry, which is defined in equation (1.15) for zero temperature. Of 
these, the leading order amphtudes in the interaction rates are tree level amplitudes. However, 
as explained in section 1.4.3, the CP-asymmetry arises as an interference between tree level and 
one- loop amplitudes and is proportional to the imaginary part of one loop diagrams. Thus, for the 
CP-asymmetry, the leading order is the one-loop level. 

The combination of SM couplings that we write as g in the HTL-corrections gT is typically of 
the order g ~ 0(10^^). Since the creation of the lepton asymmetry takes places at temperatures 
T ~ Ml, corrections of the order of gT are important to consider if the couplings are as large as 
they are in our case. We present a consistent way of calculating these HTL-corrections of order 
gT and analyse the effects they imply. Thus, we restrict ourselves to a very basic leptogenesis toy 
model which is self-consistent and contains all important features of leptogenesis even if it does 
not describe the full scenario in a quantitatively accurate way. More specifically, we look at the 
leading order interactions between the protagonists of leptogenesis, the heavy neutrinos A'^i, as 
well as the lepton and Higgs doublets i and (j). These interactions are decays Ni — t- LH, inverse 
decays LH — )■ A^i and the AL = 2 scatterings LH — )■ LH, where by L and H, we denote both the 
doublets (j),i and their charge conjugated states (p,i. Naturally, we have to include a calculation of 
the CP-asymmetry, which is a dominant quantity in the sense that there would be no asymmetry 
production without it. We include HTL-corrections for the Higgs bosons and the leptons, but not 
for the neutrinos, since the Yukawa couplings are much smaller than the SM couplings that give 
rise to the HTL effects for Higgs bosons and leptons. 

As temperatures are high, particles in the bath acquire different dispersion relations as explained 
in section 2.5.2, which can be translated into effective thermal masses. Due to these masses, it may 
happen that some interactions are kinematically forbidden, while other processes that would not 
be possible in vacuum are allowed. In leptogenesis, there is a temperature where the processes 
A^i -H- LH are forbidden, but the processes H -H- LNi are possible. The leptons always have a 
lower thermal mass than the Higgs bosons, so L -H- NiH is never allowed. These Higgs decays and 
inverse decays are new processes that govern the neutrino and lepton evolution at high temperature, 
so they need to be calculated and likewise the CP-asymmetry in these decays since it leads to the 
production of a lepton asymmetry. 

Thermal corrections to leptogenesis of various kind have been studied before, both in the 
Boltzmann-picture [15,16,18,22,101-105], as well as in the Kadanoff-Baym-picture [17,19-21, 
23,95-100]. Our self-consistent calculation is a novel approach and we discuss the differences to 
the aforementioned works, with a special focus on the extensive and important work by Giudice 
et al. [16]. In this chapter, we calculate the leading order HTL corrections to decays and inverse 
decays, in chapter 4, we do the same for the CP-asymmetries and in chapter 5, we put our results 
into the appropriate Boltzmann equations and evaluate them. 
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Figure 3.1: N decay via the optical theorem with dressed propagators denoted by a blob 

3.2 Discontinuity of the Fermion Self-Energy in Yukawa Theory 
at Finite Temperature 

We consider a leptogenesis-inspired model with a massive Majorana fermion N coupling to a mass- 
less Dirac fermion i and a massless scalar (j). The interaction and mass part of the Lagrangian 
reads 

Ant.mass = gN^Pi - ^MNN'' + h.C. , (3.1) 

The HTL resummation technique has been considered in [106] for the case of a Dirac fermion with 
Yukawa coupling, from which the HTL resummed propagators for the Lagrangian in equation (3.1) 
follow directly. We calculate the interaction rate F of A^ -f-T- i(p. 

We cut the A^ self-energy and use the HTL resummation for the fermion and scalar propagators 
(figure 3.1). According to finite-temperature cutting rules [107,108], the interaction rate reads 

r(P) = --^ tr[(f + M) Im S(po + ie, p)]- (3.2) 

2po 



At finite temperature, the self-energy reads 

d^A; 

ko=i {2n+l)nT 



nP) = -g^T Y. J J^PLS*{K)PnD*{Q), (3.3) 



where Pl and Pr are the projection operators on left- and right-handed states and Q = P — K. 
The HTL-resummed scalar propagator is 



1 
Q2 _ m'^ 



D*iQ)= ^, „2 ^ (3-4) 



where m^ = g^T"^ /12 is the thermal mass of the scalar, created by the interaction with fermions, 
and can be calculated analogously to the (/>^ self-energy in section 2.3. Due to the reduced Majorana 
degrees of freedom, m^ differs from the Dirac-Dirac case by a factor 1/2 [106]. 

The effective fermion propagator in the helicity-eigenstate representation is given by equa- 
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tion (2.64) [88,89], 








S*{K) = - 


--A+(K)(7o - k • 7) - -A_(K)(7o + k • 7 




A^{K) = 


" / _L ; , "^? /^_Li ±^0-^1 ko + k\ 
[ ^°^^+ k (^^ 2k ^"fco-fcJJ 


and 




mj = y'T\ 



(3.5) 



(3.6) 

This again differs from the Dirac case in equation (2.58) by a factor 1/2 [106]. 
The trace can be evaluated as 

tr[(i^ + M)PlS*{K)Pr] = -A+{po - pv) - A_(po +pr?), (3.7) 

where r] = cos 9 is the angle between p and k. We evaluate the sum over Matsubara frequencies 
by using the Saclay method [109]. For the scalar propagator, the Saclay representation from 
equation (2.20) reads 

D*{Q) = - f dre*'--^{[l + fniu^^Jje-^"^ + /B(^,)e-'^}, (3.8) 

where Ug = q^ + m'^. For the fermion propagator, it is convenient to use the spectral representation 
as explained in section 2.5.2 [90], 

A±{K) = - / dr'e'^"^' / dco p±(^, k)[l - Moj)]e"^^' , (3.9) 

where p± is the spectral density [88]. 

Since the quasi-particles are our final states, we will set K such that l/A-|-(i^) = 0. Thus, we 
are only interested in the pole contribution 

pl^'^iiv, k) = -'^^^—^{6{u; - a;±) + 6{cv + lo^)), (3.10) 

2m| 

where uj± are the dispersion relations for the two quasiparticles, i.e. the solutions for ko such that 
1/A lb (w-t, k) = 0, shown in figure 2.4. The analytic solutions for oj± are explained in appendix B. 
One assigns a momentum-dependent thermal mass m±{k)'^ = io±{k)'^ — k"^ to the two modes as 
explained in section 2.5.2 and shown in figure 2.5 and for large momenta the heavy mode m+ 
approaches \^mi, while the light mode becomes massless. 

In order to execute the sum over Matsubara frequencies, we write /cq = ^^n with aj„ = (2n+l)7rT 
and remember that, when evaluating frequency sums, also po = '^^m = i (2m + 1)ttT can be written 
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as a Matsubara frequency and later on be continued analytically to real values of pQ [74, 110, 111]. 
In particular e^^^ = e''^™^ = — 1. We write 



^E^'""'= E ^ir-n'f3) 



T Y^ e^P" 



-ko)T kor' _ POT r/ / 



gPor j(^/ _ ^) 



(3.11) 



since —(3<t' — t</3. After evaluating the sum over /cq and carrying out the integrations over r 
and t' , we get 



/oo 1 



PO-U-UJg 



PQ-UJ + UJq 



Integrating uj over the pole part of p± in equation (3.10), we get 



"^ 2a;„ 1 2rnj 



1 {ujl-k^ \l + fB-fF , fB + fp 



+ 



+ - 



2mj 



PO-U}±-UJq Po-U}±+UJg 
fB + fp ^l + fB-fp 



_P0 + CJzp - W<j po + W=p + Uq 

where fs = /^(wg) and fp = fp{uj±) or fp{uj^), respectively. 



(3.12) 



(3.13) 



The four terms in equation (3.13) correspond to the processes with the energy relations indicated 
in the denominator, i.e. the decay N — )■ (f)i, the production N(f) — )■ i, the production Ni — t- (p and 
the production of Nicf) from the vacuum, as well as the four inverse reactions [107]. We are only 
interested in the process N -H- (pi, where the decay and inverse decay are illustrated by the statistical 
factors 



l + fB-fp = il + fB)il - fp) + fBfp . 



(3.14) 



The decay is weighted by the factor (1 + fB){^ — fp) for induced emission of a Higgs boson and a 
lepton, while the inverse decay is weighted by the factor fsfp for absorption of a Higgs boson and 
a lepton from the thermal bath. Our term reads 



r^Z)*A± 



4 



k^ l + fB-fF 



N^<pe 



1 OJ^ 

2ijjq 2mi po — UJ± — OJq 



(3.15) 
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For carrying out the integration over the angle r], we use 

Im -— = -7rS{po -u}±- ujg) = -TT--^5{r] - r/±), (3.16) 

Po-uj±-ujq + ie ^ kp 

where 

''^ " 2r [^^°'^^ -M^- (4 - k^) + ml] (3.17) 

denotes the angle for which the energy conservation po = oj + cjq holds. The integration over r] then 
yields 

/ dT? Im(T j;Z)*A±) = --^ ^^ [1 + /b(^,±) - /F(a;±)], (3.18) 

j_i ^^^ ZKp zm^ 

where ujq± = po — uj± . It follows that 

r(P) = - -^ tT[i^ + M)Im S(P)] 
1 I ^-^ f H^A- 

= - ^^ Im i T^ / d/cd?? A;22)*[A+(po -p??) + A_(po +P^)] 



(3.19) 



fell 

,2 /■ ,.,2 _ 7,2 






327rpoP — j_i<^±<i '^nq 

X [2po(A: T oj±) ± M^ ± (4 - fc^) ^ m|], 

where we only integrate over regions with — 1 < r/ < 1. 



Using finite temperature cutting rules, one can also write the interaction rates for the two modes 
in a way that resembles the zero-temperature case [107] 

r±(P) = r^ / dfcdg i27T)HHP -K-Q) \M±{P, K)\'' 

2po J (3.20) 



where 
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and dq analogously and the matrix elements are 

2 1,2 

\M±{P,K)\^ = g^ '^^^~, u;±{poTPV±)- (3-22) 



Now that we have arrived at an expression for the full HTL decay rate of a Yukawa fermion, we 
would like to compare it to the conventional approximation adopted for example by reference [16], 
which we refer to as one-mode approximation. To this end, we do the same calculation for an 
approximated fermion propagator 

'-'approxl-'^) = w- _ ' (3.23) 

This yields the following interaction rate: 

rapprox(P) =^£^ // ^^ ^[^ + ^^^'^^^ " M^^^^M^ + "^' - "^'] 

= J_ f dk dq {2^f5\P -K- Q)|7W|2 (3-24) 

2po J 

x[1 + /bK)-/fH], 
where w^ = /c^ + m?, ojq = pQ — uj and the integration boundaries 



ki. 



1 



2 



2M2 ^0 V ^^' + ""' " ""^^^ " (2^"^^)' T P(^ + "^1 - "^0) (3-25) 

ensure — 1 < ?? < 1, where 

^ = ^ [2pow - M2 - m^2 + ^2] _ (326) 

We see that the matrix element is 

l-MP = y(M2 + m2-m2). (3.27) 

In addition to the dispersion relations and the phase space boundaries for fc, there are two major 
differences to the two-mode matrix element in equation (3.22). In the two- mode approach, we 
integrate over the residue Z±{k) = (w^ — k'^)/{2m'j). For the plus-mode, this residue is mostly 
close to unity, but can be as low as 1/2 for low momenta. For the negative mode, the residue is 
close to zero for most momenta and only up to 1/2 for low momenta. The rate for the plus-mode 
is slightly suppressed compared to the one-mode approach, while the rate for the minus-mode is 
considerably suppressed. 
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Another difference is the momentum product 

P-K = pouj-pkr] = - {M^ + mj-ml) = -. (3.28) 

For the two-mode approach, we can introduce a chirahy invariant four-momentum for the lepton 

Kl^ = ojh{l,hk), (3.29) 

where /i = ±1 denotes the hehcity-over-chirahty ratio. Then 

y 

P ■ Kh= pqu - hpujTjh = h— - po{hu} - k) , (3.30) 

where 

T.h = M^ + ml{k)-ml. (3.31) 

We see that for the plus- mode, the term po('^ ~ k) has to be subtracted from the right-hand side 
of equation (3.28), which also suppresses the rate compared to the one-mode calculation. For the 
minus-mode, the momentum product is also different from the one-mode approach. This difference 
in the momentum products is closely linked to the fact that the modified dispersion relation of the 
two-mode calculation leaves the chiral symmetry unbroken. 

Continuing our discussion of the one-mode rate, we note that it resembles the zero temperature 
result 

rT=o(P) = ,^ r dfe - [m2 + m} - ml] (3.32) 

with zero temperature masses m^, m^. The missing factor 

l + fB-fF = {l + fB){l-fF)+fBfF (3.33) 

accounts for the statistical distribution of the initial or final particles. As pointed out in more detail 
in reference [1], we see that the approach to treat thermal masses like zero temperature masses 
in the final state [16] is justified for the decay rates, since it equals the HTL treatment with an 
approximate fermion propagator. However this approach does not equal the full HTL result. 

Concluding this calculation, a caveat has to be added: In this general calculation, the external 
Majorana fermion will also acquire a thermal mass of order gT. Thus, if its zero temperature 
mass is smaller than its thermal mass, the external fermion also needs to be described by leptonic 
quasiparticles to be consistent. However, in the leptogenesis study, the Yukawa coupling giving 
rise to the Majorana neutrino decay is much smaller than the couplings giving rise to the thermal 
masses of the Higgs boson (scalar) and the lepton (Dirac fermion) and thus the thermal mass of 



3.3 Decays at High Temperature 45 

the heavy neutrino can be neglected, as pointed out in the previous section. 

We have calculated the decay rate assuming a Majorana particle, but the result can be very 
easily generalized to the case of two Dirac fermions by inserting the appropriate factors of two in 
the decay rate and the thermal masses. 

3.3 Decays at High Temperature 

When the temperature is so high that m(f,{T) > M, the scalar can decay into the Majorana fermion 
and the Dirac fermion^. The calculation can be done in the same way as for the Majorana fermion 
decay. The only difference is that in equation (3.13), we take the imaginary part of the factor 
l/(po + '^=F ~ '^9)' which corresponds to the scalar decay. The frequency sum corresponding to 
equation (3.15) reads 



fco 



2ujq 2mj pQ + uj-h - Wg ■ 



where /i = ±1 denotes the helicity-over-chirality ratio. In this case, the angle t] is given by^ 

^° = ^ [-2po^h + S^] , (3.35) 

where 

Y.^ = m\-M^ - {ojI - k^) . (3.36) 

In order to clarify the momentum relations, we revert the direction of the three-momenta q and p 
so that they correspond to the physical momenta of the incoming scalar and outgoing Majorana 
fermion. The matrix element can be derived as 

\Mh{P,K)\'' = g^ZhU^hiPo -hp-k), (3.37) 

where the momentum flip p — > — p compensates the helicity flip h ^ —h and 

Zh = \^ (3.38) 

is the residue of the modes and the angle for the reverted physical momenta reads 

Vh = ^ i^Pou^h - S^] . (3.39) 



^Note that in our model calculation, m^{T) > mi{T) 

^Note that in this notation, — q and — p are the three-momenta of the initial-state scalar and the final-state 
Majorana fermion, since their roles have been inverted. 
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3.4 Application to Leptogenesis 

When turning to leptogenesis with 

bC = iN.d^-i^'Ni - Xu,ia^i(pUa - ^MiN.N^ + h.c, (3.40) 

we sum over the two components of the doublets, particles and antiparticles and the three lepton 
flavors. Thus we need to replace g^ by 4(AlAjy)ii. Integrating over all neutrino momenta, the decay 
density in equilibrium is 

^^' = / (lf3/^"(^)rD = ^,J^dEEpf^-rn, (3.41) 

where E = pQ, f^{E) = [exp(£'/3) — 1]^"*^ is the equilibrium distribution of the neutrinos and 
Td = [1 — /"^(E)] r. We drop the subscript 1 for the neutrino in this section since it is the only 
occuring neutrino. The decay density is also present in the Boltzmann equations, 

7(7V ^HL) = J dpNdpHdpL{27T)^6\pN - PL - Ph) \Mh\^ /^^(l - 1^(1 + f^) • (3.42) 

For the Higgs boson decay, the density reads 

7(F ^ NL) = j dpMdpHdpL{27:)H\pN+PL-PH) \Mh\^ (1 - /^^)(1 - r^)Q, (3.43) 

where the matrix elements for both decays are related and given by equations (3.22) and (3.37). 
The thermal masses are given by [77-80] 

mJ(T)=(i..| + i..?. + i.? + iA)T^ 

™?(T) = (4.1 + ^4) T-'. (344) 

The couplings denote the SU(2) coupling 52, the U(l) coupling gy, the top Yukawa coupling yt and 
the Higgs self-coupling A, where we assume a Higgs mass of about 115 GeV. The other Yukawa 
couplings can be neglected since they are much smaller than unity and the remaining couplings are 
renormalised at the first Matsubara mode, 27rT, as explained in reference [16] and, in more detail, 
in reference [112]. 

In figure 3.2, we compare our consistent HTL calculation to the one-mode approximation 
adopted by reference [16], while we add quantum-statistical distribution functions to their cal- 
culation, which equals the approach of using an approximated lepton propagator l/(^ — m^) as in 
equation (3.23) [1]. In addition, we show the one-mode approach for the asymptotic mass \/2m£. 
We evaluate the decay rates for the Mi = 10"*^*^ GeV and normalise the rates by the effective neu- 
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Figure 3.2: The decay densities for the neutrino and the Higgs boson decay. We show the one-mode 
approach with the thermal mass as 7™^ and with the asymptotic mass as 7^„ ; Also the T = 
rate 70 and our two modes "f±. The temperature thresholds are explained in the text. 



trino mass rrii, which is often taken as mi = 0.06 eV, inspired by the mass scale of the atmospheric 
mass splitting. 

In the one-mode approach, the decay is forbidden when the thermal masses of Higgs boson and 
lepton become larger than the neutrino mass, Mi < m^ + m^ or Mi < \plmi + m^. Considering 
two modes, the kinematics exhibit a more interesting behavior. For the plus-mode, the phase space 
is reduced due to the larger quasi-mass, and at Mi = m_(_(oo) -|- m^, the decay is only possible 
into leptons with small momenta, thus the rate drops dramatically. The decay into the negative, 
quasi-massless mode is suppressed since its residue is much smaller than the one of the plus-mode. 
However, the decay is possible up to Mi = m^. Due to the various effects, the two-mode rate differs 
from the one-mode approach by more than one order of magnitude in the interesting temperature 
regime of z = T/Mi > 1. The v2 ?7i^-calculation is a better approximation to the plus-mode, but 
still overestimates the rate, which is mainly due to the difference between the momentum products 
in equations (3.28) and (3.30), that is the helicity structure of the quasiparticles. The residue also 
reduces the plus-rate, but the effect is smaller since Z+ is usually close to one. 

At higher temperatures, when m^ > Mi + m±{k), the Higgs can decay into neutrino and lepton 
modes and this process acts as a production mechanism for neutrinos [16]. The decay (p — t- Ni^ 
is possible when m^ > Mi, while the decay into ^-|_ is possible when m^ > Mi + m£. As for 
low temperature, the rate 7+ is unsuppressed only when ?ti^ > Mi + \'zmi. Our decay density 
approaches the decay density of reference [16] at high temperatures, but is about a factor two below. 
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This can be explained by the fact that the phase space is smaher due to the larger mass of the lepton, 
rrii < m(:{k) < \jlmi. Again, the asymptotic mass calculation is a better approximation but still 
gives a larger rate due to the momentum product and, to less extent, the residue. We see that the 
decay rate rises as ~ T^, instead of ~ T^ as for the vacuum rate 70. In the vacuum calculation, 
the squared matrix element is proportional to M^. In the finite temperature calculation, it is 
proportional to S^ = m\ — rnj^k) — M^, so the dominant contribution is proportional to m| ~ T^ 
and the rate rises by a factor T^ faster than the vacuum rate 70. 

Summarising, we can distinguish five different thresholds for the thermal decay rates we dis- 
cussed. Going from low temperature to high temperature, these are given by the following condi- 
tions: 



rpN 
rpN 



Tc 



rn<P 



Ml = v2 m£ + rriff, , 

Ml = mi + m^ , 

Mi=m^, 

m^ = mi + Ml, 

m^ = V2me + Mi. (3.45) 



We will refer to these thresholds in the following chapters. 

As discussed in detail in [1], we confirm by employing HTL resummation and finite temperature 
cutting rules that treating thermal masses as kinematic masses as in [16] is a reasonable approxima- 
tion. However, quantum statistical functions need to be included as they always appear in thermal 
field theory. Moreover, the non-trivial two-mode behaviour of the full HTL lepton propagator is 
not accounted for by the conventional one-mode approach. We have calculated the effect of the two 
modes in a general way, which is applicable to any decay and inverse decay rates involving fermions 
at high temperature. Thus, this calculation is a valuable tool for other particle interaction rates in 
the early universe, as other leptogenesis processes, the thermal production of gravitinos or the like. 

The behaviour of the decay density of the two lepton modes can be explained by considering the 
dispersion relations uj± of the modes and assigning momentum-dependent quasi-masses to them. 
The thresholds for neutrino decay reported in reference [16] are shifted and the decay density shows 
deviations of more than one order of magnitude in the interesting temperature regime T/Mi ^^ 1, 
which has implications for the dynamics of leptogenesis as we will see in chapter 5. 



CHAPTER 4 



CP- Asymmetries 



While the decay rates calculated in the last chapter are the first important ingredient of the Boltz- 
mann equations, the CP- asymmetry is the second quantity that enters the Boltzmann equations. 
This chapter is devoted to the calculation of the CP-asymmetry at finite temperature for our dif- 
ferent scenarios. At low temperature, we calculate the CP-asymmetry in neutrino decays, at high 
temperature the CP-asymmetry in Higgs boson decays 



4.1 Preliminaries 

We are calculating the CP-asymmetry in A^i decays. We denote the decaying Ni by N and the N2 
in the loop by N' . At T = 0, the CP-asymmetry is defined in equation (1.15), 

_ r(JV ^ </,£) - T{N ^ ^i) 

where F are the decay rates of the heavy A^s into Higgs boson and lepton doublet and their CP- 
conjugated processes. As we will see in chapter 5, we have to calculate the CP-asymmetry via the 
integrated decay rates at finite temperature, 

'^ ' 7^>^{N ^ 04) + 7^>°(Ar ^ 04) ' ^ ' ' 

where we define the CP-asymmetry for each lepton mode, denoted by h, as introduced in equa- 
tion (3.29). We have 



7^>°= i^fNiPN)r^>'iPj:,), (4.3) 
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where /at is the distribution function of the neutrinos and Pj^ the neutrino momentum. At T = 0, 



we write 



Ml 

r{pn = — rrf, 

Po 



(4.4) 



where Mi and pQ are the mass and the energy of the neutrino and Frf is the decay rate in the rest 
frame of the neutrino. The integration over the momentum cancels out and the CP asymmetry 
via 7 is the same as via T. At finite temperature, however, the thermal bath breaks Lorentz 
invariance and the preferred frame of reference for calculations is the rest frame of the thermal 
bath. The momentum dependence of the decay rate cannot be formulated as in equation (4.4) and 
the CP-asymmetry as defined in equation (4.1) is momentum dependent. Since we want to decouple 
the CP-asymmetry from the momentum integration in the Boltzmann-equations, the definition in 
equation (4.2) is the appropriate one. 

The CP asymmetry in equilibrium can be written as 



eq 



(T) 






(4.5) 



where Fd = T{N ^ tcfy and Td = T{N 
rate. 



are the decay rate and the CP-conjugated decay 



We see from equations (3.19), (3.22) and (3.41), that the decay density is written as 

^ j dEEpQTDh = ^7^ j dEdk^jNZD \Mhf ■ (4.6) 



iDh 



where 



Zd = {i- Q){i + C - fD = (1 + /r)(i - /D 



(4.7) 



is the statistical factor for the decay, with Bose-enhancement and Fermi-blocking. In the denomina- 



tor of the CP-asymmetry, it is sufficient to take the tree- level matrix element, 
The CP-asymmetry reads 



Mt 



Mt 



^ih{T) 



J dE dk ^ fN Zd {\M, 



\M, 



1 



_k_ 
1 



2 JdEdk^fNZD\Mh\ 



-fh{N ^ LH) 4(27r) 



dEdk—Zo 



Ml 



My, 



(4i 



The CP-asymmetry arises as the interference between tree-level and one-loop diagrams in the 
decay, so we write M = A^o + A^i) where M^ is the tree- level amplitude and M\ the sum of 
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Figure 4.1: The momentum assignments for the vertex contribution to the CP asymmetry. The 
sohd lines without arrows are neutrinos, the ones with arrows the Icptons and the dashed hues the 
Higgs bosons. All momenta are flowing from left to right and P' as indicated. 



all one-loop amplitudes. The matrix elements can be decomposed as Mi = Xili such that the 
CP-conjugated matrix element is M.i = \*Ii- Here, Aj includes the couplings and Ij accounts for 
the kinematics. Thus, 

\M\'^ - \M\'^ = -4 ImAcp Im/cp, (4.9) 

where Xcp = AqA* and Icp = -^o-^i- 



4.2 The Vertex Contribution 

4.2.1 Going to finite temperature 

Calculating the imaginary part of the kinematic term Ivalcp amounts to calculating the imaginary 
part of the one-loop diagram since the tree-level diagram is real. As explained in section 1.4.3, 
there are two one-loop diagrams for the neutrino decay, the vertex diagram and the self-energy 
diagram. The vertex diagram is shown in figure 4.1, along with the momentum assignments. The 
coupling is 

\cP = \o\\ = [{\''\)jk?gsu{2). (4.10) 

where gsu{2) = 2 denotes the sum over the Higgs and lepton doublets, A is the Yukawa coupling 
between neutrino and Higgs and lepton doublet, j = 1 is the decaying neutrino family, k = 2 is 
the family of the neutrino in the loop and we have summed over all fermion spins and the lepton 
families, both external and in the loop. Moreover, 

IV = -^ (2^ [Mk^N'^<j>'{uePRUN)iuNPRSe'PLUi)] * , (4.11) 



where p', q' and k' are the neutrino, Higgs boson and lepton momentum in the loop, Ajv 



(P'^ — M|) ^ is the denominator of the loop neutrino propagator, A^/ accordingly for the loop 
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Higgs, i-*R,L are projection operators, S^i is the loop lepton propagator and uj^ and U£ are the 
external neutrino and lepton spinors. 

The external fermions are thermal quasiparticles and can be written as spinors u^ [74] which 
are eigenstates of (70 =F k • 7) and have modified dispersion relations as explained in section 2.5.2. 
From chapter 3, we know that 



IY.\MUP,K)\' = 9'\^^± iP0Tpri±) , (4.12) 



2^' ^^ ' " -^ 2mi 

s '^ 

where s denotes the spin of the neutrino. We can also write the matrix element as 

^-Y,\MUP,K)\' =^-Y,9\ufPnul,){ul,PLuf). (4.13) 

s s 

From equations (4.12) and (4.13) we derive a rule for multiplying the spinors of the lepton states/ 

uf{K)uf{K) = Z±uj±{jo T k • 7), (4.15) 

where 

is the quasiparticle residuum. 

The HTL lepton propagator is given in equations (2.64) and (3.5) 

5*(i^) = -iA+(K)(7o-k-7)-^A_(K)(7o+k-7), 



2 / , , , , , , N -,-1 



and m£ is the thermal lepton mass. The Higgs boson propagator is 



(417) 



where m^ is the thermal Higgs boson mass. At finite temperature, we sum over the Matsubara 



^For the antiparticlo spinors v, we replace K by ~K and get 

v^KpfiK) = -Z±cu±ijo ± k ■ 7) . (4.14) 
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modes, 

^-i^E' (4-19) 

where 

k'Q = (2n + l)7riT, (4.20) 

since we are integrating over a fermion momentum. 
The spin and hehcity sum are evaluated as 

^{u'^Pru%){u%PrS^: Plu'I) = -Y1 ZhU^hM, Ah>{l - hh'k • k'), (4-21) 

s,h' h' 

where h and h' are the ratios of hehcity over chirahty for the external and the loop lepton. The 
integral reads 

/d'^k' 
^-^MuM.ZhOOh [^N'^^'K'YH^y , (4.22) 

where H- = 1 — /i/i'kk' . 

In order to carry out the Matsubara sum, we use the Saclay-representation for the propagators. 
For the Higgs propagator it is given in equation (2.20), 

K = -iy^ "'"' 2^^[^ + 4(a;,0]e-'^^'^ + n'^(a;,Oe"^'^}, (4.23) 



where UJq/ = \ q''^ + rn?, is the on-shell Higgs energy with the thermal Higgs mass m^ and f^i is 
the Bose-Einstein distribution for the Higgs bosons with energy ooqi. For the lepton propagator the 
Saclay representation is given in equation (3.9), 

A'± = A\h') = - dr' e'^o^' / doj'ph'iu', k')[l - n'^{uj')]e-^'^' , (4.24) 

Jo J-oo 



where 



2m}, 



Ph'ico', k') = :^r^[S{J - co',,) + 6{u;' + coU^)] (4.25) 



is the pole part of the lepton spectral density with the two solutions oj^/ = uj'j. for A-i-(a;^, k')~^ = 
0. We are only interested in the pole part since it corresponds to the two quasiparticle modes. 
Moreover, rrii is the thermal lepton mass and fi' the lepton Fermi-Dirac distribution. The neutrino 
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propagator reads 



An' = - [ dr" ePo-"^{[i _ f^,{ujp,)]e-^p'-" - /^K^p')e"-'""}, (4.26) 



where 



UJr,' 



p'^ + Ml (4.27) 



p' - yy ^ ^"fc 



is the neutrino on-shell energy, which is unaffected by thermal corrections since the couphng to the 
bath is neghgible and /n' is the Fermi-Dirac distribution for the neutrinos. As usual, we can write 
Po = i {2m + 1)ttT as Matsubara frequency and later on continue it analytically to real values of 
Pq. In particular e^o^ = — 1. 

4.2.2 Frequency sums for HTL fermion propagators 

In order to deal with the HTL lepton propagator, we derive frequency sums for the propagator 
parts A-i-(i^) of a fermion propagator. We write the propagator in the Saclay representation as 

AUK) = - f dTe'=°-A^(r,k), 
Jo 

/oo 
doj phfFi-u>)e''^^ . (4.28) 

Since we are only interested in the pole contribution, we write the corresponding spectral density 
as 

^poie ^ -z^[5{uj - ujh) + 5{uj + uj_h)] = -^ Zsh6{u - suJsh) , (4.29) 

s 

where sh in Zgh and ojsh denotes the product of s and h, that is, Zgh = .^+ for s = h = —1 for 
example. We have for the propagator 

/oo 
d^ ^Z,t,5ioj - sa;,^)/p(-a;)e-"^ 
-oo 3 

= -^Z,,/p(-.a;,,)e— ^ = J^ ApJ^ (r, k) , 

s s 

Af^{K) = V ZsHfFi-sUsh) f dr e(^'— )- = V Kl^':{K) , 

Kf;{K) = Z,UF{-suJsh) [ dTe(^°-^'^-)- 

Jo 

= -Zsh-, ^ , (4.30) 

ko - sujsh 
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where Zgh = (w^ — k'^)/{2rn^) is the quasiparticle residuum. 

In dealing with frequency sums of bare thermal propagators, it is very convenient to write 



A,(K) = A_,(-K). 



(4.31) 



As explained in section 2.4, replacing a boson by a fermion amounts to replacing /b(w) by —fp{uj). 
Moreover, calculating a frequency sum of fco times the propagators amounts to replacing ko with 
suj as in 



TY,koAsAK)As,{P - K) = siu:TY,^sAK)As,{P - K), 

fco fco 



(4.32) 



where co = \/k'^ + m^ and m is the mass of the first boson. The same holds for fermions. 

It is straightforward to work out the frequency sums for the resummed lepton propagator 

S2 1 - fpisiUJs^h) + fB{s2Uj) 



T^Al°l''^{ko,uj)As2{po - ko,uj) = Zg.h 



2a; pq - siUJs^h - ^2^ 



(4.33) 



where the other necessary frequency sums can be derived from this by making the appropriate 
substitutions. 



4.2.3 The frequency sum for the vertex contribution 

We calculate the frequency sum of the three propagators in the vertex loop by partial fractioning 



We are using 1^n'{P') = ^N'{-P') and 

Cs4>N' = 

The frequency sum is given by 



Sa,' r 



A P°l'^ A _ ^' A P°''= A 



2UJ^, ''''' '^ 



2a; jv' 



'■s,h' ^s^i 



ko — S(p'io^> + sj^'iuj^/ 



T 2_^ ^^"hJ ^S(l>' ^sN' = Zsh' 






-c, 



s<t)N' 



Z 



shN' 



■^shif)' 



^shN' Nsh<f>' 



where 



ZshN' = 1 - fpisOJsh') - frisN'^N') 
Zsh<f>' = 1 - fpisOJsh') + f{S(f,'^(j,') , 



(4.34) 



(4.35) 



(4.36) 



(4.37) 
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and 

NshN' = QO - SUJsh' — Sn'^N' , 

Nsh4>' =Po- SUJsh' - s^'UJ^' . (4.38) 

Summing over all propagator parts and the helicity-over-chirality ratios, we get 

^ E E Aiv' A^. A'i7_ = Y, J^^ {E-H^ + E+H+} , (4.39) 



rvQ 



where 



C = k • k' . (4.40) 



The coefficients E± are given by 



E^ = F't'^'Af - F^'^'Ai - F^^'A'l + F'''' Af' , (4.41) 

E+ = F^'^'Af - F^'t''A[ - F^'OylO + F'^'^Af, (4.42) 



and the coefficients F^^ read 



Ft>^ = Bl-B^, F^'<^'=Bf-Bt;, 

F''^ = Bl-B^, F'^'=Bl-Bi, 

Ft>''=Bl-Bf,, F^'' = B^'-B%, 

F^'^' = Bl - B^;, fOO = 5° B%. (4.43) 



The factors Bj^^^ and Ai are given by 






where the numerators and denominators read 



Nn =P0-UJ -UJq', 



Nl 


= ko — UJq' 


-UJpf, 


iVO 


= ko+ UJqi 


+ UJp>, 


Nf 


= ko- UJq' 


+ UJp,, 


< 


= ko+ UJq, 


-UJp,, 



]^N =P0 +Uj' +UJq,, 
^N =P0 -Uj' + UJq,, 
Nff = PO + Uj' - UJq, , N^' = ko + UJq, - UJp' , N^' = Qo + Uj' - UJp, , (4.45) 



n; 


= qo-uj - ujp,, 


N', 


= qo +uj' + ujp,, 


< 


= qo - uj' + ujp,, 


<' 


= qo + uj' - ujp, , 
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Figure 4.2: The monicntuni assignments for the self-energy contribution. The sohd hnes without 
arrows are neutrinos, the ones with arrows the leptons and the dashed hnes the Higgs bosons. 



and 






1 - A' + f4>' 
^N = -{fi' + f4>') 



yO 

^N 



z 



N 



fe + 



We can write 



'^h' 



AijJqiljJpi 



Z^ = 1 - fe - In' , 
Z^ = -(1 - fe' - In') , 
Z^ = -{h' - Jn') , 
Zd> = fe' - In' ■ 



F'i'^Af - F^'^'Ai - F't'^'A') + F'''' Af 



(4.46) 



H. 



+ 



F^'-^'^f _ po^'^e _ pN'o^o ^ pOO^N'^ ^^1 



k'g h' h' 

or, more explicitly, 

k'g h' 

= E i^ { [{K - Bn) Af - « - B^) A[ - [bI - 5^) ^0 + (4' - Bi) Af 

h' IP 

Bl - <) A[ - (<' - B%) A? + {Bl - B%) Af 



(4.47) 



+ 



Bf - Bt; ] A'i 



(4.48) 



4.3 The Self-Energy Contribution 



For the self-energy contribution, the integral Is = ly is the same as for the vertex contribution, only 
the momentum relations are different (cf. figure 4.2). The left diagram does not give a contribution 
since the combination of couplings, |(ATA)jfcP, does not have an imaginary part. We can use the 
Saclay representation for the Higgs and the lepton propagator A^/ and A^ as in Eqns. (4.23) and 



58 4. CP-Asymmetries 

(4.24), remembering the different momentum relations. Tlie neutrino propagator simply reads 

Atv' = -^ o, (4.49) 

since the internal neutrino momentum P' is the same as the external neutrino momentum P. 
We can calculate the frequency sum directly using e^"'^ = — 1, 



T Y^ e^o^e^o^' = eP°^5(T' - r) . (4.50) 

and get 

Tj2^^'^'ih') = - J^^dc'p'ih')^ [b§ - B^;) . (4.51) 

Alternatively, we can use equation (4.33) and write 

T-Y^APolc^,/ /^A i 1' \ 7 V 1 - /^(sWsh') + /( V^g') /„ rro^ 

^ Z^^h' sih^^ }^s^,[Po - ko,u}g') = Zsh'- v^^^_ (4 52 

^ ' ^ ^l^q' PO - SUJsh' - S^'UJg' 

fey 

Both calculations lead to 

T^^Y.^^'^'^- = E ^Z,,[{B^ - Bf,)H^ + (< - B%)H+] . (4.53) 



k'g h' h' 



2uj. 



4.4 Imaginary Parts 

The terms B'fr,, and A^ in the vertex contribution correspond to the three vertices where the 
denominator fulfills certain momentum relations when set to zero: the i?Ar-terms correspond to the 
vertex with an incoming Ni and {i' , (f)'} in the loop, the iJ^-terms to the vertex with an outgoing 
(f) and {^"2,^'} in the loop, and the ^^-terms to the vertex with an outgoing i and {N2, 0'} in the 
loop. As an example, the term 

B^ = ^-f^' + U' (454) 

corresponds to the incoming neutrino decaying into the lepton and Higgs boson in the loop. Thus, 
the terms correspond to cuttings through the two loop lines adjacent to the vertex, however, a 
correspondence with the circlings of the RTF [23, 74] is not obvious. Among these cuts, only the 
ones which correspond to a A^i or N2 decaying into a Higgs boson and a lepton are kinematically 
possible at the temperatures where neutrino decay is allowed, that is where Mi < raj). These terms 
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Figure 4.3: The cuts through the vertex contribution at finite temperature. The cuts are closed 
to form circles and the line that denotes the decaying particle in the corresponding 1—^2 process is 
indicated by a blob. 



are B^, Af and Bf 



4.4.1 Regarding the N2 cuts 

The diagrams develop an imaginary part when one of the denominators of the relevant terms Bj^, 
A^ and B^ vanishes. The contributions from these denominators, N^ , N^ and A'^^ , correspond 
to the three possible cuts shown in figure 4.3. As explained in section 1.4.3, the contribution from 
N^ is the only possible cut at zero temperature. At finite temperature, the other two cuts corre- 
spond to exchanging energy with the heat bath. When choosing the imaginary parts corresponding 
to these two cuts, the loop momentum K' is of the order of M2. Since we assume a strong hier- 
archy M2 ^ Af 1 , the thermal factors f^t , /^z and f^t are suppressed by the large loop momentum 
and the contributions become very small. In fact, they turn out to be numerically irrelevant in 
the hierarchical limit. The physical interpretation of this is as follows: Consider for example the 
cut through {i' , N2}, which is given by a vanishing denominator A'^^ . The corresponding thermal 
weighting factor is the numerator Z^ = f^i — f^^ = /^/(l — Jn^) — (1 — fi')fN-2- It corresponds to 
two processes: absorption of a neutrino from the thermal bath and induced emission of a lepton, or 
absorption of a lepton and induced emission of a neutrino. The phase space distribution of the A'^2S 
in the bath is suppressed due to their large mass and also the distribution of £s that have momenta 
large enough to fulfill momentum conservation in the process is suppressed, so the process is sup- 
pressed. Therefore, the thermal factors suppress the contribution from the A'^2-cuts. Only when we 
have degenerate masses M2 > Mi, these cuts will give a contribution similar to the one from N^ . 
In this case^, the energy and temperature scales that correspond to A^i and A'^2 processes are not 



^Note that there is a mass range for M2 where we have a contribution from the A'^2 cuts but no resonant enhance- 
ment by the self-energy contribution, which becomes relevant when AAf = M2 — M\ ^ T. This mass range is at 
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clearly separated and one has to account for the possibility of an asymmetry creation by N2 as well. 
Implications of these cuts were discussed in reference [22]. We do not consider the influence of this 
cuts, since we are in the hierarchical limit, but we present the analytical expression in appendix D. 



4.4.2 Vertex cut through {£', 0'} 

The imaginary part from N^ , which implies cutting through the lepton and Higgs boson in the loop, 
is the only cut that is also possible at zero temperature and the only vertex cut that contributes 
in the hierarchical limit'^. We denote the angle between p and k' with r/'. 



Then 



Im (I ^ dr,'^) =-"" l\ d'?''^(<) = -^ j\ ^'^(^' - ^0) (^-^^^ 



where the angle is 



and 



pk'^ 



"^0 " 2^ ^^^°^' ~ ^™''' ^^'^"^^ 



E^2 = Mf + (w'2 - k'^) - ml- (4.58) 



We get 



Im I T ^ / ^An>A^>A'H^ j =^Im I ^^ E /" dA^'fc'^dy?' T d,^'AN' A^^A'H 

N 

= - TTT^ E fdk'd<p'^Z,,Z^{4 - Af)H^. (4.59) 

It is sufficient to perform the integration over (p' from to vr since cos (p' is the only quantity that 
depends on cp' . 



Ml ~ AM > r 

^The corresponding CP-asymmetry has been calculated in reference [16], but with a thermal factor 1 — fi' + f^> 
2fii f^i instead of the correct 1 — ft' + f^,'. For details, see reference [23]. 
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We note that we can write 

2a;„ 



(fco - UJq>Y - Ul, 



4 - < = 71 " .^2 , .2 ^ 2^P'A]^^ , (4.60) 

where A]^j^ can be viewed as the propagator of the internal neutrino, since we can interpret the 
contribution we are looking at as putting the internal Higgs boson on-shell and thus we have 
ko - u}q' = fco - 90 = P'o- 



4.4.3 Self-energy cut 

For the self-energy diagram, only Nj^ contributes. Taking rj' as the angle between p and k', we get 

I- (^ E / ^ E ^N'A,,A'H^ J =^Im (t^J^^ dk'k'^dr,' f^ d<A'AN.A,,A'H_ J 

^ ^ Y. I ^^'^'t^'-^^'Z^H^. (4.61) 

h' •' ^ 



167r2 M'j - Ml 



Comparing this expression with the contribution from N^ in equation (4.59) , we see that calculating 

.y^ by A^^ = (M| - Ml 



the self-energy contribution amounts to replacing A]y, by A^, = (M? — M|) ^ in the N^-vertex 



contribution. If M^ ^ Mj, we get 

aVN _ A^SN _ 1 (-. fjnX 

SO the self-energy contribution is twice as large as the vertex contribution, e^ « 2e\/, where the 
factor two comes from the fact that we have two possibilities for the components of the SU{2) 
doublets in the loop of the self-energy diagram. This resembles the situation in vacuum. 



4.5 Analytic Expressions for the CP-Asymmetries 

4.5.1 Vertex cut through {£', cp'} 

We simplify the analytic expression for e-yiT) in equation (4.2). For ly in equation (4.22) we get 
'^^'y^NS = ^^E r^^' r<i^'-Z,'Z^{4-Af)H., (4.63) 

N IQ-K^ P ^ Jo Jo ^P' 



62 



4. CP-Asymmetries 



where it is sufficient to integrate (j)' from to vr. Tlie difference of tlie matrix elements reads for 
the vertex contribution 



\M{N ^ 40)1' - \M{N ^ lh4>)\ = - 55C/(2)Im 



AtA 



jk 



MjMk ZhOJh 



47r^ p 
T7 h h ^p' 



(4.64) 



Correspondingly, the difference in decay rates reads 



7(iV ^ 40) - l{N -^ 40) = - tet/(2)Im ■ 



aU 



jfe 



MjMk 
4(27r)5 



fc' 



J] / di^dMfc' / d0'A;F^^Z;, ;^iJ^^/.'(4 - Af )i^- , (4.65) 



where FhN = /iv (1 + /</,)(! ~ /^^) is the statistical factor for the decay. 

We know from chapter 3 that 

Y, \MUN ^ LH)\^ = gsu{2)9c{X^^)nZMPo - hpv) , 

s 

where gc = "^ indicates that we sum over N ^ (pi and N — )■ (pi. Thus 



T{N ^ LhH) = gsui2)9< 



(AtA), 



J3 



W-Kppo 



dkkZDZh{po - hpr]) 



and 



7(iV -^ LhH) = gsu{2)9i 



4(27r)3 



dEdkkfNZoZhipo - hpr]) 



where we have summed over the neutrino degrees of freedom. 



(4.66) 



(4.67) 



(4.68) 



We arrive at 



e/i(r) = - gsu{2) 



lm{[{X^X),k?}MjMk 



k' 



gc{X^X)jj 47r2 



Y, I dEdkdk' I d(f>'kF^^Zh^Z^Zh'{A'^e-Af)H- 

. (4.69) 



7(iV ^ LhN) 4(27r)5 ^ 
Iu,{[{X^XU'} M,M,EhJ dEdkdk' !,- d(P'kF^^lZ,^Z^Z,,{4- Af)K 



10 ■' P^p' 



JdEdkkfNZoZhipo - hpr]) 
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4.5.2 Self-energy cut 

For the self-energy contribution, we get 

MjMk 1 ZhOJ 



^'^k -" ^ hh' -^0 



The difference in decay rates reads 



^ / dfc' / d(t)'k'Zh'Z^H^ . 



7(iV ^ 4</>) - l{N ^ 4^ 



■ 55(7(2) 



Ini 



AtA 



jfc 



MjMk 



^ (M2 - M|) (27r)5 
^ f dEdkdk' I dcp'kF^^Zh—Z^Zh'H^ 



The GP-asymmetry reads 

Im{[(AtA),fc]2} MjMk 1 



eh{T)=-gsu(2)^:^ 



LhN) M] - Ml (27r)5 ^ 



(4.70) 



(4.71) 



^ / d^dfcdA;' / dcp'kF'f^^Zh—Z^Zh'H^ 



P 



Im{[(AtA),fc]2} M,Mk 1 E/.'/di?dM/c'/;d0'/cF^<;^Z,|z;^Z,,i7 



(AtA),,- Mf - Ml 27r2 



J dEdkkf'NZoZhiPQ - hpr]) 



(4.72) 



4.5.3 Symmetry under lepton-mode exchange 

We can use equation (4.7) and collect all factors that depend on k and k', 

(! + /</.- A)(l + U' - MZhZh'kk'A]:,f'H. , 



(4.73) 



where we have suppressed the indices for helicity-over-chirality ratios h and h' . The internal 
neutrino momentum p = k + k' — p' is symmetric under a replacement of k and k' and likewise 
the difference co — ojqi = oj + oj' — po. The Higgs boson momenta q = p — k and q' = p — k' are 
also exchanged when we exchange k and k'. Thus, the CP-asymmetry for the vertex contribution 
is symmetric under an exchange of the internal and the external lepton. This can be understood 
as follows: Taking the imaginary part of MqM\ by putting the internal lepton and Higgs boson 
on-shell corresponds to calculating the product of the amplitudes of two decays and one AL = 2 
scattering with a neutrino in the u-channel, as shown in figure 4.4. It can easily be checked that 
this symmetry also holds for the self-energy diagram, where the corresponding AL = 2 scattering 
has a neutrino in the s-channel. 
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X 



Figure 4.4: The product of diagrams that corresponds to the vertex contribution of the CP asym- 
metry at low temperature. It is symmetric under the exchange of the leptons f. and f.' . 



4.6 The CP- Asymmetry at High Temperature 



At high temperature, where we have the decays of Higgs bosons, the CP-asymmetry on amphtude 
level is defined as 



^^_ |A^(<^^iV4)r-|.M('/'^iV4)| 
^^ ~ \M{4> -^ Nlhf + \M[(I) ^ Nh)\ 



(4.74) 



as in equation (5.76). The external momenta are now related as go = Po + ^o- The momentum 
assignments are shown in figure 4.5. We take q and p as the three-momenta of the initial-state 
Higgs boson and the final-state neutrino as in section 3.3, this way we can directly use the results 
from the CP-asymmetry in neutrino decays. The matrix elements are the same as for the low 
temperature case, so A^(i;^ — ;■ Nih) corresponds to Ai{N — )■ (t)ih)i just the energy relations are 
different and the CP-asymmetry is defined with a minus sign relative to low temperature. The 
self-energy contribution from the external neutrino line is the only CP-asymmetric self-energy, the 
other self energies do not exhibit an imaginary part in the combination of the couplings. The 
couplings read 



Im 



{xtxt} 



55(7(2) 



Im 



AtA 



kj 



-9SU(2) 



Im 



AtA 



jk 



(4.75) 



The integrals for the vertex and the self-energy-contribution are 



Ipt = -TY. I -^M,M,Zf,UhA^,A^,A'f,,H'^''' 



where we remember that A^v' = 1/(M? — M|) for the self-energy graph 



(4.76) 
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Q 





K ^ K 

Figure 4.5: The vertex and the self-energy contribution for the cj) decay. 



The frequency sum for the vertex diagram reads 

fc^ h' 

i:4f:b{[K-<)<-K-'^")^? 



jqi ujp 



01 ^01 



+ 



Bf - B^;,\ A^' 



?o _ u'i>'\ 40 



Bl - 5^ j A\ + \B% - Bj,) Al 
N' TyO\ a£ , fnO rO ^ a<I>' 



B',-B%] A° - (B^' -B%)4 + [Bl - B%) A 



-}■ 

(4.77) 



Since we have M2 ^» Mi, we also have M2 ^ w.^ in the relevant temperature range, so the 
possible contributions are from A''^ , Nf and A''^ ^. Again, the N2 cuts can be neglected because 
they are kinematically suppressed. When taking the discontinuity of the diagrams, we get for the 
angle between p and k'. 



where 



'^'^fi = ^i^Po^'-^4>)^ 



(4.78) 



so we arrive at 



.A' ^,N 



hieh) 



■yh lehj y 



ImAcp 

'4(27r)5 



^^ = m'- ico" - k") - Ml , 



(4.79) 



MjMk V [ dEdkdk' r d(t>'kF^hZh-—Zh'Zf,{Af - Af)H^ , (4.80) 
,, J Jo P^p' 



If nitf, S> M2, we would have contributions from N^ and Nf instead 
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where we can write 

Contrary to the CP- asymmetry in neutrino decays, this expression can not strictly be seen as the 
propagator of the neutrino since the contribution does not correspond to a zero temperature cut but 
is a pure thermal effect induced by the presence of leptons and Higgs bosons in the thermal bath. 
This is illustrated by the factor Zfj- = f^i + fe' = /</,'(! — /^') + (1 + f<p')fi'^^ which describes the 
absorption of a Higgs boson and the stimulated emission of a lepton and the opposite process, the 
absorption of a lepton and the stimulated emission of a Higgs boson. Compared to low temperature, 
we have replaced A^^ Z^ by Aj^fzfj- . 

For the self-energy diagram, the frequency sum is given by 



^E E Ajv' A^,A'F_ = A^,^ ^ Zh> [h^ (< - <) + H+ (i?^ - 5]^)] , (4.82) 

k'g h' ^l' h' 

after taking the discontinuity, the CP-asymmetry reads 

{e!^Hll)s = - '|^tyM^-\^ E/d^dMfc-£d0-fcF,.Z,lfc-Z,,(-Z^')g, , (4.83) 



i k h 

where 



F^, = f;\l-r,;^)il-2r^). (4.84) 

3 have replaced Z^ by Zfj- . The self-energy contrib 
by replacing A^f by A^*^ = (M^ — M^) in the vertex case. For Mj. ^ Mj, we have 



Compared to low temperature, we have replaced Zj^ by Zf^ . The self-energy contribution is given 



A];f«A^f«-^, (4.85) 

so the relation e^ ~ 2€v also holds for the Higgs boson decays. 

Using /^(l — fe) = {f(p + fe)fN, the terms that depend on the lepton momenta k and k' are 

(/<^ + mU + ffj)ZhZh'kk'A^jy'^H. . (4.86) 

where now p' = k + k' -|- p and uj + ooqi = uj + uj' + pq, so the CP-asymmetry in Higgs boson decays 
is symmetric under exchanging the internal and external lepton as well. 



^Reference [16] obtains a different factor f^i — f^i — 2f^ifii due to an incorrect cfioice of cutting rules as explained 
in reference [23]. 
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4.7 One-Mode Approach 



We also calculate the CP-asyninietry within the one-mode approach where we treat the thermal 
mass like a kinematical mass and use lepton propagators ( 
section 3.2. The spin sum corresponding to equation (4.21) then reads 



m() ^ or 



\riini) ^ as in 



S^\u\PRU%){u%PRSt>PLn\) 






2M^I^f,K^K'^, 



iji^i, 



(4.87) 



where A^/ = [k'^ — '^\i)^'^ ■ Iii the frequency sums in equations (4.36) and (4.52), we replace A^|^,^ 
by the usual decomposition A^ £/ in equation (2.49), which means replacing Zf-^i by —s/{2uj') on the 
right-hand sides. One can check that in the final expression for the CP-asymmetry, this amounts 
to replacing the sum of the helicity contributions 



hh' 



ZhZh'{l-hh'0 by ^ 



1 



kk' 

LOkUJk' 



<■ 



(4.88) 



This means that in the two mode treatment, it is forbidden for the external and internal lepton to 
be scattered strictly in the same direction if they have the same helicity or in the opposite direction 
if they have opposite helicity. For the one-mode approximation this is not the case since oJk^k' is 
always larger than kk' . This result illustrates that the leptonic quasiparticles still behave as if they 
are massless in terms of the helicity structure of their interactions, while the one-mode approach 
is not able to describe this behaviour. 

For the CP-asymmetries in the decay densities we get 



(A7^)V ^ [lra{N ^ (t>i) - 7™(iV -, ^.,,y 

dEdkdk' 



m'S 



(At. 






ImAcp 



(27r)5 

ImAcp 

' 2(27r)5 

ImAcp 

■ (27r)5 



h' 



,,FNhZ 



kk'Al^ 



VN 



P 



N' 



MjMk Y^ / dEdkdk' 



^^E^lH^kk'A''^^ 



P 



^N' 



MjMkY f dEdkdk' r dcp'^J^hhikk'A^f 



MjMkV] / dEdkdk' 



b'^J^h^kk'A^t 



K-K' 

LOkCOk' 

K-K' 
K-K' 

UJkiOk' 

K-K' 

LOk^k' 



(4.89) 



where F^h = /Af(l - /a')(1 + U- h), F^h = In^^ - fN)if4, + h) - 

Let us examine the high temperature behaviour of the one-mode approach by calculating the 
CP-asymmetry in the matrix elements of a Higgs boson at rest, where we assume that Mj,m£ <^ 
rriff, <^ Mk- For simplicity, we calculate the self-enegy contribution. The integral that corresponds 
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to equation (4.76) reads 

.■J 

Kg 

The part that contributes to the imaginary part of the diagram is 






W.*« = 2/|^M.M,A.,^<(K.K') 



A' . . , ,A' s l^g' ry<i>' . 



where ^ = {kk')/{kk'), 



We parameterise m^ as m^ = g^T and obtain 

vr 1 — X 
Using the expression 



(4.90) 



lmB% = -7TZ%5iN^ ) = -n^Z^6{^ - 6) , (4.91) 



m^ — k' , , 

^0 = -^^ (4.92) 



and we have neglected Mj and m£. We get 

Im(/o/ntl = -^I^ /" dfe'Z^'(2fe' - m,^) , (4.93) 

where x = M^/M'^ and k = m^/2. For simphcity, we make the approximation 

4' = U + h ^ ^"""'^ + e""^'^ = (1 + e-*=^)e-'='^ , (4.94) 

so 

Im(/o/ntl = -^T^(l+e-'=^) r dk'e~^'P{2k'-m^). (4.95) 

The integral gives 



Im(/o/i*)t: = -^r^(l + e-'=^)r2e-'=^ . (4.96) 

47r 1 — X 



(4.97) 



\Mtot? = 4(AtA)iiK • P = 2{X^X)nglT^ (4.98) 
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we arrive at 

__ A\M\^ ImXcP 1 V^ 1 _, .2,, , „_,,/2^ 8 



IMotP {X^X)uTTl-xgl 
Assuming that g^, <^ 1, we get 



90/2(1 + g-g0/2) ^ ° e~3^/2(;L ^ e-3^/2)gp^ (4^99) 



^.^""^ 32 



eo ff^ 



(4.100) 



Taking (7^ = m^/T ~ 0.42 for T = 10^^ GeV and using the more accurate term in equation (4.99), 
we get e/eo ~ 70, while we get e/eg ~ 90 for equation (4.100). We view this result as a rough 
approximation of the value of the CP-asymmetry in Higgs boson decays at high temperature. 
Both our approximation and the exact numerical solution in the next section give a factor of 100 
difference to the CP-asymmetry in vacuum. 

4.8 Temperature Dependence of the CP- Asymmetries 

We show the temperature dependence of the CP-asymmetries in neutrino decays in the full HTL 
calculation and in the one- mode approach for m^ and ^/2m^ in figure 4.6. We choose Mi = lO^'' GeV 
and normalise the asymmetries by the product of the CP-asymmetry at zero temperature and the 
total decay density in vacuum, eo7o°*- ^^ discussed in sections 4.4.3 and 4.6, the vertex contribution 
and the self-energy contribution have the same temperature dependence for M2 ^ Mi. Moreover, 
as discussed in sections 4.5.3 and 4.6, the asymmetries are the same when we exchange the internal 
and the external lepton, therefore the asymmetry for a plus-mode external lepton combined with a 
minus-mode internal lepton is the same as the asymmetry for a minus-mode external lepton with a 

plus-mode internal lepton, in short, A7^ = A7 y. We see that generally, the thresholds are the 

ones we expect from our analysis of the decay rates in section 3.4. For the one-mode calculations we 
have the expected thresholds at Tq for m^ and at T^ for \/2m^. For all asymmetries where a plus- 
mode lepton is involved, that is A7_|-+, A7^ and A7 y, the phase space is reduced similar to the 

\f2mg case below T^ and an additional reduction of the phase space sets in between T^ and T^ 
since large momenta k or A;' that correspond to a large mass m(fc) become kinematically forbidden. 

Between these two thresholds, T^ and Tq , the asymmetry for A7_| / |_ becomes larger than the 

asymmetry for A7_|__|_. This effect occurs because in the (-|-+)-asymmetry the phase spaces of both 
the internal and the external lepton are suppressed, while for the mixed modes, (H — ) or ( — h), only 
the phase space of one momentum is suppressed, while the phase space of the other momentum 
is still large. The effect is similar to the observation that 7_ becomes larger than 7+ above T^^. 
Relying solely on phase-space arguments, one would expect that A7 /2m, is a good approximation 
for A7-1-+. The fact that A7_|-+ is clearly smaller than A7 72^ is due to two suppressing factors: 
One factor is the effect of the two residues Zfi{k) and Zhi{k'), which suppress the rate somewhat for 
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Figure 4.6: The CP- asymmetries in neutrino decays normalised by the CP-asymmetry in vacuum 
and the total decay density in vacuum, A7/(7Q°*eo). We choose Mi = 10^° GeV and M2 ^ Mi. 
The term A'yhih2 denotes the difference between the decay rate and its CP conjugated rate, which 
is proportional to the CP-asymmetry. Here, hi denotes the mode of the external lepton, while /12 

denotes the mode of the lepton in the loop. For example, A7^ = j{N -^ <l>i+) — jiN ^ 4>i+), 

where a minus- mode lepton is present in the loop. Ajme and A7 /jm denote the rate differences 
for the one-mode approach with a thermal mass me, and an asymptotic thermal mass V2me. 



small momenta k and k' . The other, more important factor is the fact that the helicity structure 
and angular dependence of the integrals are different for the (++)- and the v^'TT-^-case as explained 
in section 4.7. Since neutrino momenta are of the order ~ Mi ~ T for our temperature range, the 
lepton momenta will be of the same order, that is fc > mi, and the leptons and Higgs bosons will 
preferentially be scattered forward. Thus also the angle ^ between the two leptons will be small and 
the factor H^ defined in equation (4.40) is suppressed, while the corresponding one-mode factor 
1 — ^{kk')/{(jjk(jjki) is larger than H^ for small angles and still finite if both leptons are scattered 
strictly in the same direction, that is ^ = 1. We have checked numerically that this is the main 
reason why M^„^^ > A7++ in the range 1/2 Afi <T <T^. 

Since the CP-asymmetries follow the corresponding finite-temperature decay rates that are in 
figure 3.2, it is very instructive to normalise them via these decay rates, that is 7+, 7_, 7^^ and 
7 /2 m ■ This also gives a more intuitive definition of the CP-asymmetries at finite temperature. 
These asymmetries are shown in figure 4.7, normalised by the zero temperature CP-asymmetry. 
Compared to the normalisation via 70 in figure 4.6, we see that the (-|— |-)-asymmetry does not fall 
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Figure 4.7: The CP-asymmetries in neutrino decays normalised by the CP-asymmetry in vac- 
uum and the corresponding total decay density at finite temperature, that is A7+_|_/ (7^*60)) 
A7-+/(7reo), A7+_/(7!°t6o), A7— /(7l°'eo), ^l^J{ll°leo) and M^^Jil'^^^^^), where the 
CP asymmetries A7 are explained in figure 4.6. We choose Mi = lO^'' GcV and M2 ^ Mi. 



as steeply as the corresponding decay rate 7+ between T^ and T^, so the ratio A7_|__(_/7+ is dented 
at the threshold T^ . This illustrates that the (+-|-)-CP-asynLnLetry shows a stronger suppression 
below the threshold r:|^, since it suffers from two phase space reductions and two residues that 
are smaller than one. Therefore, the (-|-+)-asymmetry is not affected as strongly as 7+ by the 
additional suppression above T^ when large momenta k and k' are forbidden and the transition 
over this threshold is smoother than for the decay rate 7+. So 74. falls more steeply than A7_|__|_ 
above the threshold and the ratio of the two rates has a dent at T^ . For the (+— )-asymmetry, 
this effect is even stronger, since it is less suppressed than the (++)-asymmetry above Tj^, so the 
ratio A7_| /(eo7-) rises up to a value of 0(0.1). 

The CP-asymmetries in Higgs boson decays at high temperature are shown in figure 4.8, nor- 
malised to €07*°* and we have assumed that M2 ^ Mi,T. The behaviour is similar to the neutrino 

decays, where the ( )-asymmetry is strongly suppressed, while the (H — )-, the ( — h)- and the 

(-l--l-)-asymmetries have a strict threshold at Tq and are suppressed due to the reduced phase 
space between Tq and T^, as expected. The (H — )- and ( — h)- asymmetries are the same and are 
somewhat less suppressed than the (-|— |-)-asymmetry between the thresholds Tq and T'^. In our 
approximation in the end of section 4.7, we have seen that the difference of the matrix elements 
A|A^P rises as T^, so A7 rises as T^, which can be seen in the plot for all finite-temperature 
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Figure 4.8: The GP-asymmetries in Higgs boson decays normalised by the CP-asymmetry in 
vacuum and the total decay density in vacuum, A'j / {-Jq"^ eo) , where the asymmetries A7 are explained 
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Figure 4.9: The CP-asymmetries in Higgs boson decays normalised by the CP-asymmetry in 
vacuum and the corresponding total decay density at finite temperature, that is A7_(-+/(7^'eo)j 
A7_+/(7^°*eo), A7+_/(7l°'eo), A7__/(7l°'eo), A7™,/(7;„°^'eo) and A7v2™,/(7^^/o), where the 
CP asymmetries A7 are explained in figure 4.6. We choose M\ = 10^° GeV and M^ ^ M\. 
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asymmetries. Tlie one-mode asymmetries A^„i^, Aj/^m ^^^ ^^^ (++)-asymmetry are very close 
to each other at high temperature. 

We also normalise the asymmetries to the corresponding decay rates in figure 4.9 and find that 
they all approach a constant value at high temperature, as is expected since the decay rates also 

rise as T^. The dents in the ratios A7_|_+/(eo7+) and A7_| /(eo7+) with an external plus-mode 

lepton are similar to the ones for the neutrino decays and due to the very strong suppression of 
7+ below the threshold T'^. The numerically dominant asymmetries A'-fmei ^Ij^m ^^*^ ^7++ ^ 
settle at a rather high CP-asymmetry, two orders of magnitude higher than at zero temperature, 
as we expect from our approximate calculation for a Higgs boson at rest at the end of section 4.7. 
This is partly due to a suppression of |A^p which rises as m}, = g'iT'^, but mainly due to the larger 
difference in matrix elements A|A^p for Higgs boson decays. 



CHAPTER 5 



Boltzmann Equations 



Preliminary 

In this chapter, we calculate the Boltzmann equations for leptogenesis. We include decays and 
inverse decays involving neutrinos, leptons and Higgs bosons. Concerning the scatterings, we argue 
that we can neglect all of them expect the on-shell contribution of the AL = 2 scatterings. We take 
into account thermal dispersion relations, but assume distributions close enough to equilibrium that 
we can use Boltzmann equations. For the distribution functions, we use the full quantum statistics, 
that is, Fermi-Dirac and Bose-Einstein statistics, but assume the kinetic equilibrium approximation 
ft = "■j/^'i'^/r'^- It has been shown by reference [105] that this is a good approximation. 

5.1 Particle Kinematics 

The Boltzmann equations describe the time evolution of the distribution function of a particle 
species ^l^. We assume an isotropic and spatially homogeneous universe described by the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) metric [72], 

ds^ = dt^ - a{tf I ^^^ + r^de^ + r2 sin^ M<^2 1 ^ (g^^) 

where a{t) is the cosmic scale factor, which describes the expansion of the universe, k = ±1,0 
specifies the curvature, and {t,r,9,(j)) are the comoving coordinates. 

The trajectory of a particle ip with mass m^ > moving in a gravitational field is given by the 
geodesic equations of motion [113]: 



' + ^ ua Pii Pi> — 

dx^,. 

■■?;■ (5.3) 



^+r(;,p^p^ = o, (5.2) 

dx^i. 



dr 



dr 
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Since s = m^T is the eigen-time of the particle, r is fixed and p^ is the momentum of a particle ip. 

In the FLRW metric the fi = component of Eq (5.2) is given as 

dp L a r. ^ da , ^ 

-^ + -p2=0, With a=— 5.4 

(IT a ^ ot 

Writing p° dp° = |p^| d|p^|, this leads to: 

|p^|a + d|p^| = 

<^|p^| = const. X -. (5.5) 

a 

Therefore the 3-momentum scales as 1/a. 

In general, the Liouville operator describing the evolution of a point particle's phase space in a 
gravitational field is given by 

With this operator the equations of motion (5.2) and (5.3) can be written for the momentum as 

^ = L[p^], (5.7) 



and for the space-time as 

dx^ 



L [x^'] . (5.8) 



dr 
Furthermore, the time derivative of the phase space distribution of a non-interacting gas vanishes, 

df{x,p) 



dr 



0. (5.9) 



Using the equations of motion for the particle we obtain the Boltzmann equations for the non- 
interacting particle species ip, 

L[U{x,p)]=0. (5.10) 

Since we are assuming a Robertson-Walker universe which is isotropic and homogeneous, the 
distribution function /^ depends only on t and |p^|. Therefore, the Boltzmann equation can be 



5.1 Particle Kinematics 77 



written as [72] 



L[U] = E^^-H\P,\'I^ = 0, (5.11) 



where we have omitted arguments for the sake of notational clarity. 

Since p3, = rn'^ and because of the spatial isotropy of the Robertson- Walker-Metric, we have 



,2 9/^ _ „ I I df^ 



C-C/v c)\Pi)\ 

After dividing by E^, equation (5.11) has the form 

Interactions are introduced on the right-hand side by a collision term C[f^], which drives the 
distribution function towards its equilibrium value. The complete Boltzmann equation reads 

i'l/*l = f -i^W^ = C[/,l^ (5^14) 

Thus, the Boltzmann equation in a Robertson- Walker universe has the form of a partial differential 
equation. However, in the radiation dominated phase of the universe, in which leptogenesis takes 
place, equation (5.14) can be written as an ordinary differential equation by transforming to the 
dimensionless coordinates z = m^/T and y^ = |p^|/T. Using the relation dT/dt = —HT, the 
differential operator dt — |p^|ii^3|p.| is written as zHdz, and consequently [114] 

with H (m^) = H W^.^^ . In this form, the Boltzmann equation can be easily solved numerically 
on a grid for specific rescaled momenta y. For the right hand side, we have to sum over the collision 
terms of all processes which involve the particle V ^-iid change the phase space distribution. The 
collision term for a process tp + a + ----'r^i + j + --- is given by [72] ^ 



g^C['>P + a + --- ^i+j + ■■■] = -^^ Yl dpa{27r)^6^{p^ + Pa + ■ ■ ■ 



Pi - Pj 



2 



\M{iP + a + ---^i + j + ---)\'Ufa---{l±mi±f, 



3) 



|2 



\M{i +j + ---^ij + a + ---)\'fify{l±U){l±fa)-- 

(5.16) 



We have chosen a nomiaUsation different from Kolb and Turner, so Cherc = o^Ckt 
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where a = (a,- • • ,i,j,- 



^^~" = J^We-. ^'-'^^ 



The terms (1 ±/j) hold for fermions (— ) and bosons (+) and are interpreted as Fermi-blocking (— ) 
and Bose-enhancement (+). In practice, we will only look at processes which involve three or four 
particles, that is, decays, inverse decays and scatterings. We have included the internal degrees of 
freedom, g^,ga,- ' ' ^9i^9ji' ' ' ^ in the matrix elements, therefore we need to put (7^ in front of the 
collision term since it is not included in the phase-space density /^. 

We integrate equation (5.15) over the phase space of the incoming particle with g^ f d^p^/(27r)'^ 
and arrive at 

"'^ - ^ Yj [7(^ + + >i + J + ---)-l{i + J + >ip + a + ■■■)], (5.18) 



dz H(mj,) 

^ ^' processes 



where 



7(V' + a^ >i+ j ^ ) = -^V / TTTTv? ^[^ + "^ >i + 3 ^ ] 



(2vr)= 



= I W dp/3(27r)'^(5'^(p^ +Pa^ Pi-Pj ) 

X |A^(V + a + • • • ^ i + i + • • • )|' Wa • • • (1 ± /*)(! ± /,) 



(5.19) 



where we now integrate over p^ as well, that is, /3 = (^",0, • • • ,hj,''')- The analogous equation 
holds for 7(i + j + ---^ip + a + ---). 



5.2 Low Temperature 

5.2.1 Neutrino evolution 

The Boltzmann equation for the evolution of the lightest right-handed neutrino is 

+ 7(iVi ^(t)i^) + 7(7Vi ^4)1-)- -/{(pi- -^ Ni)- 7(0^_ -^ Ni)] (5.20) 

where i± denote the two lepton modes. We neglect scatterings since they are of higher order in the 
coupling constant. We will from now on omit the subscript 1 for the neutrino and write N. The 
CP-asymmetry in the matrix element is not relevant for neutrino decay, so we calculate the matrix 



5.2 Low Temperature 79 

element, which is the same for the above processes and define 

\Ml\^ = \M{N -^ HL±)\^ = \M{HL± -^ N)\^ , (5.21) 

where now HL± denotes the sum of leptons and Higgs doublets, i± and 0, and their charge 
conjugated states i± and 0. The matrix elements are, however, different for the different lepton 
modes and also the momentum-conserving delta functions differ from each other. The subscript it 
means (+) or ( — ), not the sum. When summing an expression A± that is dependent on the kind 
of lepton dispersion relation over the lepton modes, we write ^_|_ A±. We have 

For each of the two lepton modes, we have now 

7(iV ^ HL±) - jiHL± ^N)= j dpNdpL±dpH{2T:f5\pN - Ph - Pl±) 

'\MiN ^ HL±)\^ Ml + /h)(1 - fL±) 
- \M{HL± ^ iV)|2 (1 - /jv)/h/l±] . (5.23) 

The term in square brackets in equation (5.23) reduces to 

\Ml\^ [fN{l + /h)(1 - fL±) - (1 - fN)fHfL±] = \Ml\^ [CN^HL± " ChL±^n] , (5.24) 

where 

CN^HL± = /Ar(l + ///)(! - /l±) , 

CHL±^N = (1 - fN)fHfL± ■ (5.25) 

Throughout this chapter, we make the kinetic equilibrium assumption, that is, the phase space 
densities can be written as 



/. = ^/r = ^^fP, (5-26) 



where 



^i — ~eq 1 (5-27) 

where n'^'^ and /'^'^ are the equilibrium number densities and distributions. For the neutrino evo- 
lution, we can assume that the Higgs bosons are in equilibrium since they couple very strongly to 
the thermal bath, fn = f]^- The lepton distributions of the two modes are, strictly speaking, out 
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of equilibrium since leptons and antileptons are created asymmetrically. However, the leptons are 
much closer to equilibrium than the neutrinos, so the neutrino evolution is not influenced by the 
lepton asymmetry. Therefore we approximate the lepton densities with their equilibrium density, 
fL± = /2±- We will relax this assumption in the section on the lepton asymmetry evolution. 
Using the relation 

r^i + /h')(i - fii) = (1 - ri^)fiifZ (5.28) 

we can write 

CN^HL± = XAf/^'^(l + f'j^){l - fl\) = {XN - XNf'f^)fHfl\ , 

CHL±^N = (1 - XNOf^fll ' (5-29) 

and 

CN^HL± - CN^HL± = {xn - l)/|^/2±- (5.30) 

The decay densities are 

7(7V ^ HL±) - ^{HL± ^N)= I dpNdpL±dpH{2nY5\pN - Ph - Pl±) 

= {xN-lhS±, (5.31) 



where 

7^± = J dpNdpL±dpH{27T)*SHpN -PH- Pl±) \MI\^ fj^fH- (5.32) 

Note that 7^_|_ is not the same as the equilibrium decay density in equation (3.42), but differs 
from the latter through the thermal factor fufh- It is an effective decay density, which enters the 
Boltzmann equations. The Boltzmann equation for the neutrinos reads 

where 7^ = 7+ + 7-) or, in analogy to equation (1.29), 

^ = -D^{nM-n'^), (5.34) 

where 

^^ = 4t^ (5.35) 

n^ Hz 
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and we have used H{M]\^) = Hz^ . Most conveniently, the number densities are normahsed by the 
entropy density s in order to factorise their dependence on the expansion of the universe. The 
entropy density scales as 

27r2 

s = g,—T\ (5.36) 

where 5* counts the total number of effectively massless degrees of freedom and is defined in 
equation (1.32). We define all number densities in terms of the entropy density as 

Y. = ^, (5.37) 



s 



then 



The Boltzmann equation reads 

AYn 



Xi = Y^. (5.38) 



,N 



dz sHi 



{xN-im, (5.39) 



where 



^,.^,X.M.).;i^^. ,5.40) 



and the Planck mass is 



Mpi = 1.221 • 10'^ GeV. (5.41) 



5.2.2 Lepton asymmetry evolution 

We set up evolution equations for the two different lepton modes separately and define the phase 
space density of the lepton asymmetry in the respective mode as 

fch = fih - hh- (5-42) 

where /i = ±1 denotes the helicity-over-chirality ratio of the leptons. The final lepton asymmetry 
is then n^° = n^, + n^'i after evaluating the Boltzmann equations for each mode separately. The 
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Figure 5.1: The s-channel contribution to tlie AL — 2 scattering £</) —> 



Boltzmann equations for leptons and antileptons read 
dnehi _ z 






d?^a 



h2 

+ Yl b^^h^^ ^ ^^2</') - 7(4. </> ^ 4i<^)] }, (5.43) 

where we have {hi, /12) = ±1 to account for the second lepton involved in the scatterings and we 
have only included AL = 2 scatterings since the other scatterings involving neutrinos are negligible. 
For the evolution of the lepton asymmetry, we have 

^^ = - -i^{7(4i</' ^ iV) - 7(4i<A ^ iV) - 7(iV ^ 4i</') + 7(iV ^ 41-^) 

+ E ^(^i'^ ^ ^"^^2^) - ^(^"'^1'^ ^ 4^'^) + ^(^^'Z' ^ 4i^) - lihj ^ 4i</')] } . (5.44) 

h2 

At leading order in the couplings, the AL = 2 scatterings are computed at tree level and are 
consequently CP-conserving. However, from these scatterings we must subtract the CP-violating 
contribution where an on-shell A'^i is exchanged in the s channel, shown in figure 5.1. This is 
because in the Boltzmann equations the process is already taken into account by inverse decays 
with successive decays, i(p ^ N ^ i(j) [16, 115, 116]. We must therefore replace the scattering rate 
by the subtracted rate, 

7(4i</' ^ ihj) ^ 7'"^(4i</' ^ 4.0) = 7(4i</' ^ ihj) - 7°""''''^4i0 ^ ihj) (5.45) 

and 7(4i^ ~^ ^h2^) accordingly, where ^o'^-'^hcU is the on-shell contribution. The Boltzmann 



5.2 Low Temperature 83 

equations then read 

^^ = - ;^{7(4i<A ^N)- ^{ih,4> ^N)- 7(iV ^ icp^,,) + 7(iV ^ 4,<^) 
+ E [7'"'(4i'^ ^ ^hj) - f'^^ih,^ ^ 4,.^) 

+ 7^"'^(4.<^ ^ 4,0) - I'^'^ikJ ^ 4i</')] } . (5.46) 

Since the CP-asymmetry in neutrino decays is defined on amphtude level as 

^_ \M{N^cPh)\^-\M{N^^h)\' ,^._, 

^h ~ n ; ^^ ?) V"3-^' j 

|A^(Af^</>4)|' + |A^(A^^04)| 
and |A^(iV -^ (t)ih)? + \M{N -^ 04)^ = \Mohf, we write 

1 I JV 

\MiN ^ 04)1' = \Mi^ih ^N)\^ = —^ \Moh\^ , 

1 _ AT 

\MiN ^ 04)1' = |A^(04 ^ N)\' = -^ \Moh\' ■ (5.48) 

It is useful to write the decay rates for the above 10 2 processes as in section 5.2.1, 

7(process) = / I ldpj(27r) 5 ( / ^Pj] c(process), (5.49) 

where pj denotes the relevant momenta pj\f, pih = Pi^ and p^p = pj and S'^{J2Pj) the momentum 
conservation 5'^{pj\f —Peh ~P(f))- The information about the specific process is encoded in c(process) 
and we have 

CN^e^^ = \M{N^ 04)1' /jv(l - feh){l + U) 

Ci^^^N = |A^(04 ^ N)\^ (1 - fN)fihf^ 

c^^i,-^ = \M{N^ 04)1' fN{l - fih){l + U) 

ci^-^^N = 1-^(04 ^ N)\^ (1 - fN)hhf-4,- (5.50) 

Since we are looking at the lepton asymmetry, the lepton distributions have to be out of equilibrium, 

f _ _ /.eq 

Jmh — ^e/ehJfh ' 



li/ih — -^e/ehJeh 

fcq 

I eh 

llh '^ ^Jih 



JCh — ^C.hJifi 1 

fih + hh^'^fZ (5-51) 



while the Higgs bosons can be assumed to be in equilibrium. 
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As explained in appendix E, the scattering rates can be written as 

= E [^'"'(^/^/'^ ^ 4'^) - 7'"'(4,.^ ^ 4»" = 

= jdpr,dhHAP<p{^^f5\PN-PiK -P^K |-^°j'/,\Vr(l - fN)- (5-52) 
so we define^ 

cr^ = 2e^KfOr(l-/^^) (5.53) 

and calculate the integrand for the right-hand side of the Boltzmann equation (5.46), 

c{N ^ 40) - c(40 ^ ^)-c(iV ^ 40) + c(40 -^N) + 4"b = 

= ^chmr;" + ^^/^") - 2^;:/;,^;'^ (xn - 1) (i - 2/^^) . (5.54) 

We can easily check that this term vanishes when the neutrinos are in equilibrium, xn = 1, and 
there is no previous lepton asymmetry, X£h = 0. 

The Boltzmann equation reads now 

-^ = -Jj^ [-^Ihleh {XN - 1) + — [-iwh + XNlNh)\ , (5-55) 

where 7^^^ = 7^^ is defined in equation (5.32) and 

7^1 = J dpNdpihdp427r)^6\pN -p^- Pm) \Mo\' f^'Qi'^ " 2/^^) 
iNh = j dpNdp(,hdp<j,{2-KY5'^{pN -P4>- Pm) \Mof flh^N^ 
e^H = — I dpNdpihdp^{27:f5\pN - p^ - Pm)e^ \Mo\' ffQ^l - ^f^). (5.56) 

We see that the rates and the CP-asymmetries that enter the Boltzmann equations have slightly 
different thermal factors than the equilibrium rates and CP-asymmetries in equations (3.42) 
and (4.69), which employ the factor /Af(l — /£)(! + f<p) for N decays. 



^Note that our factor d"^ differs from reference [105], where they have the out-of-equihbrium distribution (1 — /jv) 
instead of (1 — /^). However, as derived in appendix E, we must employ /^, even if we had only one lepton mode, 
which also results in a Boltzmann equation for (I — I) which is slightly different from theirs. 
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We may also write 



dYch 



dz sHi . 

or, corresponding to equation (1.30), 



<.7.^ (XN - 1) + ^ (7^. + :.^7;^.)1 (5-57) 



dz 



e^^D^^iuN - n^^) - {WS, + WS^XN)nch, (5.58) 



where 





1 ll 

Hzn''^ 


W^, 


1 -v^ 
-■- iWh 


yyoh — 


Hz 2n2 


V,.,, — 


1 <. 


^ Nh 


Hz2n2- 



(5.59) 



5.3 High temperature 

As discussed in chapter 3, the neutrino processes N -^ i(j) are forbidden when the thermal masses of 
the Higgs bosons and leptons become too large, that is, when tti^ > M^. However, new processes 
with the Higgs as single initial or final state are then allowed, (p -H- Ni. These are the dominant 
contributions to the neutrino and lepton evolution and they can be CP-violating as well, so they 
contribute to generating a lepton asymmetry. We derive the Boltzmann equations for this high 
temperature regime in the following. 



5.3.1 Neutrino evolution 

We derive the Boltzmann equation analogously to section 5.2.1, 



dz Hz 

h 



We have 



jiNLh ^H)- 7(F ^ NLh) = I dpNdpLh^pH{'iT^f5\pH - Pn - PLh) 

X \M{NLh^H)\^fNfLh{l + fH) 

- \M{H ^ NLhf (1 - /jv)(l - fLh)fH] ■ (5.61) 
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The tree-level matrix elements \M^\ are the same at high temperature for the Higgs-processes, 
just the kinematics differ. So we have 

\Ml\^ = \M{NLh ^ H)\' = \M{H ^ NLh)\' • (5.62) 

Again, we assume the Higgs bosons and leptons to be in equilibrium. We write 

lA^ol' [fNfLhil + Ih) - (1 - /7v)(l - fLh)fH] = \Mof [c{N L^ ^H)- c{H ^ NLh)] . (5.63) 

Using the relation 

fN^llil + Ih') = (1 - /^^)(1 - flDa (5.64) 

we get 

ciNLh ^H)- c{H ^ NU) = (XN - 1)(1 - flDfH^- (5.65) 

The Boltzmann equation then reads 

or 

^ = -D^{nN-n'^), (5.68) 

where 7^ = 7^_^+7^_, 

^th = ldpNdpLhdpHi2n)^6\pH-PN-PLh) \MI\' f'^^l - fH) (5.69) 

and 



7l 1 

eq 



D'f' = ^—. (5.70) 
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Figure 5.2: The u-channel contribution to tlie AL — 2 scattering £(j> 



5.3.2 Lepton asymmetry evolution 

The Boltzmann equations for leptons and antileptons read 



d. -mb<''"'' 



7W^4,iV) 



dn 



£hi 



dz 



H{M, 



N 



-{7(4iiV^0)-7(</'^4iA^) 



'n^)] }, 



+ 



^ [7(4i<A ^ 4^0) - 7(4, </) ^ ihM }, 

/l2 



combined we get 
dnchi 



dz 



(5.71) 



(5.72) 



+ 



-^{7(4iiV ^ (^) - 7(4iA^ ^ </>) - 7(<^ ^ 4iiV) + 7(</' ^ 4iiV) 

j; [7(4i0 ^ ihj) - 7(4i<^ ^ 4,'/') + 7(42<^ ^ 4i<^) - lih,^ ^ 4i</')] } . (5.73) 

h2 



At high temperature, there can be no on-shell neutrino in the s-channel of the AL = 2 scatterings, 
but there can be an on-shell neutrino exchange in the ti-channel as shown in figure 5.2. Again, we 
need to subtract the AL = 2 rates since the n-channel on-shell neutrino exchange corresponds to 
a Higgs decay followed by an inverse decay, icj) — )• £N£ — )• (j)i. We replace 



7(4i<^^42 



7r(4,<A^4,0)^7(4,<A 



7r-^"^"(4,</'^4. 



(5.74) 
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and get 

^^ = - ;^{7(4,iV ^^)- ^{h,N ^cP)- 7(<^ ^ ih,N) + 7(0 ^ 4iiV) 



We define tlie CP-asymmetry in Higgs decays as 



^ _ \M{^ ^ iV4)|' - \M{ct> ^ iV4)|' 



|A^(0^iV4)r + |A^(0^iV4) 



2' 



thus 



+ 2^ [7^"^4,<A ^ 4.0) - 7'"^4i0 ^ 4.0) + 7'"^4.0 ^ 4i0) - 7^"^42'/' ^ 4W 

(5.75) 



(5.76) 



\Mi^^Nih)r = \M{ihN 



iMicP ^ iV4)r = \M{ihN ^ 0)1 



1 + 6 



1-e 



\Mh\, 



\Mhr ■ 



(5.77) 



As explained in appendix E, the scattering rates are written as 






,,sub 



7^"^(4>^4, 



^ [7^-^4,0 ^ 40) - 7^"'(4,0 ^ 4, 



1 2 /.cq 



dp^dp,;,^dp427r)^5^(p^-pa^-p^)6^|A^^JV,X(l + V)^A^' (^•'^8) 



SO 



sub 



•^^h V^^h\ hh\^ + /<A )JN 



(5.79) 



and we get for the right-hand side of equation (5.75), 

c(4Ar ^ 0)-c(4iV ^ 0) - c(0 ^ 4iV) + c(0 ^ 4iV) + df 

= xchQUT + xnQ) + 26^(1 - /,T)/r (^^ - 1) (1 - 2/^') 



(5.80) 



The Boltzmann equation reads now 



dz 



ii'z 



a;£fe 



S/^7;?, (a^iv -^) + — [lwh + ^Nllh 



(5.81) 
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where 

A = J dpNdp,,dp427r)''5\pN -p^ + p,,) \Mlf fl\l - f2){l - 2Q) 
^Wh = j dpNdp(>hdp^{2TT)'^5'^{pN -P4> + Pih) \MI\ f^flh 
^Nh = j dpNdpe,hdp^{2T:Y5'^{pN -p<p+Peh) \Mlf fthQ^ 

4 = \ [dpNdpeHdp^{27:)'6\pN-p^+Pih)et |-M°.|'/;"(1 - /,T)(1 - 2/^^). (5.82) 

Analogous to equation (1.30), we may also write 

^ = -4l^Diir^N - nj^) - (< + <.^iv)n^., (5.83) 

where 

C"^ Hz 2n2 



Hz 2n„u 



Normalised by the entropy density, the equation reads 
dYch z 



"S 



'hlth i^N - 1) + ^ {itvh + ^Nlth) ] ■ (5-85) 



dz sHi 

5.4 Interacting Modes 

The Boltzmann equations in the previous sections were derived under the assumption that the 
only interactions in which the leptons take part are the Yukawa interactions with Higgs bosons 
and heavy neutrinos, which have very small coupling constants A. This scenario would imply that 
the two modes only interact with each other via intermediate neutrinos or Higgs bosons, where 
the distributions and also the asymmetries in each mode are to first approximation decoupled. In 
a more realistic model, the lepton modes will couple to each other via the SU{2) and U{\) gauge 
bosons W^ and B^ in processes like i.± — t- Izs^B. While it is conceptually interesting to consider 
the case that the two modes are completely decoupled, it might be more realistic to study the 
scenario where the interactions between the lepton modes are fast enough to keep them in chemical 
equilibrium. 

Chemical equilibrium implies that for species that interact via processes a + h — t- i + j, the 
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corresponding cliemical potentials are related as 

Ma + Mfe = Mi + Mj • (5-86) 

When the processes which create or annihilate the particles and antiparticles of some species are 
fast, for example via a + a ^ i + j, where i and j are in equilibrium and their chemical potentials 
vanish, then the chemical potentials of a and a behave as 

fJ-a + fJ-a = fJ-i + fJ-j = , 

=^ fJ-a = -fJ-a ■ (5.87) 

In order to derive the corresponding Boltzmann equation, we introduce a chemical potential ^/i 
for the lepton mode ih- For simplicity, we approximate the distribution with Maxwell-Boltzmann 
statistics, an approximation which is sufficient to derive the final Boltzmann equations. The dis- 
tribution functions in kinetic equilibrium are 

Mk) = e-'^("'^-^'^) , 
/,-JA;) = e-'^(-'^+^'') , 
fihik) - fihik) = e-^^'^ie^^'^ - e"'^'^'^) « 2(3 ^f,Q , (5.88) 

for Hh <^ T. We assume chemical equilibrium between the plus- and the minus-mode, 

fi+ = 11^ = fii. (5.89) 

Moreover, we can make the approximation that the equilibrium densities are about the same since 
the thermal mass m£ ~ 0.2 T is too small to affect the momentum integration considerably in 



(2^ 
i\^ ~ n^"^ ss nfi , (5.90) 



""=/(W-f"w 



where ti^^q is the distribution for massless leptons. With these approximations, we have 



nc+ 


= 2Pii,nl = 


nc 


— t 


nc± 


= nc+ + nc_ 


= 


2n/:, 


xc± 


— eq 5 







(5.91) 

where the subscript it indicates that we sum over the two modes, contrary to its use in the previous 
sections. We can now add the Boltzmann equations for the two modes in equations (5.57) and (5.85) 



5.5 One-Mode Approximation 91 



and arrive at 



dy£± 



dz sH^ 



A7± {XN - 1) + — ;— i'-fw± + XN7N±) , (5.92) 



where 



Yc± = Yc+ + Yc- , 

A7± = £7+7^+ + ^7-76- , 
7w± = lw+ + 7VK- > 

7Ar± = 1N+ + IN- ■ (5.93) 

The factor 1/4 comes from the fact that xch = xc±/'^- Depending on the temperature regime, we 
either have to employ the Higgs boson or the neutrino rates in the Boltzmann equations. 

5.5 One-Mode Approximation 

As we did in sections 3.2 and 4.7 for the decay rates and the CP asymmetries, we also employ 
the one-mode approach for the Boltzmann equations. The equations are derived in analogy to 
sections 5.2 and 5.3 and read 

dlAT z 

-{XN - InDme, 



dz sHi 

dYc ^ z_ 

dz sHi 



X r 

-A-^rne {xN - 1) + -^ {iWme + XN'yNme) , (5.94) 



where jDmn IWmi-, iNmi and A7m^ are defined in equations (5.32), (5.56), (5.69) and (5.82) and 
one has to make the appropriate replacements for the matrix elements and the lepton dispersion 
relations of the one- mode approach for m^ and \/2m^. 

5.6 Evaluation of the Boltzmann Equations 

We solve the Boltzmann equations for five different scenarios: 

1. the zero temperature case with Maxwell-Boltzmann statistics, 

2. the two-lepton-mode approach where the two modes do not interact with each other, 

3. the two- mode approach where the modes couple strongly to each other, 

4. the one-mode approach for a thermal mass m£, 

5. and the one- mode approach for an asymptotic thermal mass \/2m(. 
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In the decoupled case, the lepton asymmetries for the plus- and the minus-mode evolve separately 
from each other. When solving the equations, one has to specify the initial conditions for the 
neutrino abundance and the lepton asymmetry. We assume a vanishing initial lepton asymmetry 
and distinguish between three cases for the neutrino abundances: 

1. Zero initial abundance: this is the case, for example, when an inflaton field decays mostly 
into SM particles and not into the heavy neutrinos. 

2. Thermal initial abundance: this can be realised when some additional interactions keep the 
neutrinos in equilibrium at T S> Mi, for example via a heavy Z' boson related to ^©(lO) 
unification [117]. 

3. Dominant initial abundance: this is the case, for example, when an inflaton decays predomi- 
nantly into A'^i. 

The coupling (A^^A)ii, which enters the neutrino decay rate, is parameterised by the so-called decay 
parameter K^ defined as 

K^%. (5.95) 

where rh\ and m* are the effective neutrino mass and the so-called equilibrium neutrino mass, 
defined in equations (1-27) and (1.28). The case i^ > 1 is called strong washout regime and the 
case K < 1 weak washout regime. 

We want to analyse the evolution of the neutrino abundance and lepton asymmetries for the 
weak and strong washout regimes and different initial abundances. To this end, we write the 
Boltzmann equations for the different scenarios in the form of equations (5.34), (5.58), (5.68) and 
(5.83), 



^ = -DiY. - y-) 



where 



dYr 

-^ = eoD,(YN - Y^'') -{W + Wnxn)Yc , (5.96) 



i^ = -^, D.- ' ^^ 



W = —^^^^ Wn = —^^^ (5 97) 

One usually refers to D^iY^ — Y^) as source term since this term is responsible for the production 
of a lepton asymmetry. The term {W + Wnxn)Y^ is called washout term since it usually has the 
opposite sign as the source term and reduces the production of the lepton asymmetry. The terms 
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D, D^, W and Wn are different for tlie different scenarios. Note tliat for tlie finite temperature 
cases, D^ is not the same as D and there is an additional washout term Wn due to the quantum 
statistics. Our analysis closely follows the arguments and explanations in reference [116] and the 
interested reader will find a comprehensive explanation of leptogenesis dynamics in the vacuum 
case therein. 

5.6.1 Weak washout for zero initial abundance 

Let us start with the weak washout regime and zero initial abundance. We define a value Zeq by 
the condition 

y^(^cq) = ^^^''(^cq) . (5.98) 

For 2: <C 1, the neutrino abundance is negligible compared to Y^"^ 



N ' 



dz - ^^^ ' 



Dy^^ , (5.99) 



where Y^ is approximately constant for z <^ 1. The entropy density s is proportional to z~^ 
and 7/5 is proportional to z~'^ in vacuum and z^^ for the Higgs boson decays at high temperature 
in the finite temperature cases. Thus D r^ z'^ in the vacuum case and D ~ const, at finite 
temperature. Neglecting y^iti^-i and 2;'"'*'*^^ the integration yields Yat — zD{z)/3 '^ z^ for the 
vacuum case and Yn ~ zD ~ z for the finite temperature cases. We show the numerical results 
for K = 0.005 and zero initial abundance in figures 5.3 and 5.4, where these power laws for Y/v(z) 
can be observed for z < 0.1. We also see that Y^* > Y^ "^ > Y^ » Y^, which reflects 
7d* > 7^ ™'^ > 7^, ^ 70- Between the thresholds z'!^. and z^ , the finite-temperature abundances 
do not evolve much, which reflects that the decay rates are very low or vanishing in this regime. 
The neutrino abundance for the asymptotic mass v2 w-£ does not rise at all at high temperature, 
while YJ^^ rises slightly between Zq and z'\_ , where the rate is non-zero. The two-mode rate 
7-1- is, though very suppressed, present over the whole threshold range between z^ and z^ due to 
the minus-modes, so Y^ rises slightly. At low temperature, z > Zcq, the neutrino abundances of 
the different scenarios are very close to each other, since for z > 2, the rates are very close to the 
vacuum rate, 7^)'^^^ — A'-f'^^^/eo — 70. In this regime, Yj\j is much larger than Y^ since the 
coupling is too small to keep the abundance close to equilibrium. 

Having outlined the evolution of the neutrino abundance at high temperature, we can un- 
derstand the evolution of the lepton asymmetry at high temperature. The washout term in the 
Boltzmann equations proportional to Yc is much smaller than the source term D^{Y]\f — Y^), since 
Yc/€q <^ Y^^, and we have for z <^ Zcq, 



Yj>o ~ -zeoD, 



94 



5. Boltzmann Equations 



0.01 




0.03 



U.J- 2;_|_ Zq Zc Zq Z_|_ J- 



^ N 
^ N 

Yc/^0 



10 



70 



Figure 5.3: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K = 0.005 
and zero initial neutrino abundance. We show the the two-mode cases and the vacuum case. 
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Figure 5.4: Evolution of neutrino density and lepton asymmetry for K = 0.005 and zero initial 
neutrino abundance. We show the one-mode cases and the vacuum case. 
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lf=°--.eo^. (5.100) 

In the high temperature regime, the lepton asymmetry is negative and follows the neutrino abun- 
dance in its absolute value, 

Yc{z) ~ -eo^liv(^) = -—Yn{z) . (5.101) 

D JD 

For the vacuum case, D^ = D and Yc/eo — —Y^, while for the finite temperature cases, D^/D = 
A7/(eo7D) ~ 0(10^).^ Both behaviours can be observed in figures 5.3 and 5.4. The lepton 
asymmetry in the minus-mode obeys the same power law as the other finite-temperature modes, 
but is about a factor 100 lower due to the lower rates. The combined (ib)-abundance closely follows 
the plus- abundance since the influence of the minus-mode rates is very suppressed and also the 
different washout term with factor 1/4 instead of 1/2 can be neglected. 

Before turning to the intermediate temperatures z ~ Zgq; let us discuss the low temperature 
regime. For z > Zcq, Y]\f S> Y^, so the source term dominates and washout can be neglected. Since 
in this regime, D^^^ ~ D ^^ ~ Z) , we can write 



^^^ T^tV V-cqN d^Ar 



eoD{YN - Y^'') = -e— ^ . (5.102) 



To first order, the negative lepton asymmetry created below Zcq and the positive contribution from 
above Zeq have the same magnitude and cancel each other. For the remaining asymmetry that did 
not cancel, the washout contribution up to ^cq and the exact behaviour of the abundances around 
Zeq are crucial. 

Assuming that Y]\f{z = oo) = 0, we get 

y|° ~ eoY^(^eq) - \Yc{Zeq)\ , (5.103) 

SO we see that the evolution of the difference AY(z) = Y/v(z) — 1 1^(^)1 /eg below Zcq is crucial for 
the final lepton asymmetry. For the regime 1 ^ ^; ^ z^q, we can write 

z 1 ( A7\ 

miYf i™ - -j 

For the finite-temperature cases, the CP-asymmetry in the decay rates, A7/eo, is considerably 



"-''• = Hllrffy-'--^)- <^-™> 



^We saw in chapter 4 that A7/(eo7D) ~ O(IO^) at high temperature instead of 0(10^). The discrepancy is 
due to rates and asymmetries that occur in the Boltzmann equations and that are sUghtly different from the usual 
rates and CP-asymmetries in chapters 3 and 4 due to the different statistical factors they employ, for example 1 or 
(1 — 2/^) instead of the usual factor (1 — /^). See the discussion in section 4.5.3 and compare the thermal factors 
in equations (3.43) and (5.69). 
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smaller than the decay rate ^d in the range Zq ^ z <2, which can be seen in figure 4.7. Moreover, 
for the one-mode cases, A7m/eo approaches jm faster than A7+_|- approaches 7+ for the plus-mode. 
Above z^, the ratio A7/(eo7D) for the two one-mode cases is about the same. Thus, the difference 
AD = D — D^ is largest for the two-mode approach, smaller and about the same for the two 
one-mode approaches and vanishing for the vacuum approach, AD^ > AD"^^ ~ /\^£)V2mi ^^^^ 
AD^ = 0. As a result. Ay"'" > AY"^^ ~ ^yv2m^ ^^^ ^^^ ^^^^ therefore the final asymmetries are 
related as Y^ > Y^'^ ~ Y^ "^^ ^ Y^. Note that the final asymmetry is non- vanishing for the 
vacuum case, since the washout at higher temperature is larger than at lower temperature due to 
the larger decay rate. For the finite temperature cases, the difference D — D^ bin the crucial regime 
z ~ Zcq governs the final asymmetry. 

The evolution of the decoupled minus-mode at low temperature is very different but not hard 

to understand. The asymmetry rises at z > z^ , because this is the regime above Zc where A7 |- 

is maximal and therefore D~ is maximal as well. Between the thresholds z'^ and z^ , A7 |- 

is suppressed by the internal plus-lepton and A7__ is suppressed by the residue of the internal 

minus-lepton, so D^ is negligible and Y^ does not rise. Above z ~ 1, A7 ^ falls due to the residue 

of the external minus-mode and Y^ does not change. The final lepton asymmetry therefore does 
not change its sign above z > ^cq and keeps the value it achieves at around 1 < z < 2 where the 
CP-asymmetry A7 |_ becomes small. 

The combined (ib)-mode does not evolve differently from the plus-mode since the infiuence 
from the 7^ rates can be neglected and also the washout term with the additional factor 1/2 is 
not noticeable since washout is very small in all temperature regimes. Summarising, there are 
four differing lepton asymmetries in this regime: the vacuum case, the m^-case, the (+)-case and 
the (— )-case. The v2?7i^-case yields the same asymmetry as the m^-case and the (ib)-case yields 
the same asymmetry as the (-l-)-case. We show the four differing lepton asymmetries together in 
figure 5.5. 



5.6.2 Strong and intermediate washout for zero initial abundance 

For strong washout, the evolution of the neutrino abundance is analogous to the weak washout 
regime, with Y^ ~ z and Y^ ~ z^, as shown in figures 5.6 and 5.7. The couplings are stronger, 
therefore the abundances rise faster and meet Y^ earlier at Zcq ~ 1- For larger z, the couplings are 
strong enough to keep Yj^ close to equilibrium. The evolution of the lepton asymmetry is nicely ex- 
plained in reference [116] for the vacuum case with some rather accurate analytical approximations. 
In this work, we are only interested in the difference of the vacuum case to the finite temperature 
case. In the strong washout regime, the lepton asymmetries rise rather fast and the washout term, 
which competes with the source term, becomes larger than the latter at some temperature Zmin; 
where the lepton asymmetry reaches its most negative value. The source term becomes small when 
Yat approaches its equilibrium value, so the washout term drives the asymmetry evolution back to 
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Figure 5.5: Evolution of neutrino density and lepton asymmetry for K = 0.005 and zero initial 
neutrino abundance. We display the four modes from figures 5.3 and 5.4 that give different final 
lepton asymmetries, that is, the plus-mode the minus-mode, the Tn^-mode and the vacuum case. 
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Figure 5.6: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K = 100 and 
zero initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 5.7: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K 
zero initial neutrino abundance. We show the one-mode cases and the vacuum case. 



100 and 



zero. At z > Zeq, the neutrino abundance slightly overshoots Y^, so the source term changes sign 
and adds to the washout term until the lepton asymmetry becomes positive, where the washout 
term changes its sign as well and is competing again. At low temperature, source term and washout 
term have the same magnitude when the lepton asymmetry reaches a maximum at Zmax- Above 
-Zmaxj the lepton asymmetry is again driven to zero by the larger washout. At very low temperature, 
washout and source term become very small and do not influence the asymmetry further, which 
settles at a final value ^£". We see that in the strong washout regime, the dynamics are governed 
by the washout term. The evolution of the finite temperature lepton asymmetries is analogous 
to the vacuum case, but they settle to a different final value. For the finite temperature cases, 
the equilibrium density of leptons is smaller than in the vacuum case due to the thermal mass 
rrii ~ 0.2 T. The washout term is effectively larger than for the vacuum case and competes with 
the source term in a stronger way. Therefore, the asymmetry evolution appears slightly damped 
compared to the vacuum case and the final asymmetry is marginally lower. The evolution of the 
minus-mode is analogous to the evolution at weak washout, rises fast below z^ and does not change 
above the thresholds since Y^ ~ Y^ in this regime. The combined (ib)-mode tracks the plus-mode 
until washout becomes relevant at z > Zeq- The washout term for the (ib)-mode is always about a 
factor two smaller than for the plus-mode, since we add the minus- and plus-washout rates, where 
the minus-rate is always negligible compared to the plus-rate. Thus, the (it) -abundance is less 
affected by washout, so the dynamics are affected by the source term in a stronger way and the 
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Figure 5.8: Evolution of neutrino abundance Ym{z) and lepton asymmetry Yc{z) for K — 1 and 
zero initial neutrino abundance. We show the two-mode cases and the vacuum case. 



final asymmetry is larger. We can view this behaviour as always distributing half the asymmetry 
in a mode ^_ which couples strongly to £+ and is not affected by washout. The final asymmetry is 
about a factor two larger than for the other scenarios in the strong washout regime. 

The case of intermediate washout is shown in figure 5.8, where we only show the two-mode cases 
since in this regime, the final lepton asymmetries of the one-mode cases are the same as for the 
plus-mode. The dynamics can be viewed as an interpolation between the strong and weak washout 
regimes and the final asymmetries are very similar to each other. 



5.6.3 Non-zero initial abundance 

We also present the dynamics for thermal and dominant initial abundance in figures 5.9-5.14. For 
weak washout and thermal initial abundance, Y/v ^ l^^l for low temperatures, and according 
to equation (5.103), l^" /eo ~ Y^^ . For weak washout and dominant initial abundance, this 
equation holds as well, as can be seen in figures 5.9 and 5.10. For intermediate washout iT ~ 1 
and dominant abundance, shown in figures 5.11 and 5.12, the lepton asymmetry production is 
stopped between the thresholds for the thermal cases and the production above z ~ 1 does not 
succed in producing an asymmetry as high as in the vacuum case. For the (ib)-case, the asymmetry 
production is larger since it is not as much affected by washout. 

For the strong washout regime and large initial neutrino abundances, the dynamics at high 
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Figure 5.9: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K = 0.005 
and tiiermal initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 5.10: Evolution of neutrino abundance Y]s[{z) and lepton asymmetry Yc{z) for K = 0.005 
and dominant initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 5.11: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) ioi K = 1 and 
dominant initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 5.12: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for if = 1 and 
dominant initial neutrino abundance. We show the one-mode cases and the vacuum case. 



102 



5. Boltzmann Equations 



0.01 



0.001 



0.0001 



Yx 10-5 




^ N 

y'N — 
y'cho — -- 

^ N 

nh^ — 

>7Ao 



Figure 5.13: Evolution of neutrino abundance Ym{z) and lepton asymmetry Yc{z) for K = 100 
and tliermal initial neutrino abundance. We show the two-mode cases and the vacuum case. 



temperature are interesting, as shown in figures 5.13 and 5.14, but the interplay between source term 
and washout term at low temperature governs the final asymmetry as in the zero-abundance case. 
We reproduce the well-known fact that the initial conditions do not influence the final asymmetry in 
the strong washout regime, while the arguments concerning the equilibrium distribution of leptons 
with thermal mass and the reduced washout of the (=b)-mode still hold and lead to the same lepton 
asymmetry as for zero initial abundance. The decoupled minus-mode is very much affected by the 
coupling, that is the decay parameter K, and the initial conditions, since the final lepton asymmetry 
is produced at high temperatures. The stronger the coupling, the larger the asymmetry production 
of the minus-mode at high temperatures and the larger the final value. Moreover, the larger the 
initial deviation of the neutrino abundance from equilibrium, the larger the asymmetry production 
and the final asymmetry. The final lepton asymmetry in this mode is thus lowest for neutrinos 
with thermal initial abundance. 



5.6.4 Final lepton asymmetries 

The values of the final asymmetries are shown in figures 5.15-5.17 for different initial abundances. 
For zero initial abundance and weak washout, shown in figure 5.15, the asymmetry for the finite- 
temperature cases are larger than for the vacuum case by about one order of magnitude due to the 
difference of the thermal rates 7/5 and the CP- asymmetries A7 at 2; > z^ . The lepton asymmetry 
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Figure 5.14: Evolution of neutrino abundance Yn{z) and lepton asymmetry Yc{z) for K = 100 
and dominant initial neutrino abundance. We show the two-mode cases and the vacuum case. 
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Figure 5.15: Final value of the lepton asymmetry for different values of K for zero initial neutrino 
abundance. 
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Thermal initial abundance 




K 

Figure 5.16: Final value of the lepton asymmetry for different values of K for thermal initial 
neutrino abundance. 



for the plus-mode is also slightly larger than for the one-mode cases due to a suppression of the CP 
asymmetry compared to the one-mode approaches. For strong washout, the asymmetry production 
in the vacuum case is marginally more efficient than in the thermal cases due to a smaller lepton 
equilibrium distribution, while the lepton asymmetry in the (ib)-approach is by a factor 2 larger 
than in the other cases since half of the asymmetry is in the ^_-modes and not affected by washout. 
The minus-mode is completely decoupled, the lepton asymmetry bears the opposite sign as the 
other lepton asymmetries and rises with stronger couplings, that is with larger decay parameter K. 
As discussed in section 5.4, this scenario might not be realistic since the modes will couple to each 
other via gauge bosons, so an evolution similar to the (ib)-case seems more likely. 

For thermal initial abundance and weak washout, shown in figure 5.16, the final asymmetry 
equals the equilibrium abundance 1^/eo, while in the strong washout regime, it shows the same 
behaviour as in the case of zero initial neutrino abundance. The minus-mode asymmetry is very 
low for thermal initial neutrino abundance since the neutrinos are close to equilibrium at high 
temperatures. Contrary to the zero initial abundance case, it bears the same sign as the lepton 
asymmetries of the other scenarios. 

For dominant initial abundance, shown in figure 5.17, the final asymmetries assume their max- 
imal value in the weak washout regime, when the coupling is weak enough not to wash them out 
at low temperature. For larger couplings ET ~ 1, the thresholds lead to a halted asymmetry pro- 
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Figure 5.17: Final value of the lepton asymmetry for different values of K for thermal initial 
neutrino abundance. 
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Figure 5.18: Final value of the lepton asymmetry for different values of K for zero, thermal and 
dominant initial neutrino abundance. The letters z, t and d denote the curves for zero, thermal and 
dominant abundance. Note that the final asymmetry of the minus-mode has opposite sign for zero 
initial neutrino abundance, compared to the asymmetries of all other cases. 
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duction for the finite temperature cases and not as much asymmetry can be produced as for the 
vacuum case. The (ib)-case shows a larger asymmetry compared to the other thermal cases due to 
the weaker washout. At strong coupling K ^ 1, the asymmetries are the same as for thermal and 
zero initial neutrino abundance. The minus-mode asymmetry is large in all washout regimes, since 
the neutrinos are far from equilibrium at high temperatures when the £_ -asymmetry is produced. 
A summary of the several initial conditions can be seen in figure 5.18, where we have omitted 
the \plmg case since it is very close to the mi case in all scenarios. In the weak washout regime, 
the case with zero initial abundance is most strongly affected by thermal corrections which amount 
to one order of magnitude, and the plus-mode asymmetry is additionally enhanced by a factor 
of about two. In the intermediate regime, the dominant-initial-abundance case is influenced very 
much by thermal corrections. Therefore, a production mechanism for dominant neutrino abundance 
has to take into account such thermal effects. In the strong washout regime, one would naturally 
expect that thermal corrections can be neglected. We see that this is not the case since for strongly 
interacting leptonic quasiparticles, a part of the lepton asymmetry can be hidden in the £_-mode 
which is unaffected by washout, thus producing an asymmetry by up to a factor of two larger than 
at zero temperature. The effect of the thermal lepton mass on the equilibrium distribution of the 
leptons is an interesting feature, but very small and can be neglected for all practical purposes. 



Conclusions 



We see this work as a contribution to the revived quest for calculating and capturing effects that the 
finite density and temperature of the medium has on early universe dynamics. For a minimal and 
self-consistent toy model of leptogenesis, which consists only of neutrinos, leptons and Higgs bosons, 
we have performed an extensive analysis of the effects of HTL corrections. This implies capturing 
the effects of thermal masses, modified dispersion relations and modified helicity structures. We 
put special emphasis on the influence of the two fermionic quasiparticles, which show a different 
behaviour than particles in vacuum, notably through their dispersion relations, but also the helicity 
structure of their interactions. Our work is thus similar to the work done in reference [16], where the 
authors of the latter work did not include the effects of fermionic quasiparticles and get a different 
result for the CP- asymmetries, which are crucial for the evolution of the lepton asymmetry. Our 
toy model produces two lepton asymmetries stored in the two different lepton modes without the 
possibility of an equilibration of these asymmetries by SM processes. Since we expect the lepton 
modes to interact via gauge bosons in the bath, we examine a second case where the modes are 
strongly coupled to each other. As a third and fourth case, we approximate the lepton propagators 
by zero temperature propagators with the zero temperature mass replaced by the thermal lepton 
mass in one case and the asymptotic mass in the other case. We refer to these cases as one-mode 
approach. All four thermal cases are compared to the zero-temperature case. 

Our analysis proceeded in three steps: Calculating interaction rates, calculating the necessary 
CP-asymmetries and deriving and solving the corresponding Boltzmann equations. We calculated 
the decay rates in chapter 3. Thermal decay rates in the context of leptogenesis have been calculated 
using vacuum states as external states and replacing the zero-temperature mass by the respective 
thermal mass [16]. We calculated the decay rates with the full HTL-dependence via the neutrino 
self-energy and the optical theorem. While for scalars, that is the Higgs bosons in our case, this 
gives the same result as inserting thermal masses in the kinematics, the results for fermions are 
somewhat different, both conceptually and numerically. Neglecting the zero-temperature fermion 
mass, the resummation of HTL fermion self-energies results in an effective fermion propagator that 
does not break chiral invariance and is split up in two helicity modes. The external fermion states 
therefore behave conceptually different from the ones with chirality-breaking thermal masses that 
have been inserted in the kinematics by hand. Moreover, one has to take care of one additional 
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mode, which has iniphcations for the Boltzniann equations. 

In order to understand and describe the behaviour of the modes, we have assigned momentum- 
dependent quasi- masses to them in section 2.5.2. Numericahy, the presence of the two modes with 
momentum-dependent masses leads to different thresholds corresponding to the asymptotic masses 
of the two modes, v2 ^n-^ and zero mass, and the zero momentum mass m^. It is a well known fact, 
that for large temperatures, the thermal mass of the Higgs boson becomes so large that neutrino 
decays are no longer allowed, but rather the Higgs bosons decay into neutrinos and leptons. While 
reference [16] finds that decays are possible above or below their respective thresholds, we find 
that the decay into the negative mode is always possible, while the decay into the plus-mode is 
possible up to the thresholds corresponding to the thermal mass m£, but becomes suppressed by 
the reduction of possible lepton momenta at the thresholds with the asymptotic mass v2 w-f- From 
pure kinematical arguments, one would expect the one-mode rate with the asymptotic mass to 
be a much better approximation to the plus-rate than the one-mode rate for the thermal mass. 
The fact that this is not the case shows the importance of the helicity structure of the leptonic 
quasiparticles, which suppresses the rate for the plus-mode such that it is much lower than 7^„ . 
For the Higgs boson decays at high temperature, the thresholds are analogous to the neutrino 
decays at low temperature. We observe the same T^-proportionality for Higgs boson decays as 
reference [16]. 

We calculated the CP-asymmetries in chapter 4. To our knowledge, this is the first correct 
calculation of a CP-asymmetry in leptogenesis that includes HTL-corrections.^ We present rules for 
the product of spinors that are related to the fermionic quasiparticles in equations (4.14) and (4.15) 
and derive frequency sums for the HTL fermion propagator in equation (4.33) . We find four different 
CP-asymmetries corresponding to the four different choices of lepton modes both in the loop and 
as external states. We find the CP-asymmetry to be symmetric under an exchange of the lepton 

mode in the loop and the external lepton mode, such that A7^ = A7 ^. At finite temperature, 

there are three possible cuttings for the vertex contribution, the {i' ,(j)'}-cut that corresponds to 
zero temperature and two additional cuts involving the internal N' , namely through {N',i'} and 
{N',(f)'}, which have been found by [16,21,96] and examined more closely in [22], using the real- 
time formalism. We obtain the same cuts using the imaginary time formalism and concentrate 
on the {£' ,(j)'}-cut, assuming the hierarchical limit of M2 ^ Mi. As expected from the zero- 
temperature result, we find the vertex contribution proportional to the self-energy contribution in 
this limit. Contrary to reference [16], we find that the CP-asymmetry in Higgs boson decays is 
larger than the asymmetry in neutrino decays by about a factor of 100 if appropriately normalised. 



^Reference [16] has an incorrect calculation for the CP-asymmetry of the neutrino and Higgs decays with HTL 
corrections in the propagators and thermal masses in the external states. The discrepancy is discussed in detail in 
reference [23]. The authors of the first reference get a factor 1 ~ ft' + fcp' — '2ft>f<j>> for the neutrino decays and 
f^i — fti — 2f^i fii for the Higgs boson decays due to an erroneous choice of cutting rules in the real time formalism. 
The correct calculation gives 1 — f(i + f^i for the neutrino decays and /^z + fii for the Higgs boson decays. This 
discrepancy is also responsible for our CP-asymmetry in Higgs decays being a factor ten larger than their result. 
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This is due to a suppression of the Higgs boson decay rate and a thermal enhancement of the 
CP-asymmetry by the distribution functions of the Higgs bosons and leptons. We compare the 
CP-asymmetries in the two-mode approach to the CP-asymmetries in the one-mode approach. We 
find that for the two-mode approach, the hehcity structure of the modes prohibits the two leptons 
to be scattered strictly in the same direction while for the one-mode approach, this direction is 
only mildly suppressed. Notably this fact is responsible for suppressing the (++)-CP-asymmetry 
compared to the asymmetries of the one-mode approach, as well as the residues of the plus-modes 
to less extent. 

We derive and evaluate the Boltzmann equations for our approaches in chapter 5, performing 
the crucial subtraction of on-shell intermediate states in appendix E^. We compare the results 
of the Boltzmann equations for our five cases, that is, decoupled lepton modes, strongly coupled 
lepton modes, the one-mode approach with m£, the one-mode approach with ^/2m£, and the 
vacuum case. We assume three different initial values for the abundance of neutrinos: zero, thermal 
and dominant abundance, motivated by different scenarios for the production of heavy neutrinos 
after inflation [16]. In the weak washout regime, we find that using thermal masses enhances 
the final lepton asymmetry by about one order of magnitude for zero initial neutrino abundance. 
This is due to the fact that the CP-asymmetry and the decay rate evolve differently at z > 1 
when using thermal masses, since the CP-asymmetry suffers from an additional suppression by 
thermal masses through the leptons and Higgs bosons in the loop. Due to the helicity structure 
of the modes, the CP-asymmetry of the plus-mode is additionally suppressed, which results in an 
additional enhancement of the plus-mode lepton asymmetry compared to the final asymmetries of 
the one-mode approaches. The enhancement we find is similar to the one found in reference [16] 
in this regime, but hard to compare quantitatively due to their different approach which includes 
scatterings and the discrepancy in the CP-asymmetry. In the strong washout regime, thermal 
masses do not show an influence as expected^. However, when we couple the plus- and minus- mode 
strongly, we observe an enhancement of the lepton asymmetry by a factor of about two, since we 
stored half of the asymmetry in a mode that essentially does not interact with the neutrinos and is 
therefore not affected by washout. For intermediate washout, that is X ~ 1, we find that the lepton 
asymmetries with thermal masses are about one magnitude lower than in the vacuum case when we 
assume dominant initial neutrino abundance. This is due to the fact that the lower CP-asymmetry 
does not succeed in producing as much lepton asymmetry at z > 1 when using thermal masses. A 
decoupled minus-mode would show a behaviour completely different from the other thermal cases 
and the vacuum case for all initial values of the neutrino abundance. The lepton asymmetry in 
such a decoupled mode is produced mainly at high temperature and only slightly affected by the 
development at z > 1, where it decouples from the evolution of the other abundances. Therefore, 



^The thermal factor (1 — /^) we find is difTerent from the factor (1 — /at) that reference [105] uses, where they 
employ the out-of equilibrium distribution function for neutrinos without deriving their result explicitly. 

^There is a slight suppression of the lepton asymmetry for thermal masses, since the thermal mass suppresses the 
equilibrium distribution of the leptons somewhat and thereby enhances the washout term. 
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the washout parameter K, which determines the coupling strength and thereby the asymmetry 
production at high temperatures, is crucial for the final value of the lepton asymmetry in this 
mode, as is the initial neutrino abundance. 

Summarising, we argue that for an accurate description of medium effects on leptogenesis, 
the influence of thermal quasiparticles, notably the effects of the two fermionic modes, cannot be 
neglected. Similar to reference [16], our study shows that thermal masses have a strong effect in the 
weak washout regime, while the effect of fermionic modes has an additional influence on the final 
lepton asymmetry in this regime. We also showed that notably in the strong washout regime, the 
presence of a quasi-sterile lepton mode that is not affected by washout can have a non-negligible 
effect on the final lepton asymmetry. Future studies should clarify the dynamics of the interaction 
between the two fermionic modes and determine whether the evolution of the asymmetries in the 
two modes is closer to the decoupled or the strongly coupled case. 

Another important aspect that might be studied in future works is the influence of the finite 
width of the fermionic modes [118], notably the minus-mode. Such effects could be studied using 
Kadanoff-Baym equations [17,19-21,23,95-100] or some other formalism that takes into account 
non-equilibrium quantum effects. In the quest for a unified description of finite-temperature effects 
on leptogenesis, it is important to include SM interactions in the Kadanoff-Baym studies that are 
under way. To this end, quasiparticle excitations of fermions and gauge-bosons should be taken 
into account. Last but not least, the fermionic modes might have an influence on other related 
dynamics in the early universe that involve fermions, such as thermal production of axions, axinos 
or gravitinos. Notably our introduction of sum rules for quasiparticle spinors in equations (4.14) and 
(4.15), as well as the derivation of frequency sums in equation (4.33) can be a valuable contribution 
to future calculations. 

Future experimental results from the neutrino sector and from high-temperature dynamics, such 
as RHIC or observations at the LHC, will provide interesting input for these theoretical efforts. 
Whatever the outcome of such future experiments and calculations, the quest for a sufficiently 
complete and unified description of the influence of the hot, dense medium on early-universe dy- 
namics, such as leptogenesis or the production of dark matter particles, is guaranteed to remain an 
attractive and vibrant field of research with much insight to be gained both from theoretical and 
experimental side. 
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Green's Functions at Zero Temperature 



We consider the case of a scalar field (p{x) with mass m, where x = {xq,x.), xq = t, we are using 
natural units /i = c = /c^ = 1 and are in the relativistic regime where we use the Minkowski metric 
x^ = x^xf^ = gf_iyx'^x^ = Xg — x^. The amplitude of a particle to travel from x to y if xq < yo 
or from y to x if yo < ^0 is called Feynman propagator and is equal to the two-point correlation 
function defined as 

\^f{x -y)^ {Q\T{(t>{x)(t>{y)m , (A.l) 

where T denotes the time ordering operator which acts on two fields as 



[ 4>{y)4>{x) xo < 2/0 



We can express the scalar field by its Fourier transform 

d^A; 1 



/ (2t)'/.(2.!)V. I°We-"-- + at(l.)e 



(A.3) 



h)=Wk 



where K = (fcojk) is the four-momentum of the field, w^ = vk^ + m^ the energy and the Fourier 
coefficients represent creation a(k) and destruction operators a^(k), which create or destroy a boson 
with momentum k. They act on a state |n(k)) with n particles of momentum k as 



a(k)|n(k)) = V^l"(k)-1) (A.4) 



at(k)|n(k)) = V^(k) + 1 |n(k) + 1) , (A.5) 

in particular a(k)|0) = for all k. 

Using the Fourier representation of (/>, the Feynman propagator can be written as 

J [zTip K'^ — m'^ + \e 
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where the term +i e accounts for the so-called Feynman prescription, that is in the complex fco-plane 
we integrate above the pole at k^ = Uk and below the pole at /cq = —ojk- When we carry out the 
fco-integration, we find for xq > yo 



Af{x -y) = -[ 



d^k 1 



'iK-(x-y) 



(A.7) 



ko=uik 



APPENDIX B 



Analytic Solution of the Dispersion Relations for HTL 

Fermions 



The dispersion relations of the two lepton modes are given by the poles of the corresponding 
propagator. Hence, we seek the zeros of 



D±{K) = A±(i^)-i 



771? / , _ zb/cn — /c , kn + k 



.k,±k+--^[±l-=^^:^\n' 



k \ 2k kf) — k 

The equations D± = can be transformed by the substitutions 

_ kQ + k 
ko — k 
ko — k 1 



(B.l) 



(B.2) 
(B.3) 



ko + k x+ 
c := ^. (B.4) 



mt 



This yields 



k 1 

D±=± (-2c-l + x± -lnx±). (B.5) 

c x± — 1 

s :=-exp(-2c-l) (B.6) 

=f2A: 1 

^± = i^T r ^ ^ [x±+H-s) -lnx±] . (B.7) 

1 + in(— sj x± — 1 

Since the prefactor does not have poles for the values of K we are looking at, solving D± = 
amounts to solving 

x± + In(-s) - In x± = 0, (B.8) 



Further introducing 



leads to 
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which in turn means 

s = -x±e-^±. (B.9) 

This is the defining equation of the Lambert W function [119, 120], thus the solution reads 

x± = -W{s). (B.IO) 

According to the definition in Eq. (B.6) 

-l/e<s<0, (B.ll) 

thus the two real branches of the Lambert function. Wo and W-i, correspond to the two solutions 
we seek. In the range given by Eq. (B.ll) Wq > —1 and W-i < —1. For ko > k we have x+ > 1 
and a;_ < 1. Hence, the physical solutions for x± read 

x+ = -W-iis) and x_ = -Wois). (B.12) 

The corresponding results for uj± are then given by 

, , W---W-1 ,B.i3) 

+ H'_i(s) + 1 ' ' 

Wo{s) + 1 ^ ' 

Making use of the relations [121] 

WoM^) +HWo,-i{z)) = lm, (B.15) 

one can directly prove the result by plugging Eqs. (B.13) and (B.ll) into Eq. (B.l). 
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Quantities at Zero Temperature 



C.l Decay Rate at Zero Temperature 

We know that at T = 0, the decay rate of a heavy neutrino with definite spin at rest is given by 

r,,{N^^£)=^-^M,, (C.l) 

where we have not summed over neutrino spins. Then 



E E levr 

and the decay density for definite neutrino spin is 



r„ = M.r„ = ^i^M, (C,2) 



1 f^ 

jd{N^^£) = —, dEEpf^TD (C.3) 



i^'>^)n ..2 



oo 



Mf / dEpfN (C.4) 



4(27r)3 ^ j^ 
where Ki{z) is a Bessel function of second kind and z = Mj/T. 

C.2 CP- Asymmetry at Zero Temperature 

For the CP-asymmetry at zero temperature, the integrals are the same as for T > 0, just the 
thermal masses are zero and the lepton propagator is 

S'i = |i, (C.6) 
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so the spin sum is 



^Y.^uePRnN){uNPRS'ePLUe) = M.K^'K'^A',, 



(C.7) 



where 






(C.8) 



and we have averaged over the neutrino degrees of freedom. At zero temperature, we can use the 
Cutkosky rules to compute discontinuities, so we replace propagators by (5-functions. There are 
three possible cuttings, however, the cuttings which involve A'2 correspond to a process where one 
massless particle decays into N2 and the other massless particle which is not possible kinematically. 
The only possible cutting is through the Higgs boson and lepton propagator, so 



A' 



-27ri(5(Q'2) 



(C.9) 
(CIO) 



We get 



Im(/o/i*) 



2 

MjMk 
167r 



Disc(/o/i*) 



MjMfc Disc 



1 - (1 + x)ln 



1 + x' 



4 1,/ 



d'^k 



(2vr) 



jA'^A'^A'.K.K'^^ 



(C.ll) 
(C.12) 



where 



Mi 
Mf 



(C.13) 



Analogously, we get for the self-energy 



Im(/o/i*) 



MjMk 1 
IGvr 1 — X 



(C.14) 



The CP-Asymmetry at zero temperature reads 



1 Im[(AtA)2 



jki 



Svr (AtA)jj- 



g{x) 



(C.15) 



where 



g{x) = Vx 



1 , , 1+x 
+ 1- 1 + x In 



1 — X 



(C.16) 
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C.3 Boltzmann Equations at Zero Temperature 

We can derive tiie Boltzmann equations for the neutrino and lepton evolution at zero temperature, 
approximating the phase-space densities with Maxwell-Boltzmann distributions, 

fi{Ei) = eM-ErP) , (C.17) 

where energy conservation in scatterings and decays implies 

fN = IlIh (C.18) 

and there are no Higgs decays at high temperature. For the neutrino evolution, we get, analogous 
to equation (5.39), 

where 



70 = y dpNdpLdpH{2Trf6''{pN -PH- PL) \MqV 17- (C-20) 

The matrix element evaluated at zero temperature reads 

\Mo\^ = Ax2Pn-Pl, (C.21) 

where the factor 4 comes from summing over i and l, as well as over the doublets (e^,0+) and 
(^,'/'°). 

We can express the decay rate 70 in terms of the total decay width rj;^* in the rest-frame of the 
neutrino, 

where g^ = 2 accounts for the internal degrees of freedom of the neutrino, the two spins, and 

T'^\N^HL)= ^^^^^^'^' (C.23) 

describes the decay of a neutrino with a definite spin into {(j)t) and {(j)t). 
Evaluating equation (C.22), we get 

lo = gN^TT\fK^{z), (C.24) 
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where z = M/T and Ki{z) is a Bessel function of second kind. For the equihbrium density of the 
neutrinos, we get 



eq / d?PN 






(C.25) 



so that 



70 _-ntot:^l(£l 



_ tot- 

eq -■- rf 

N 



(C.26) 



We can write the Boltzmann equation as 



Y^ = -D{Ym-Y^% 



(C.27) 



where 






(C.28) 



D = ^:il = zk""'^'^ 



Hin^ 



Koiz) 



(C.29) 



^cq 45 ^AT 2t^ ( \ 



- N 



47r4 5, 



(C.30) 



and 



■^tot 



H\ m* 



(C.31) 



is called decay parameter. 



For the lepton evolution, the subtraction of on-shell propagators can be performed analogously 
to the finite temperature case, so that 



sub/ 



^sub(^^ ^ ^^) _ ^SUb(^^ ^ ^^) ^ ^^^^ 



(C.32) 



where 



eo = 



r(iV -^ i(p) - T{N 

r{N -^ e(p) + r(iv 



(C.33) 
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is calculated in appendix C.2. Analogous to equation (5.57), we get 

'^^'^ ' ' eo(x^-l) + ^)7o, (C.34) 



dz sHi\ ^ ' 2 

which can be rewritten as 



where 



F^ = e^D{YN - Y'^) - WYc, (C.35) 



w^4r^- (C.36) 



i7i2n7' 



We have 



eq _ / d?Pi .eq _ T^ 



^t = 9e TTT-y^fe = 9e^T, (C.37) 



where gi = 2 accounts for the lepton doublet components, so we get 

1 ffJV .3 i 



W = -^^z-'KK2(z). (C.38) 
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APPENDIX D 



The Other Cuts 



D.l Imaginary Parts 

D.1.1 Vertex cut through {N2,(p'} 

We use the conventions for the vertex contribution in A^-decays in chapter 4. For N^ , we shift 
integration variables to dc^q' after carrying out the Matsubara sum over /cq. We consider the angle 



Vq' 



kq' 

kq' 



(D.l) 



between k and q' and write 



I-( //^«']^)=--g^' 



(D.2) 



where the angle is 



'^'"'''^ ^ 2kq' (~'^^^<i' + ^fc) ' 



(D.3) 



and 



Sfc = Mi - {oj' - k') - m, 



(D.4) 



The imaginary part reads 



/d^k' \ 1 / r°° 

— ^A^,A^,A,,F_) , = --gim [tY, dq'q'^drj,, / dcP,,A^,A^,At,,H^ 
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lib E / ^'(^^ir/" [(4' - <) «- + (Bl - BX) H, 



D.1.2 Vertex cut through {A^a,^'} 

For the A'^^ -term, we integrate over k', and choose the polar angle between q and k', 

qk' 



llqk' 



qk' 



We write 



'"^W_/^'"='iv;^ 



-vr 



qk'' 



where the angle is 



and 



^^^'° " 2^ (-2u;qtj' + Sgfc/) , 



S,fc, = M| - (a;'2 - fe'2) - m^. 



The imaginary part is given by 






(D.5) 



(D.6) 



(D.7) 



(D. 



(D.9) 



^AAf/A^/A/^/77_ 



-| / /"OO /"TT 

=— ^Im r ^ / dk'k''^driqk' / d(/)qA:'AAr/A<^/A/j/iJ_ 



(D.IO) 



D.2 Analytic Expressions for the CP- Asymmetries 

D.2.1 Vertex cut through {N2, </>'} 
We get for Nf 

MjMk ZhOj 



Im(/y; 



Nf 



167r2 k 



00 /•TT 



V] / dg'' / d0,/ — Zh 



S^-i?^)i7_ + (B°-<W4 



(D.ll) 



D.2 Analytic Expressions for the CP-Asymmetries 
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The difference in decay rates reads 



AtA 



jk 



J 4(27r)5 



Y, / dEdkdq' / d(l)g,kF'^lZh^Zh 



kcOn 



> - ^N 



Brh - Bat ] H- + { Bj, - Bat H. 



?0 ncl>' 



'(t> ^N 



(D.12) 



and the CP-asymmetry reads 

Im{[(AtAU]2} M,Mfe 



(■h{T) - - gsu{2)' 



-f{N -^ LhN) 4(27r)E 



^ / dEdkdq' / d0,,fcF^^Zft-^Z,,, 



fctJo 



B^) - B%\ H- + [B^- Bj^ ] Hj, 



Im{[(AtA),fc]2}M,Mfe 



5c (AtA),. 



n 



47r2 



E,. /di?dMfc7;d0-fcF^';^z,^z,,[(i3; - Qg, + (i?g - i3^)ff+] 

J dEdkkfNZoZhipo - hpr]) 



(D.13) 



D.2. 2 Vertex cut through {N2,£'} 

For N^ , we get 



Im(/^ 



V)^N 



MjMk ZhUJ 



V/ dk' d(t>,k'—Zh,Z^'{Ai-A'l)H+. (D.14) 



The difference in decay rates reads 

7(iV ^ 4'/') - l{N ^ 4<A) = - ff5l/(2)Im ■ 



AtA 



n2 



jk 



J 4(27r)5 



k' 



V / d^dMfe' / dd^qykF'^^Zh Z^'Zh'{Af - A^DH^ . 



(D.15) 



The CP-asymmetry reads 

Im{[(AtA),fe]2} M.Mfe 



(^h{T) = - gsu{2) 



7(iV ^ LhN) 4(2 



^Y. I'^Edkdk' r dd^^k'kF'i.lZh^Z^' ZyiAf - A^,)H. 



Im{[(At A),,]2} M,Mk T.H' I ^Edkdk' £ d0,,,fcP^<;^Z,^Zf Z,,(Af - A<^,)H^ 



gciX^X)jj 47r2 



JdEdkkfNZDZhivo - hpr]) 



(D.16) 
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where we integrate over (f)qki and we take the coordinate system differently than for N^ . 



APPENDIX E 



Subtraction of On-Shell Propagators 



E.l Low Temperature 

We verify the relation in equation (5.52). The scattering rate ^{Icj) — )■ Icj)) can be split up into four 
scatterings with different kinematics, corresponding to the four possibilities of combining the in- 
and outgoing lepton modes. The scattering rates read 

4x4/ 



7(^/,.(/> -^ ih}(t>) = / (^PthAP't'^PihjP^i'^T^) ^ (Pih, +Pct>- Pihf - Pi) 

X |X(4^0^4^,^)|2/,,,^(l-/,-,p(l + /^), (E.l) 

where {hi,hf) = ±1 denote the helicity-to-chirality ratio of the initial- and final-state leptons (or 
antileptons). We will drop the subscript for this appendix part, unless it is necessary, and all 
equations are valid for one specific mode for each involved lepton, unless otherwise noted. With 
this simplified notation, each of the four matrix elements is evaluated as 



Y^\Miih,<P^ihf^)\^= [(AU)n 



DNr2 



Kpn ■ Pe.hi){PN ■ Pihf) - (PN ■ PN){pih, ■ Pehf] 



(E.2) 



where we sum over the lepton spins Si and sg and the lepton flavours and Dj\f = l/[-P^ — Mfj + 
W%^n{Pn)] is the neutrino propagator in the narrow- width approximation and Tn{p%) the total 
width of the neutrino, which equals the total interaction rate, including both lepton modes. Putting 
the propagator on its mass shell, P^ = M^, we get 



si,si 



|A^°^(4,</'^4,<^)|'= [(AU)n 



I nos|2 , 



2 (pN ■ Pih, ) {PN ■ Pihf ) - Mn {Pth, -Pehf) ' (E • 3) 
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where 



\D' 



N\ 






(E.4) 



In vacuum without thermal masses, this reads 



EK 



(AtA) 



11 






si.si 



<, 



1 + 7?) 



(E.5) 



where the dependence on the angle r] between the external leptons cancels out in the integration 
for symmetry reasons, so we can neglect it and write 



^ \M°\i^ -^ #)f = J^ \M{1^ -^ N)\'^ |D^f \MiN 



(E.6) 



At finite temperature with quasiparticle dispersion relations, we can not derive equation (E.6) 
accurately, but in the narrow-width approximation [16], one assumes that the influence of the 
angle between the external particles is negligible and equation (E.6) holds. 
Using the relations in equation (5.48), we derive 

|X-(4>. ^ ihfh)\' hhJU^ - /ft,)(i + 4/) - \M°%lh,4>^ ^ 4,0/)|' /,-,,4.(i - fihf)ii + Uj) 



1 

4 






|2 



Ml 



2 r 



/£/ii(l j£hf) ~^ hhiJ^hf ^^hJlhM hhf) J<bA^ Jtb.f'^ 



(E.7) 



where we have neglected terms of order e^ and x^- and added the subscripts i and / in the Higgs 
boson distributions to clarify which momentum to use. 



f<t>,i = f<j>{^<j>,i) = f<i>{^N - i^eh,) 



(E.^ 



and likewise for /</,,/. 

For the tree-level, CP-conserving amplitude, we have 

iM'^'^ih^ci) ^ 4,0)1' = iM'^'^ih,^ ^ 4/ 



\M 



AL=2\h^hf ■ 



(E.9) 



For the full amplitude |A^al=2| > the on-shell part is also dominant. Since it is CP-conserving, we 
write 



I A /( 1 2 ,-., I A /(OS 1 2 I 7-iOS|2 '- I • vQ I 2 

\MAL=2\h,hf - \MAL=2\h,hf = \^n\ -^ l-^fej 



Ml, 



(E.IO) 



and we get 



A^tree(^^^^ ^ 4,</>) P /ft J0(1 " hh,){^ + /<?) 



E.l Low Temperature 
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= \MAL=2\h,hf 
Subtracting equations (E.7) and (E.ll), we derive 



JC,hA^ ~ hhj) + hhiJC-hf 



(E.ll) 



f"\th^4> 



7^"'(4> 



dm,dp0dp^-^/<^(2vr) 5 {pun, +P<p-Pihf - P4,) 

2 






'h,\ 



Ml 



e^Hl^liLh^H ^ L,,^H) 



(E.12) 



Using the relations 



(i-/.T)(i+/r) 



(i-/^^)(i-/.T+/r)' 



feq-^eq _ ^cq/-, _ ^cq pcqN 

pcq 



and /-/r.^(l - /;v^) = (1 + /r)(l - Q)fN^ 
which hold for w^v = ^eh + ^ip, it is straightforward to derive 



.sub/ 



(E.13) 
(E.ll) 

(E.15) 



7^"^4,<A ^ 4,<A) - y"'(4,0 ^ 4,0) = 7'"^4,<A ^ 4,0) - f ""(4,0 ^ 4,0), 

7-(L;,^i/ ^ L;,^i/) = j^liLh^H ^ L^^i7) 



(E.16) 



Inserting 1 = / d^pAr/(27r)^(5^(p7v —Pehi — P</>) into equation (E.16), again using the first relation 
from equations (E.12) and the expression for the total neutrino width, as derived in chapter 3, 



rjv(p^) = ^ E I^PLhM^'')'^'(PN-PLh,-PH)\Ml\\l-fll^+f; 



hf=±r 



we arrive at equation (5.52), 

E [7^"'(4,0 ^ 4,0) - 7^""(40 ^ ihfC 



Y, [7^""(4,0 ^ 40) - 7^""(4,0 ^ 4, 



(E.17) 



dpNdpeHAP<t>i'^^rnPN-pm,-p^)e',: \Mix f,lf;^{l - /^^) 
= ef!je^iLf,^H ^ iV). (E.18) 
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E.2 High Temperature 

For the u-channel resonance at high temperature when Higgs bosons decay into neutrinos and 
leptons while the neutrinos are stable, we can derive a relation similar to equation (5.52). The 
width in the on-shell neutrino propagator is then not the decay rate but an interaction rate which 
accounts for the processes where the neutrino interacts with the medium, that is, H — )■ NL and 
NL — 7- H. This width acts as a regulator of the u-channel resonance. 
In the narrow-width approximation, the on-shell amplitude reads 



E l-^°'(4.^ 



E \M{(p ^ iV4,)|' |Z?^f |A^(iV4, 



(E.19) 



Si,Sl 



Sl,Sl 



where the on-shell propagator is the same as in equation (E.4), but the width T^ is given by the 
kinematically allowed processes, H — )■ NL and NL — )■ H. 
Using the relations in equation (5.77), we derive 



|X°^(4>. ^ ihfh)\' fihJW^ - fmf)a + Uj) - \M°'{K.^^ ^ 4,</'/)|' fmj^.il - fihf){l + Uj) 



2 r 



D^r I \Mi\^ Mi^ fcHAi - ni) + QjcK, + Kfiii^ - Q,) m - /: 



1^' 



(E.20) 



Analogous to equation (E.ll), we derive 



A^«-(4> ^ thj4>)\ hhMl - /,-,,)(! + /^) 



\M 



^L=2\h,h 



irif 



fchiC^ ~ fmJ + fmJch 



I eh J J -r Jih^jL.nf 



(E.21) 



so that 



7™^4,</' ^ 4,</>) - 7^"'(4> ^ 4,</') = -etlTc^iLh^H ^ L,^H) 



(E.22) 



Using the relations 



pcqs 



and mi - Q){^ - fn'') = (1 + OQfZ 



(E.23) 

(E.24) 
(E.25) 



which hold for w^ = ojih + ojn, it is straightforward to derive 



7^"^4.'/' ^ ihf<P) - 7'"'^(4> ^ 4,0) = 7'"'^(4,</' ^ 4» - 7'"'^(4,'/' ^ ^, 



(E.26) 



E.2 High Temperature 
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Inserting 1 = / d*^pAf/(27r)^(5'^(p0 —pihi —Pn) into equation (E.26), again using the first relation 
from equations (E.23) and tlie expression for the total neutrino width at high temperature, as 
derived in chapter 3, 



rAf(PAr) = YW ^ / (^PLhfPHi'iTT)'^S^iPH - PLhj 



PN) 



hf=±l 



Ml 



(/S, + /?) 



we arrive at equation (5.78) , 

E [7'"'(4,0 ^ 4,0) - ^^"(^^ ^ ihfC 



Y^ [f "^4,0 ^ 4<^) - f "^4,<^ ^ 4,^ 



(E.27) 



J dpNdpe^.^dp^{27r)''6\pN - Pih. - P^)et \Mlf /.^^l + Z^^^)/^^ . (E.28) 
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